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Preface 



This book is focused on the recent developments on problems of probability 

model under uncertainty by using the notion of nonlinear expectations and, in 
particular, sublinear expectations. Roughly speaking, a nonlinear expectation 
E is a monotone and constant preserving functional defined on a linear space of 
random variables. We are particularly interested in sublinear expectations, i.e., 
E[X + Y]< E[X] + E[Y] for all random variables X, Y and E[XX] = XE[X] if 
A > 0. 

A sublinear expectation E can be represented as the upper expectation of 
a subset of linear expectations {Ee ■ G 0}, i.e., E[JC] = supggQ Ee[X]. In 
most cases, this subset is often treated as an uncertain model of probabilities 
{Pg : 9 G 0} and the notion of sublinear expectation provides a robust way to 
measure a risk loss X. In fact, the sublinear expectation theory provides many 
rich, flexible and elegant tools. 

A remarkable point of view is that we emphasize the term "expectation" 
rather than the well-accepted classical notion "probability" and its non-additive 
counterpart "capacity" . A technical reason is that in general the information 
contained in a nonlinear expectation E will be lost if one consider only its 
corresponding "non-additive probability" or "capacity" F{A) = E[1a]- Philo- 
sophically, the notion of expectation has its direct meaning of "mean" , "av- 
erage" which is not necessary to be derived from the corresponding "relative 
frequency" which is the origin of the probability mcasiirc. For example, when a 
person gets a sample {xi, ■ ■ ■ jXn} from a random variable X, he can directly 
use X = j^^Xi to calculate its mean. In general he uses ip(X) = jj"^ f(xi) 
for the mean of ip{X). We will discuss in detail this issue after the overview of 
our new law of large numbers (LLN) and central limit theorem (CLT). 

A theoretical foundation of the above expectation framework is our new 
LLN and CLT under sublinear expectations. Classical LLN and CLT have been 
widely used in probability theory, statistics, data analysis as well as in many 
practical sitiiations such as financial pricing and risk management. They provide 
a strong and convincing way to explain why in practice normal distributions are 
so widely utilized. But often a serious problem is that, in general, the "i.i.d." 
condition is difficult to be satisfied. In practice, for the most real-time processes 
and data for which the classical trials and samplings become impossible, the 
uncertainty of probabilities and distributions can not be neglected. In fact the 
abuse of normal distributions in finance and many other industrial or commercial 
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domains has been criticized. 

Our new CLT does not need this strong "i.i.d." assumption. Instead of 
fixing a probability measure P, we introduce an uncertain subset of probability 
measures {Pg G &} and consider the corresponding sublinear expectation 
E[X] = supg^Q Ee[X]. Our main assumptions are: 

(i) The distribution of Xi is within a subset of distributions {Fg{x) : 6 G Q} 
with 

Jl = E[Xi]>ix = -E[-Xi]- 

(ii) Any realization of Xi, • • • , X„ does not change the distributional uncer- 
tainty of Xn+l- 

Under E, we call Xi, X2, • • • to be identically distributed if condition (i) is 
satisfied, and we call is independent from Xi, - ■■ ,X„ if condition (ii) is 

fulfilled. Mainly under the above weak "i.i.d." assumptions, wc have proved 
that for each continuous function if with linear growth we have the following 
LLN: 

g 

lim E[<^( — )] = sup ip(v). 

n^oo n tJ.<v<TI 

Namely, the uncertain subset of the distributions of Sn/n is approximately a 
subset of dirac measures {5y : fj, < v < Jl} . 

In particular, if /i = /I = 0, then Sn/n converges in law to 0. In this case, if 
we assume furthermore that = E[Xf] and = —E[—Xf], i = 1,2, - ■ ■ , then 
we have the following generalization of the CLT: 

lim E[^(S„/V^)] =E[^(X)]. 

Here X is called G-normal distributed and denoted by -/V({0} x The 
value E[(/?(X)] can be calculated by defining u{t, x) := E[(p(a; + ViX)] which 
solves the partial differential equation (PDE) dtu = G{uxx) with G{a) := 
i(CT^a"'" — CT^a~). Our results reveal a deep and essential relation between the 
theory of probability and statistics under uncertainty and second order fully 
nonlinear parabolic equations (HJB equations). We have two interesting situa- 
tions: when (f is a. convex function, then 

E[^{X)] = -y== / v{x) exp(- -^)da;, 
V2na^ J-00 ^cr 

but if is a concave function, the above o"^ must be replaced by o;^. liq_ = a = a, 
then A^({0} x [a;^,^^]) — A^(0,(t^) which is a classical normal distribution. 

This result provides a new way to explain a well-known puzzle: many practi- 
tioners, e.g., traders and risk officials in financial markets can widely use normal 
distributions without serious data analysis or even with data inconsistence. In 
many typical situations E[(/5(X)] can be calculated by using normal distribu- 
tions with careful choice of parameters, but it is also a high risk calculation if 
the reasoning behind has not been understood. 
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We call iV({0} x [a^, ct^]) the G-normal distributio n. T his new ty pe of s ublin - 
ear d i strib utions was first introduced in Peng (2006) |l00l | (see also 102 1, 103 1, 



104j . |105| ) for a new type of "G-Brownian motion" and the related calculus of 



Ito's type. The main motivations were uncertainties in statistics, measures of 



risk and superhedging in finance (see El Karoui, Peng and Quenez (1997) 44 1 



Artzner, Delbaen, Eber and Heath (1999) [1], Chen and Epstein (2002) 
Follmer and Schied (2004) [Hlj). Fully nonlinear super-hedging is also a possi- 
ble domain of apphcations (see Avellaneda, Levy and Paras (1995) [5|], Lyons 
(1995) [12, see also Cheridito, Soner, Touzi and Victoir (2007) where a new 
BSDE approach was introduced). 

Technically we introduce a new method to prove our CLT on a sublinear 
expectation space. This proof is short since we have borrowed a deep interior 
estimate of fully nonlinear partial differential equation (PDE) in Krylov (1987) 
[tgI ]. In fact the theory of fully nonlinear parabolic PDE plays an essential 
role in deriving our new results of LLN and CLT. In the classical situation the 
corresponding PDE becomes a heat equation which is often hidden behind its 
heat kernel, i.e., the normal distribution. In this book we use the powerful notion 
of viscosity solutions for our nonlinear PDE initially introduced by Crandall 
and Lions (1983) [i^. This notion is specially useful when the equation is 
degenerate. For reader's convenience, we provide an introductory chapter in 
Appendix C. If readers are only interested in the classical non-degenerate cases, 
the corresponding solutions will become smooth (see the last section of Appendix 

C). ^ ^ 

We define a sublinear expectation on the space of continuous paths from K+ 
to M'' which is an analogue of Wiener's law, by which a G-Brownian motion 
is formulated. Briefiy speaking, a G-Brownian motion {Bt)t>o is a continuous 
process with independent and stationary increments under a given sublinear 
expectation E. 

G-Brownian motion has a very rich and interesting new structure which non- 
trivially generalizes the classical one. We can establish the related stochastic 
calculus, especially Ito's integrals and the related quadratic variation process 
{B). A very interesting new phenomenon of our G-Brownian motion is that its 
quadratic variation process {B) is also a continuous process with independent 
and stationary increments, and thus can be still regarded as a Brownian motion. 
The corresponding G-Ito's formula is obtained. We have also established the 
existence and uniqueness of solutions to stochastic differential equation under 
our stochastic calculus by the same Picard iterations as in the classical situation. 

New norms were introduced in the notion of G-expectation by which the cor- 
responding stochastic calculus becomes significantly more fiexible and powerful. 
Many interesting, attractive and challenging problems are also automatically 
provided within this new framework. 

In this book we adopt a novel method to present our G-Brownian motion 
theory. In the first two chapters as well as the first two sections of Chapter III, 
our sublinear expectations are only assumed to be finitely sub-additive, instead 
of "cr-sub-additive" . This is just because all the related results obtained in this 
part do not need the "cr-sub-additive" assumption, and readers even need not to 
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have the background of classical probability theory. In fact, in the whole part of 
the first five chapters we only use a very basic knowledge of functional analysis 
such as Hahn-Banach Theorem (see Appendix A). A special situation is when 
all the sublinear expectations in this book become linear. In this case this book 
can be still considered as using a new and very simple approach to teach the 
classical Ito's stochastic calculus, since this book does not need the knowledge 
of probability theory. This is an important advantage to use expectation as our 
basic notion. 

The "authentic probabilistic parts", i.e., the pathwise analysis of our G- 
Brownian motion and the corresponding random variables, view as functions of 
G-Brownian path, is presented in Chapter VI. Here just as the classical "P-sure 
analysis" , we introduce "c-sure analysis" for G-capacity c. Readers who are not 
interested in the deep parts of stochastic analysis of G-Brownian motion theory 
do not need to read this chapter. 



This book was based on the author's Lecture Notes |l02l | for several series of 
lectures, for the 2nd Workshop Stochastic Equations and Related Topic Jena, 
July 23-29, 2006; Graduate Courses of Yantai Summer School in Finance, Yan- 
tai University, July 06-21, 2007; Graduate Courses of Wuhan Summer School, 
July 24-26, 2007; Mini-Course of Institute of Applied Mathematics, AMSS, 
April 16-18, 2007; Mini-course in Fudan University, May 2007 and August 2009; 
Graduate Courses of CSFI, Osaka University, May 15- June 13, 2007; Minerva 
Research Foundation Lectures of Columbia University in Fall of 2008; Mini- 
Workshop of G-Brownian motion and G-expectations in Weihai, July 2009, a 
series talks in Department of Applied Mathematics, Hong Kong Polytechnic 
University, November-December, 2009 and an intensive course in WCU Center 
for Financial Engineering, Ajou University. The hospitalities and encourage- 
ments of the above institutions and the enthusiasm of the audiences are the 
main engine to realize this lecture notes. I thank for many comments and sug- 
gestions given during those courses, especially to Li Juan and Hu Mingshang. 
During the preparation of this book, a special reading group was organized with 
members Hu Mingshang, Li Xinpeng, Xu Xiaoming, Lin Yiqing, Su Chen, Wang 
Falei and Yin Yue. They proposed very helpful suggestions for the revision of 
the book. Hu Mingshang and Li Xinpeng have made a great effort for the final 
edition. Their efforts are decisively important to the realization of this book. 
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Chapter I 



Sublinear Expectations and 
Risk Measures 

The sublinear expectation is also called the upper expectation or the upper 
prevision, and this notion is used in situations when the probability models 
have uncertainty. In this chapter, we present the basic notion of sublinear ex- 
pectations and the corresponding sublinear expectation spaces. We give the 
representation theorem of a sublinear expectation and the notions of distribu- 
tions and independence under the framework of sublinear expectations. We 
also introduce a natural Banach norm of a sublinear expectation in order to get 
the completion of a sublinear expectation space which is a Banach space. As 
a fundamentally important example, we introduce the notion of coherent risk 
measures in finance. A large part of notions and results in this chapter will be 
throughout this book. 

§1 Sublinear Expectations and Sublinear Expec- 
tation Spaces 

Let be a given set and let be a linear space of real valued functions defined 
on fl. In this book, we suppose that % satisfies c e W for each constant c and 
\X\ & % \i X & %. The space U can be considered as the space of random 
variables. 

Definition 1.1 A Sublinear expectation E is a functional E : 'H — >■ M satis- 
fying 

(i) Monotonicity: 

E[X]>E[y] ifX>Y. 

(ii) Constant preserving: 

E[c] = c /or c e R. 
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(iii) Sub-additivity: For each X,Y gH, 

E[X + Y] < E[X]+E[Y]. 

(iv) Positive homogeneity: 

E[XX] = \E[X] for A > 0. 

The triple {il,H,E) is called a sublinear expectation space. If (i) and 

(ii) are satisfied, E is called a nonlinear expectation and the triple {^,H,E) 
is called a nonlinear expectation space . 

Definition 1.2 Let Ei and E2 be two nonlinear expectations defined on (fij'H). 
El is said to be dominated by E2 if 

Ei[X]-Ei[Y]<E2[X -Y] forX,YGH. (1.1) 

Remark 1.3 From (iii), a sublinear expectation is dominated by itself. In many 

situations, (iii) is also called the property of self-domination. If the inequality in 

(iii) becomes equality, then E is a linear expectation, i.e., E is a linear functional 
satisfying (i) and (ii). 

Remark 1.4 (iii)+(iv) is called sublinearity. This sublinearity implies 

(v) Convexity: 

E[aX + (1 - a)Y] < aE[X] + (1 - a)E[Y] for a € [0, 1]. 

// a nonlinear expectation E satisfies convexity, we call it a convex expecta- 
tion. 

The properties (ii)+(iii) implies 

(vi) Cash translatability: 

E[X + c]= E[X] +c forcGR. 

In fact, we have 

E[X] +c = E[X] - E[-c] < E[X + c]< E[X] + E[c] = E[X] + c. 

For property (iv), an equivalence form is 

E[XX] = \+E[X] + X-E[-X] for A e K. 

In tliis book, we will systematically study the sublinear expectation spaces. 
In the following chapters, unless otherwise stated, we consider the following 
sublinear expectation space (J7, 'H, E): ii Xi, - ■ ■ , Xn € "H then (p{Xi, • • • , X„) G 
% for each ip S Ci.Ljp(M") where Ci.Lip(ffi") denotes the linear space of functions 
satisfying 

\^{x)-^{y)\<C{l + \xr + \yr)\x-y\ iovx,y&W, 
some C > 0, m e N depending on ^p. 

In this case X = {X\, ■ ■ ■ , Xn) is called an n-dimensional random vector, de- 
noted by X e ■H". 
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Remark 1.5 It is dear that if X e H then \X\, X"^ € H. More generally, 
tp{X)ip{Y) € H if X,Y e U and ip,tp € Q.LipiM.). In particular, if X € H then 
E[|X|"] < 00 for each n e N. 

Here we use C(.Lip(R") in our framework only for some convenience of tech- 
niques. In fact our essential requirement is that contains all constants and, 
moreover, X E H implies \X\ G H. In general, Cj.Lip(R") can be replaced by 
any one of the following spaces of functions defined on M". 

• L°°(M"): the space of bounded Borel- measurable functions; 

• C(,(M"): the space of bounded and continuous functions; 

• C^(M"): the space of bounded and fc-time continuously differentiable func- 
tions with bounded derivatives of all orders less than or equal to k; 

• Cunifi^^)' the space of bounded and uniformly continuous functions; 

• Cb.Lipi^")- the space of bounded and Lipschitz continuous functions; 

• I,°(M"): the space of Borel measurable functions. 

Next we give two examples of sublinear expectations. 

Example 1.6 In a game we select a ball from a box containing W white, B 
black and Y yellow balls. The owner of the box, who is the banker of the game, 

does not tell us the exact numbers of W, B and Y . He or she only informs us 
that W + B + Y = 100 and W = B £ [20, 25]. Let ^ be a random variable defined 



Problem: how to measure a loss X = (^(^) for a given function ip on R. 
We know that the distribution of ^ is 



by 




if we get a white ball; 
if we get a yellow ball; 
if we get a black ball. 




with uncertainty: p /it] = [0.4,0.5]. 



Thus the robust expectation of X = ip{^) is 



E[^{i)] := sup Ep[^{0] 



= sup [^(<^(l) + <^(-l)) + (l-p)<^(0)]. 



pe[£i,M] 



Here, ^ has distribution uncertainty. 
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Example 1.7 A more general situation is that the banker of a gam,e can choose 
among a set of distributions {F{0, ^)}AeB(R),6»ee of a random variable ^ . In this 
situation the robust expectation of a risk position for some ip G Ci,Lip{^) 
is 



E[ip{^)] :=sup / (p{x)F{0,dx). 



Exercise 1.8 Prove that a functional E satisfies sublinearity if and only if it 
satisfies convexity and positive homogeneity. 

Exercise 1.9 Suppose that all elements in % are bounded. Prove that the 

strongest sublinear expectation on H is 

E°°[X] ■.= X* = supX(w). 

wen 

Namely, all other sublinear expectations are dominated by E°° [•] . 



§2 Representation of a Sublinear Expectation 

A sublinear expectation can be expressed as a supremuni of linear expectations. 

Theorem 2.1 Let E be a functional defined on a linear space T-L satisfying sub- 
additivity and positive homogeneity. Then there exists a family of linear func- 
tionals {Eg : 6 G Q} defined on H such that 

E[X] = sup Ee [X] forX&n 

and, for each X gH, there exists Ox & & such that E[X] = Ee^iX]. 

Furthermore, if E is a sublinear expectation, then the corresponding Eg is a 
linear expectation. 

Proof. Let Q = {Eg : 9 € &} be the family of all linear functionals dominated 
by E, i.e., Ee[X] < E[X], for all X €n,Ee€Q. 

We first prove that Q is non empty. For a given X e "H, we set i = {aX : a G 
R} which is a subspace of H. We define / : L -^R by I[aX] = aE[X], Va gM, 
then /[•] forms a linear functional on T-L and / < E on L. Since E[-] is sub- 
additive and positively homogeneous, by Hahn-Banach theorem (see Appendix 
A), there exists a linear functional E on H such that E = I on L and E <E 
on T-L. Thus -E is a linear functional dominated by E such that E[X] = E[X]. 

We now define 

Ee[X] := sup Ee[X] for X eU. 
eee 

It is clear that Ee = E. 

Furthermore, if E is a sublinear expectation, then we have that, for each 
nonnegative element X e Ti, E[X] = -E[-X] > -E[-X] > 0. For each c eR, 
-E[c] = E[-c\ < E[-c] = -c and E[c] < E[c] = c, so we get E[c] = c. Thus E 
is a linear expectation. The proof is complete. □ 
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Remark 2.2 It is important to observe that the above linear expectation Eg is 
only ''finitely additive". A sufficient condition for the a-additivity of Eg is to 
assume that E,[Xi] — > for each sequence {Xi}°^^ of H such that Xi{uj) I 
for each uj. In this case, it is clear that Eg[Xi\ -> 0. Thus we can apply the 
well-known Daniell-Stone Theorem (see Theorem \3.3\ in Appendix B) to find a 
a-additive probability measure Pg on {rt,a{'H)) such that 

Eg[X] = [ X{uj)dPg, X en. 

The corresponding model uncertainty of probabilities is the subset {Pg : 6 £ O}, 
and the corresponding uncertainty of distributions for an n-dimensional random 
vector X inn is {Fx{0,A) := Pg{X £ A) : A € S(M")}. 

In many situation, we may concern the probability uncertainty, and the 
probability maybe only finitely additive. So next we will give another version 
of the above representation theorem. 

Let Vf be the collection of all finitely additive probability measures on 
(ri, we consider Lg°(rj, J^) the collection of risk positions with finite val- 
ues, which consists risk positions X of the form 

N 

X{uj) = J2 ^i^A, (^)' e R, A,£F,i = l,--- ,N. 

i=l 

It is easy to check that, under the norm \\-\\^, 1^0^(^1,7-) is dense in 
For a fixed Q eVf and X £ Lg°(f7, J") we define 

N N . 

Eq[X] = EQ[y2xdAM] := Vx,Q(^,) = / X{uj)Q{duj). 

^=l ^=l •'^ 

Eq : Lg°(ri, J^) — > M is a linear functional. It is easy to check that Eq satisfies 

(i) monotonicity and (ii) constant preserving. It is also continuous under ||X||^. 

\EQ[X]\<snv\X{uj)\^\\X\\^. 

Since Lg° is dense in L°° we then can extend Eq from Lg° to a linear continuous 
functional on h°°{n,T). 

Proposition 2.3 The linear functional Eq[-] :L°°(ri, J^) — > M satisfies (i) and 

(ii) . Inversely each linear functional ri{-) :L°°(r2,J^) — )■ R satisfying (i) and (ii) 
induces a finitely additive probability measure via Q,j(A) = ri{lA), A £ J-. The 
corresponding expectation is rj itself 

r,{X) = / X{Lu)Q,,{du:). 
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Theorem 2.4 A sublinear expectation^, has the following representation: there 
exists a subset Q dVf, such that 

E[X] = sup Eq[X] for X en. 
QeQ 

Proof. By Theorem 12.11 '^e have 

E[X] = supEe[X] for X eV., 
eee 

where Eg is a hnear expectation on H for fixed 9 G Q. 

We can define a new sublinear expectation on h°°{^l,a{H)) by 

Eg[X] ■.^M{Ee[Y];Y>X,Y en}. 

It is not difficult to check that Eg is a sublinear expectation on h°^{fl,a{H)), 
where cr('H) is the smallest a-algebra generated by H. We also have Eg < Eg 
on H, by Hahn-Banach theorem, Eg can be extended from H to ]L°^{n,a{H)), 
by Proposition 12. 31 there exists Q e Vf, such that 

Eg[X] ^Eq[X] ioT X eU. 

So there exists Q CVf, such that 

E[X] = sup Eq[X] for X en. 

□ 

Exercise 2.5 Prove that Eg is a sublinear expectation. 

§3 Distributions, Independence and Product Spaces 

We now give the notion of distributions of random variables under sublinear 
expectations. 

Let X = {Xi, • • • , Xn) be a given n-dimensional random vector on a nonlin- 
ear expectation space {n,H,E). We define a functional on C;.Lip(R") by 

Wxiv] := n^{X)] : ^ e Cl,L^p{R") ^ M. 

The triple (R", C;.Lip(R"), Fx) forms a nonlinear expectation space. ¥x is 
called the distribution of X under E. This notion is very useful for a sublinear 
expectation space E. In this case Fx is also a sublinear expectation. Further- 
more we can prove that (see Remark 12.21) . there exists a family of probability 
measures {i^|-(-)}eee defined on (K",S(R")) such that 

Ex[^] = sup / ip{x)Fx{dx), for each ip e Cb.Lipi^"). 

Thus Fx[-] characterizes the uncertainty of the distributions of X. 



§3 Distributions, Independence and Product Spaces 
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Definition 3.1 Let Xi and X2 be two n -dimensional random vectors defined on 
nonlinear expectation spaces (f2i,?^i,Ei) and (02,^2,1152); respectively. They 

are called identically distributed, denoted by Xi = X2, if 
Ei[(^(Xi)] - Mv{X2)] for ^ e G.L.p(M"). 

It is clear that Xi = X2 if and only if their distributions coincide. We say that 
the distribution of Xi is stronger than that of X2 if¥.i[ip{Xi)] > E,2[(p{X2)], for 
each ip e C;.Lip(lR"). 

Remark 3.2 In the case of sublinear expectations, Xi ^ X2 implies that the 
uncertainty subsets of distributions of Xi and X2 are the same, e.g., in the 
framework of Remark ] 2. 21 

{FxM, ■):9ie 61} = {FxM, ■)-02^ ^2}- 
Similarly if the distribution of Xi is stronger than that of X2, then 

{FxAdi, ■) : ^1 e 61} D {FxM, ■)-d2^ 62}. 

The distribution of X e "H has the following four typical parameters: 

A2:=E[X], ^:=-E[-X], a"^ ■.= ¥.[X% := -¥.[-X\ 

The intervals [/x,//] and [a;^,CT^] characterize the mean- uncertainty and the 
variance-uncertainty of X respectively. 

A natural question is: can we find a family of distribution measures to 
represent the above sublinear distribution of XI The answer is affirmative: 

Lemma 3.3 Let (r2,'H,E) be a sublinear expectation space. Let X G 'H'^ be 
given. Then for each sequence {fn]'^=i C C;.Lip(R'^) satisfying (pn i 0, we have 
E[ip„{X)]iO. 

Proof. For each fixed iV > 0, 

V^{x) < C + ^i(^)I[|.|>^] < C + ^^i^ for each x E M'^^^™, 
where = max|3.|<jv <y5n (a;). We then have 

E[MX)]<k^ + ^nip,{x)\x\% 

It follows from ipn i that k^ ], 0. Thus we have lim„^oo E[(p„(X)] < 
^E[(^i(X)|X|]. Since N can be arbitrarily large, we get E[(/?„(X)] | 0. □ 

Lemma 3.4 Let (f2,'H,E) be a sublinear expectation space and let ¥x[yy] '.= 
E[ip{X)] be the sublinear distribution of X E H'^. Then there exists a family of 
probability measures {Fe}g^Q defined on {R'^,B{M.'^)) such that 

ExM-sup/ ^{x)Fg{dx), G C/,L,p(E'*). (3.2) 
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Proof. By the representation theorem, for the subhnear expectation Fjf[(/j] 
defined on (M'*, C/.Lip(]R")), there exists a family of hnear expectations {/ejeee 
on (K'^,Cz.L»p(K")) such that 

FxM = sup/eM, e C/.L,p(M"). 

By the above lemma, for each sequence {ipn}^=i in Cfc.Lip(M") such that ipn i 
on IR'', Fx[<Pn] i 0, thus feiipn] i for each 6* G 0. It follows from Daniell-Stone 
Theorem (see Theorem 13.31 in Appendix B) that, for each 6* e 8, there exists 
a unique probability measure Fe{-) on {R'^,a{Cb.Ltp{R'^)) = (K'', S(R'^)), such 
that fe[(p] = /jjd if{x)Fg{dx). Thus we have ((X^ . □ 

Remark 3.5 The above lemma tells us that in fact the suhlinear distribution 
Fx of X characterizes the uncertainty of distribution of X which is an subset 
of distributions {Fg}g^Q. 

The following property is very useful in our sublinear expectation theory. 

Proposition 3.6 Let (ri,'H,E) be a sublinear expectation space and X,Y be 
two random variables such that¥J^] = —¥\—Y], i.e., Y has no mean-uncertainty. 
Then we have 

E[X + aY] = E[X] + aE[Y] for a e R. 
In particular, ifE[Y] = E[-Y] = 0, then E[X + aY] = E[X]. 

Proof. We have 

E[aY] a+E[Y] + a"E[-r] = a+E[Y] - a"E[y] = aE[Y] for a£R. 
Thus 

E[X + aY] < E[X] + E[aY] = E[X] + aE[Y] = E[X] - E[-aY] < E[X + aY]. 

□ 

A more general form of the above proposition is: 

Proposition 3.7 We make the same assumptions as the previous proposition. 
LetE be a nonlinear expectation on (^1,1-1) dominated by the sublinear expecta- 
tion E in the sense of IfE[Y] = E[—Y], then we have 

E[aY] ^ aE[Y] aE[Y], aSM (3.3) 

OS well as 

E[X + aY]=E[X]+aE[Y], Xen,aeR. (3.4) 

In particular 

E[X + c]= E[X] + c, for c G M. (3.5) 



§3 Distributions, Independence and Product Spaces 
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Proof. We have 

-]E[r] = E[0] - E[Y] < E[-r] = -E[Y] < -E[Y] 

and 

E[Y] = -E[-Y] < -E[-Y] 
= E[0]-E[-Y]<E[Y]. 

From the above relations we have E[Y] = E[Y] = -E[-Y] and thus Still 
by the domination, 

E[X + aY]-E[X] < E[aY], 

E[X] - E[X + aY] < E[-aY] = -E[aY]. 

Thus (inH) holds. □ 

Definition 3.8 A sequence of n- dimensional random vectors {rii}'^^ defined 
on a sublinear expectation space {yi,'H,E) is said to converge in distribu- 
tion (or converge in law) under E if for each ip e Cb.Lip(R"), the sequence 
{E[ip{r]i)]}°l^ converges. 

The following result is easy to check. 

Proposition 3.9 Let {rji}'^^ converge in law in the above sense. Then the 
mapping ¥[■] : Cb.Lip(IR") — > K defined by 

¥[ip] := lim E[(^(77,)] for ip G Cfc.i,p(R") 

>-oo 

is a sublinear expectation defined on (R", Cb.Lip(]R")). 

The following notion of independence plays a key role in the nonlinear ex- 
pectation theory. 

Definition 3.10 In a nonlinear expectation space (fJj'HjE), a random vector 
Y G "H" is said to be independent from another random vector X G under 
E[-] if for each test function (p G Ci,Lip(MJ^~^"') we have 

E[ip{X,Y)]=E[E[ipix,Y)],^x]- 

Remark 3.11 In particular, for a sublinear expectation space (fij'HjE), Y is 
independent from X means that the uncertainty of distributions {Fy{9, ■) : 9 £ 
0} of Y does not change after the realization of X = x. In other words, the 
"conditional sublinear expectation" ofY with respect to X is E[ip(x^Y)\x=x- In 
the case of linear expectation, this notion of independence is just the classical 
one. 

Remark 3.12 It is important to note that under sublinear expectations the con- 
dition 'Y is independent from X " does not imply automatically that "X is in- 
dependent from Y " . 
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Example 3.13 We consider a case where E is a sublinear expectation and 
X,Y & H are identically distributed with E[X] = E[—X] = and = 
E[X^] > ^ = -E\-X\ We also assume that E[\X\] = E[X+ + X"] > 0, 
thus E[X+] = ^E[\X\ + X] =|E[|X|] > 0. In the case where Y is independent 
from X, we have 

E[XY^] = E[X+ct2 - X-ct"] = (a^ - a'')E[X+] > 0. 

But if X is independent from Y, we have 

E[XY^] = 0. 

The independence property of two random vectors X, Y involves only the 

"joint distribution" of {X,Y). The following result tells us how to construct 
random vectors with given "marginal distributions" and with a specific direction 
of independence. 

Definition 3.14 Let (51^, Hi, Ej), i = 1,2 be two sublinear (resp. nonlinear) 
expectation spaces. We denote 

Hi 0^2 := {Z{coi,ui2) = (p{X{uJi),Y{co2)) : (wi,W2) G Cli x O2, 

{X,Y) eUTxn^, a.iip(M™+")}, 

and, for each random variable of the above form Z{uji,uj2) = (p{X{uji),Y{ui2)), 

(El (g)E2)[Z] := Eii^iX)], where f>{x) ~ E2[ip{x,Y)l x G R". 

It is easy to check that the triple {Qi x Q2,'Hi 0^2, Ei ^^2) forms a sublinear 
(resp. nonlinear) expectation space. We call it the product space of sublinear 
(resp. nonlinear) expectation spaces {flijHi, Ei) and {0,2, 'H2, E2). In this way, 
we can define the product space 

n n n 

i=l i=l i=l 

of given sublinear (resp. nonlinear) expectation spaces (fii, "Hi, E^), i = 
1,2, ••• ,n. In particular, when (f2j,'Hi,Ei) = (fii,?^!, Ei) we have the product 
space of the form (0?,Hf", Ef"). 

Let X, X be two n- dimensional random vectors on a sublinear (resp. non- 
linear) expectation space {fl,'H,E). X is called an independent copy of X if 

X = X and X is independent from X. 
The following property is easy to check. 

Proposition 3.15 Let Xi be an Ui- dimensional random vector on sublinear 
(resp. nonlinear) expectation space {^li,Hi,Ei) for i = I,-- - ,n, respectively. 
We denote 

Yi{uii,--- ,uJn) := Xi{wi), i = !,■■■ ,n. 
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Then Yi, i = I,-- - ,n, are random vectors on (JliLi ®"=i ^"=1 ^«)- 
Moreover we have Yi = Xi and l^i+i is independent from (Yi, • • • ,Yi), for each 
i. 

Furthermore, if (i7i,'Hi,Ei) = (ilij'Hi, Ei) and Xi — Xi, for all i, then we 

also have Yi = Yi. In this case Yi is said to be an independent copy of Yi for 
i = 2, - ■ ■ ,n. 

Remark 3.16 In the above construction the integer n can be also infinite. In 
this case each random variable X € (S'i^i belongs to (JliLi (S'iLi ®i=i ^i) 
for some positive integer k < oo and 

oo k 

(g)E,[X] :=(g)E4X]. 

1=1 i=l 

Remark 3.17 The situation 'Y is independent from X "often appears when Y 
occurs after X , thus a robust expectation should take the information of X into 
account. 

Exercise 3.18 Suppose X,Y eW^ and Y is an independent copy of X . Prove 
that for each a G M, 6 G M'^,a + {b, Y) is an independent copy of a + (6, X) . 

In a sublinear expectation space we have: 

Example 3.19 We consider a situation where two random variables X and Y 
in % are identically distributed and their common distribution is 

Fjf M - FyM = sup / ^iy)F{e,dy) for ip G Q.i,p(M), 
flee Jr 

where for each G 6, ^)}yieZ3(R) is a probability measure on (]R,S(R)). 

In this case, "Y is independent from X " means that the joint distribution of X 
and Y is 



Vx,yH'] = sup / 
Biee JR 



sup / tp{x,y)F{e2,dy) 
e2ee 



F{0i,dx) for^£ C,.L,p(R2). 



Exercise 3.20 Let (ri,H,E) be a sublinear expectation space. Prove that if 
¥\lp{X)] = E[(^(y)] for any ip G Cb,Lip, then it still holds for any Lp G Gi^Lip- 
That is, we can replace ip G Ci^Lip in Definition \3.1\ by ip> Cz Cb,Lip- 



§4 Completion of Sublinear Expectation Spaces 

Let (ri,H,E) be a sublinear expectation space. We have the following useful 
inequalities. 

We first give the following well-known inequalities. 
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Lemma 4.1 For r > and 1 < p,q < oo with p + q = 1; have 

|a + 6r' < max{l,2'-i}(|a|'' + |5r) fora,b£R, (4.6) 

|a5|<^ + ^. (4.7) 

p q 



Proposition 4.2 For each X,Y eH, we have 

n\x + vn < 2'-i(E[|xn + E[|yn), (4.8) 

E[\XY\] < {E[\X\P])^/P ■ (Eilrl'])!/', (4.9) 

{E[\X + Y\P]y/P < {E[\X\P])^/P + {E[\Y\P])^/P, (4.10) 

where r> 1 and 1 < p,q < oo with - + - = 1. 

— ^' ^ p q 

In particular, for I < p < p' , we have (E[|X|p])1/p < {E[\X\p']Y'p' . 



Proof. The inequality (|4.8p follows from (|4.6 
For the case E[\X\p] ■ E[|r|«] > 0, we set 



X Y 



(E[|x|p])i/p' (E[|r|9])i/9- 



By (li?7)) we have 



W hi' l£|P hi' 

p q p q 

1 1 
= - + - = 1. 

p q 

Thus gH) follows. 

For the case E[|X|p] • E[|y|'] = 0, we consider E[|X|p] +£ and E[|y|'] +£ for 
e > 0. Applying the above method and letting £ — >■ 0, we get (|4.9p . 
We now prove (|4.10p . We only consider the case E[|X + Y\p] > 0. 

E[\X + Y\P] = E[\X + Y\-\X + Y\P-^] 

< E[\X\ ■ \X + Y\P-^] + E[\Y\ ■ \X + Y\P-^] 

< {E[\X\P]y/P ■ iE[\X + y|(P-i)9])i/9 

+ {E[\Y\p])^/p ■ {E[\x + r|(p-i)'])i/'. 

Since {p — l)q = p, we have (|4.10p . 

BygH), it is easy to deduce that (E[|X|p])1/p < (E[|X|p'])1/p' for 1 < p < p' . 

□ 

For each fixed p > 1, we observe that = {X e H, E[\X\p] = 0} is a 
linear subspace of H. Taking as our null space, we introduce the quotient 
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space T-L/Hq. Observing that, for every {X} e 'H/'Hq with a representation 
X G T-l, we can define an expectation E[{X}] :— E[X] which is still a sublinear 
expectation. We set \\X\\p :— (E[|X|p]) p . By Proposition S^l it is easy to check 
that 11 'lip forms a Banach norm on H/Hq. We extend "H/Hq to its completion 
Hp under this norm, then {Hp, \\-\\p) is a Banach space. In particular, when 

p 1, we denote it by ("H, ||-||). 
For each X G H, the mappings 

X+iu):H~^H and X-{iu):H^H 

satisfy 

\X+ -Y+\<\X -Y\ and {X' -Y-\ = \{-X)+ - {-Y)+\ < \X -Y\. 

Thus they are both contraction mappings under and can be continuously 

extended to the Banach space {iip, 

We can define the partial order ">" in this Banach space. 

Definition 4.3 An element X in {H, ||-||) is said to be nonnegative, or X >Q, 
0<X,ifX^ X+. We also denote by X >Y, or Y < X , if X-Y >0. 

It is easy to check that X >Y and Y > X imply X = Y on {T-Lp, \ \-\\p)- 
For each X,Y G H, note that 

\E[X]-nY]\<n\X-Y\]<\\X-Y\\^. 
We then can define 

Definition 4.4 The sublinear expectation E[-] can be continuously extended 
to ("Hp, II 'lip) on which it is still a sublinear expectation. We still denote by 

in,'Hp,E). 

Let (r2,H,Ei) be a nonlinear expectation space. Ei is said to be dominated 
by Elf 

El [X] - El [Y] < E[X ~Y] forX,YeH. 

From this we can easily deduce that \Ei[X] ~ E-^[y]| < E\\X — y|], thus the non- 
linear expectation E\\\ can be continuously extended to ('Hp,||-||p) on which it 
is still a nonlinear expectation. We still denote by (n^Hp^Ei). 

Remark 4.5 It is important to note that Xi,--- ,Xn € H does not imply 
ip{Xi, ■ ■ ■ ,Xn) e H for each ip £ C;.i,ip(R")- Thus, when we talk about the no- 
tions of distributions, independence and product spaces on (fij'HjE), the space 
C'i.Lip(]R") is replaced by Cb.Lip(R") unless otherwise stated. 

Exercise 4.6 Prove that the inequalities . ((i^ . ((i?TU|) still hold for {il, U, E). 
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§5 Coherent Measures of Risk 

Let the pair {fl^H) be such that is a set of scenarios and H is the collection 
of all possible risk positions in a financial market. 

If X e "H, then for each constant c, X W c, X A c are all in H. One typical 
example in finance is that X is the tomorrow's price of a stock. In this case, any 
European call or put options with strike price K of forms {S — K)~^, {K — 5)+ 
are in %. 

A risk supervisor is responsible for taking a rule to tell traders, securities 
companies, banks or other institutions uiidc^r his supervision, which kind of risk 
positions is unacceptable and thus a minimum amount of risk capitals should 
be deposited to make the positions acceptable. The collection of acceptable 
positions is defined by 

A = {X X ]s acceptable}. 
This set has meaningful properties in economy. 

Definition b.l A set A is called a coherent acceptable set if it satisfies 

(i) Monotonicity: 

X gA, Y>X imply Y G A. 

(ii) G A but -1 ^ A. 

(iii) Positive homogeneity 

X e A implies XX e A for \> 0. 

(iv) Convexity: 

X,Y eA imply aX + {l- a)Y eA for a e [0, 1]. 

Remark 5.2 (iii)+(iv) imply 

(v) Sublinearity: 

X,Y G A ^ HX + uY G A for fj.,v>0. 

Remark 5.3 If the set A only satisfies (i),(ii) and (iv), then A is called a 
convex acceptable set. 

In this section we mainly study the coherent case. Once the rule of the 
acceptable set is fixed, the minimum requirement of risk deposit is then auto- 
matically determined. 

Definition 5.4 Given a coherent acceptable set A, the functional p{-) defined 
by 

p{X) = pAiX) := inf{m eR: m + X &A}, X eH 
is called the coherent risk measure related to A. 
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It is easy to see that 

p{X + p{X)) = Q. 



Proposition 5.5 p(-) is a coherent risk measure satisfying four properties: 

(i) Monotonicity: If X >Y then p{X) < p(Y). 

(ii) Constant preserving: p(l) = — p(— 1) = — 1. 

(iii) Sub-additivity: For each X,Y H, p{X + Y) < p{X) + p{Y). 

(iv) Positive homogeneity: p{\X) = Xp{X) for A > 0. 

Proof, (i), (ii) are obvious. 
We now prove (iii). Indeed, 

p{X + Y)= inf {m e M : m+{X + Y) eA} 

= inf {m + n : m,n eR, (rn + X) + (n + Y) E A} 

< mi{m eR: m + X E A} + inf{n E M. : n + Y E A} 

=p{X) + p{Y). 

To prove (iv), in fact the case A = is trivial; when A > 0, 

p{XX) = inf{m E R : m + XX E A} 

= A inf {n E R : n + X E A} = Xp{X), 

where n = m/X. □ 
Obviously, if E is a sublinear expectation, we define p{X) E[— X], then p 
is a coherent risk measure. Inversely, if p is a coherent risk measure, we define 
E,[X] :— p{—X), then E is a sublinear expectation. 

Exercise 5.6 Let p(-) be a coherent risk measure. Then we can inversely define 

Ap-.^ {X eH: p{X) < 0}. 
Prove that Ap is a coherent acceptable set. 

Notes and Comments 

The sublinear expectation is also called the upper expectation (see Huber (1981) 



6l| in robust statistics), or t he up per prevision in the theory of imprecise prob- 
abilities (see Walley (1991) ll2C| and a rich literature provided in the Notes 
of this book). To our knowledge, the Representation Theorem 12.11 was firstly 
obtained for the case where 17 is a finite set by ^61], and this theorem was redis- 
covered independently by Artzner, Delbaen, Eber and Heath (1999) [s] and then 
by Delbaen (2002) 35] for the general H.. A typical example of dynamic nonlin- 
ear expectation, called ^-expectation (small g), was introduced in Peng (1997) 
[9^ in the framework of backward stochastic differential equations. Readers are 
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referred to Briand, Coquet, Hu, Memin and Peng (2000) [14], Chen (1998) 
Chen and Epstein(2002) [H, Chen, Kulperger and Jiang (2003) Chen 
and PenK(1998) ^ and (2000) Coquet, Hu, Memin and Peng (2001) ^ 
(2002) ^ , Jiang (2004) (69^, Jiang and Chen (2004) ^ Peng (1999) 
^ and (2004) [g^, Peng and Xu (2003) [107|] and Rosazza (2006) for the 
further development of this theory. It seems that the notions of distributions 
and independence under nonhnear expectations were new. We think that these 
notions are perfectly adapted for the further development of dynamic nonlinear 
expectations. For other types of the related notions of distributions and inde- 
pendence under nonlinear expe ctations or non-additive probabilities, we refer 
to the Notes of the book |120l | and the references listed in Marinacci (1999) 
and Maccheroni and Marinacci (2005) [sJl- Coherent risk measures can be 



also regarded as sublinear expectations defined on the space of risk positions 
in financial market. This notion was firstly introduced in Readers can be 
referred also to the well-known book of Follmer and Schied (2004) [5l| for the 
systematical presentation of coherent risk measures and con vex r isk measures. 
For the dynamic risk measure in continuous time, see [97| or |ll2| . Barrieu and 
El Karoui (2004) Q using ^-expectations. Super-hedging and super pricing (see 
El Karoui and Quenez (1995) [S] and El Karoui, Peng and Quenez (1997) ^) 
are also closely related to this formulation. 



Chapter II 



Law of Large Numbers and 
Central Limit Theorem 

In this chapter, we first introduce two types of fundamentally important distri- 
butions, namely, maximal distribution and G-normal distribiition, in the theory 
of sublinear expectations. The former corresponds to constants and the lat- 
ter corresponds to normal distribution in classical probability theory. We then 
present the law of large numbers (LLN) and central limit theorem (CLT) un- 
der sublinear expectations. It is worth pointing out that the limit in LLN is a 
maximal distribution and the limit in CLT is a G-normal distribution. 

§1 Maximal Distribution and G-normal Distri- 
bution 

We will firstly define a special type of very simple distributions which are fre- 
quently used in practice, known as "worst case risk measure" . 

Definition 1.1 (maximal distribution) A d- dimensional random, vector r] = 
(771, • • • .rjd) on a sublinear expectation space is called maximal dis- 

tributed if there exists a bounded, closed and convex subset F C M'^ such that 

E[ip{r])] =mBx.ip{y). 

Remark 1.2 Here T gives the degree of uncertainty of r]. It is easy to check 
that this maximal distributed random vector r] satisfies 

ar] + bfj = {a + b)r] for a,b>0, 



where fj is an independent copy of rj. 
characterizes a maximal distribution, 
case risk measure" in finance. 



We will see later that in fact this relation 
Maximal distribution is also called "worst 
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Remark 1.3 When d ~ 1 we have T — [/i, where = ^[77] and /i = —£[—77]. 
The distribution of rj is 

¥,,[ip]^E[ip{i^)]= sup ip{y) foripeCi.Lip{R). 

Recall a well-known characterization: X = N{0, S) if and only if 

aX + bX = \/a^ + b^X for a,b>0, (1.1) 

where X is an independent copy of X. The covariance matrix E is defined by 
E = E[XX^]. We now consider the so called G-normal distribution in probabil- 
ity model uncertainty situation. The existence, uniqueness and characterization 
will be given later. 

Definition 1.4 (G-normal distribution) A d- dimensional random vector X — 
{Xi,--- ,Xd)^ on a sublinear expectation space {il^T-L^K) is called (centralized) 
G-normal distributed if 

aX + bX = Vo^+^X for a, 6 > 0, 
where X is an independent copy of X . 

Remark 1.5 Noting that E[X + X] ^ 2E[X] and E[X + X] = E[V2X] = 
V2E[X], we then have E[X] = 0. Similarly, we can prove that E[—X] = 0. 
Namely, X has no mean-uncertainty. 

The following property is easy to prove by the definition. 

Proposition 1.6 Let X be G-normal distributed. Then for each A e M™^'', 
AX is also G-normal distributed. In particular, for each a £ M'^, {a,X) is a 
I- dimensional G-normal distributed random variable, but its inverse is not true 
(see Exercise ] 1.15\) . 

We denote by S{d) the collection of all d x d symmetric matrices. Let X 
be G-normal distributed and rj be maximal distributed d-dimensional random 
vectors on (fij'HjE). The following function is very important to characterize 
their distributions: 

G{p,A):^E[^{AX,X) + {p,T^)], {p,A)eR''x§{d). (1.2) 

It is easy to check that G is a sublinear function monotonic in A E E>{d) in the 
following sense: for each p,p eM.'^ and A,Ag §{d) 

G{p + p,A + A) <G{p,A) + Gip,A), 

G{Xp,XA) =AG(p,_A), VA>0_, (1.3) 
G{p,A) >G{p,A), ifA>A. 



§1 Maximal Distribution and G-normal Distribution 
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Clearly, G is also a continuous function. By Theorem 12.11 in Chapter U there 
exists a bounded and closed subset F C K'' x M'^^'' such that 

G{p, A) = sup [UrlAQQ^] + {p, q)] for {p, A) eR'' x §(d). (1.4) 
(■7,Q)Gr ^ 

We have the following result, which will be proved in the next section. 

Proposition 1.7 Let G : M"^ x S((i) -> R &e a given sublinear and continuous 
function, monotonic in A G in the sense of (jl.3p . Then there exists a G- 
normal distributed d- dimensional random vector X and a maximal distributed 
d-dimensional random vector rj on some sublinear expectation space (51,'H,E) 
satisfying (jl.2l) and 

(aX + bX, a^ri + b^fj) = {\/a'^ + b^X, {a^ + b'^)r]), fora,b>0, (1.5) 
where {X,fj) is an independent copy of{X,ri). 

Definition 1.8 The pair {X,ri) satisfying (II. 5p is called G- distributed. 
Remark 1.9 In fact, if the pair {X,ri) satisfies (|1.5p . then 

aX + bX ^ \/a^ + b^X, arj + bfj^{a + b)r] for a,b>0. 

Thus X is G-normal and rj is maximal distributed. 

The above pair {X, rf) is characterized by the following parabolic partial 
differential equation (PDE for short) defined on [0, oo) X M'' X K'' : 

dtu-G{DyU,Dlu)^Q, (1.6) 

with Cauchy condition u\t=o = (p, where G : x E>{d) — >■ R is defined by 
(fO]) and D^u = (9^,^^.w)fj=i, Du = {d^,u)f^^. The PDE (dH) is called a 
G-equation. 

In this book we will mainly use the notion of viscosity solution to describe 
the solution of this PDE. For reader's convenience, we give a systematical intro- 
duction of the notion of viscosity solution and its related properties used in this 
book (see Appendix C, Section 1-3). It is worth to mention here that for the case 
where G is non-degenerate, the viscosity solution of the G-equation becomes a 
classical G^'^ solution (see Appendix C, Section 4). Readers without knowledge 
of viscosity solutions can simply understand solutions of the G-equation in the 
classical sense along the whole book. 

Proposition 1.10 For the pair (X^rj) satisfying (|1.5p and a function ip G 
Ci.LtpiR'^ X M''), we define 

u{t, x, y) -.^ E[ip{x + ViX, y + tr])], {t, x, y) G [0, oo) x K'' x R''. 
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Then we have 

u{t + s,x,y) =E[u{t,x + y/sX,y + ST])], s > 0. (1.7) 

We also have the estimates: for each T > 0, there exist constants C, fc > such 
that, for all t,s G [0, T] and x, x,y,y G R"^, 

\u{t,x,y) - u{t,x,y)\ < C(l + \x\'' + ly]" + + \y\''){\x - x\ + \y - y\) (1.8) 

and 

\uit,x,y) - u{t + s, X, y)\ < C(l + \x\'' + \y\'')is + \s\'/^). (1.9) 

Moreover, u is the unique viscosity solution, continuous in the sense of p.8p 
and ILH), of the PDE (dH). 

Proof. Since 

u{t, x, y) - u{t, X, y) = E[.^(a; + ViX, y + trj)] ~ E[(f{x + VtX, y + ti])] 

< E[ip{x + VtX, y + tr])- tp{x + ViX, y + trj)] 

< E[Ci(l + \X\^ + M'^ + \x\'' + \y\' + + \yn 

X i\x-x\ + ly-yl) 

< C{1 + \x\'^ + \y\' + \x\'^ + \y\''){\x - x\ + \y - y\), 

we have (jl.Sp . 

Let {X,fj) be an independent copy of {X,ri). By (|1.5|) . 

u{t + s, X, y) = EYp{x + Vt + sX, y+{t + s)vi)] 

= EYp{x + VsX + ^^tX, y + srj + tf/)] 

= ]E[E[(^(a; + ^fsx + VtX, y + sy + i^)](2,y)=(x,,7)] 

= E[u{t,x + VsX, y + srj)] , 

we thus obtain (|1.7p . From this and (|1.8p it foUows that 

u{t + s, X, y) — u{t, X, y) = E[u{t, x + V^X, y + sry) — u{t, x, y)] 
< E[Ci(l + Ixl^ + + IXI" + \v\')iV~s\X\ + sm, 

thus we obtain (|1.9p . 

Now, for a fixed it,x,y) e {0,oo) xW^x W^, let ip e C^^'^([0,oo) x K'^ x R'^) 
be such that ip > u and Tp{t,x,y) = u(t,x,y). By (|1.7p and Taylor's expansion, 
it follows that, for S £ {0,t), 

< EiiPit -S,x + VSX, y + Sri)- tP{t, x, y)] 
<C{S'/^+5^)-dMt,x,y)5 

+ E[{D,^it, X, y), X)VS+ {Dy^it, x,y),7])d + ^ {Dim y)X, X) S] 

= -dt^t, X, y)6 + E[{Dym, X, y), 77) + ^ (l?^^(i, x, y)X, X)]5 + C((53/2 ^ ^2) 
- -dt^, X, y)5 + SG{Dy^, Dl^j){t, x, y) + C(j3/2 + ,52)^ 
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from which it is easy to check that 

[dti;-G{Dy^,Dli,)]{t,x,y)<0. 

Thus w is a viscosity subsolution of (|1.6I) . Similarly we can prove that u is a 
viscosity supersolution of (|1.6|) . □ 

Corollary 1.11 If both {X,ri) and {X,fj) satisfy ()1.5|) with the same G, i.e., 
Gip, A) := E[i {AX, X) + {p,rj)]=E[^{AX,X) + (p, fj)] for [p, A)eR''x §{d) , 

then {X,'r]) = (X,fj). In particular, X = ~X. 
Proof. For each ip £ Ci.Lip{W^ x R''), we set 
u{t, X, y) := ¥.[ip{x + ^tX, y + t^)\, 

u{t,x,y) := E[ip{x + VtX, y + tfj)] , (i, x, y) S [0, oo) x M"* x M"*. 

By Proposition ll.lOl both u and u are viscosity solutions of the G-equation (|1.6p 
with Cauchy condition u|f=o = u\t=o — y^- It follows from the uniqueness of the 
viscosity solution that u = u. In particular, 

E[^{X,rj)]=E[ipiX,f})]. 

Thus {X,rj) = {X,fj). □ 

Corollary 1.12 Let {X,r]) satisfy ()1.5p . For each ijj E Ci.Lip(K'^) we define 

v{t, x) E[i;{{x + ViX + trj)], {t, x) £ [0, oo) x R'*. (1.10) 

Then v is the unique viscosity solution of the following parabolic PDE: 

dtv-G{D,v,Dlv)={), w|t=o = ^. (1.11) 

Moreover, we have v{t,x + y) = u{t,x,y), where u is the solution of the PDE 
(jl.6p with initial condition u(t, x, y)\t=o — ''P{x + y). 

Example 1.13 Let X be G-normal distributed. The distribution of X is char- 
acterized by 

u{t,x)^E[ip{x + ^tX)], ip e Ci.MpiM.'^). 

In particular, E[ip[X)\ = u(l, 0), where u is the solution of the following parabolic 
PDE defined on [0, oo) x R'' : 

dtu-G{D'^u)^Q, u\t=o = ip, (1.12) 

where G = Gx{Al) : S{d) R is defined by 

G{A) ■.= ^E[{AX,X)], Ae§(d). 
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The parabolic PDE (|1.12p is called a G-heat equation. 

It is easy to check that G is a sublinear function defined on S{d). By Theorem 
\2.1\ in Chapter\^ there exists a bounded, convex and closed subset @ C S{d) such 
that 

Ie[{AX, X)] = G{A) = i sup tr[AQ], A e §(d). (1.13) 
2 2 Qge 

Since G{A) is monotonic: G{Ai) > G{A2), for Ai > A2, it follows that 

e C §+(d) = {6* e §(d) : 6* > 0} = {BB^ : B e R'^'"^}, 

where R'^^'^ is the set of all d x d matrices. If Q is a singleton: Q — {Q}, then 
X is classical zero-mean normal distributed with covariance Q. In general, O 
characterizes the covariance uncertainty of X . We denote X = A^({0} x Q) 
(Recall equation (jl.4p . we can set {q,Q) G {0} x 

When d = 1, we have X ^ N{{{)} X [ct^jCt^]) (We also denoted by X — 
N{0, [a^o-^]);, where = E[X'^] and = -£[-^2]. The corresponding G- 
heat equation is 

dtu- ^{^^{dl^u)+ ~ a^idl^u)-) = 0, u\t=o = tp. 

For the case > 0, this equation is also called the Barenblatt equation. 
In the following two typical situations, the calculation of E[(^(X)] is very easy: 
• For each convex function we have 

E[ipiX)] = -= / pia'y) exp(-^)dy. 
V zvr J -00 ^ 

Indeed, for each fixed t > 0, it is easy to check that the function u(t, x) 
IS convex m x: 

u(t, ax + (1 — oi)y) = ¥,[ip{ax + (1 — a)y + VtX)] 

< aE[ip{x + ViX)] + (1 - a)E[ip(x + ViX)] 
= au{t, a;) + (1 — a)u{t, x). 

It follows that (d'^^u)^ = and thus the above G-heat equation becomes 

—2 

dtU = —dl^u, u\t^Q = ip. 



• For each concave function Lp, we have 

1 



E[p{X)] 



'l-K J- 



y 

ip{qfy) exp(- — )dy 
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In particular, 

E[X] = E[-X] = 0, E[X^] = -E[-X2] = £2 

and 

E[X^] = 3o^, ~E[~X^] = 3^4 . 

Example 1.14 Leti] be maximal distributed, the distribution of rj is character- 
ized by the following parabolic PDE defined on [Q.co) x : 

dtu~g{Du)^0, u|t=o = <^, (1.14) 

where g = (7,,(p) : M*^ — > M is defined by 

5„(p) :=E[(p,77)], peR"*. 

It is easy to check that is a sublinear function defined on W^. By Theorem 
\2.1\ in Chapter\^ there exists a bounded, convex and closed subset Q C R"^ such 
that 

5(p) = sup(p,g), PGR". (1.15) 
qee 

By this characterization, we can prove that the distribution of rj is given by 
F^M=E[(^(77)] = sup</j(t;) = sup / ^ix)Sy{dx), G G.L,p(R'*), (1.16) 

where is Dirac measure. Namely it is the maximal distribution with the 
uncertainty subset of probabilities as Dirac measures concentrated at 0. We 

denote rj = N{Q x {0}) (Recall equation ()1.4|) . we can set {q,Q) € 9 x {0}/ 
In particular, for d = 1, 

9ri{p) ■■= Hpv] = MP^ - MP", P e R, 

where p. — £[77] and ji = — E[— ry]. The distribution of rj is given by I11.16\) . We 
denote rj = N{[ji,ji] x {0}). 

Exercise 1.15 We consider X = {Xi,X2), where Xi = A^({0} x [ct^,ct^]) with 
a > CT, X2 is an independent copy of Xi . Show that 

(1) For each a G R^ {a,X) is a 1-dimensional G-normal distributed random 
variable. 

(2) X is not G-normal distributed. 

Exercise 1.16 Let X be G-normal distributed. For each ip g Ci,Lip{R'^), we 
define a function 

u{t, x) E[^(a; + ViX)], {t, x) G [0, 00) x R^. 

Show that u is the unique viscosity solution of the PDE hi. 12]) with Cauchy 
condition u\t=o — f. 
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Exercise 1.17 Let rj be maximal distributed. For each ip E C;.Lip(R'^), we 
define a function 

uit, y):=EMy + t7])] , (t, y) G [0, oo) x M"^. 

Show that u is the unique viscosity solution of the PDE \1.14^ with Cauchy 
condition u\t=o — (p. 

§2 Existence of G-distributed Random Variables 

In this section, we give the proof of the existence of G-distributed random 
variables, namely, the proof of Proposition II .71 

Let G :M.'^ X §{d) R be a given sublinear function monotonic in A e §{d) 
in the sense of ()1.3|) . We now construct a pair of d-dimensional random vectors 
{X,r]) on some sublinear expectation space (r2,H,E) satisfying (|1.2p and (|1.5p . 

For each cp g C;.Lip(R^''), let it = it''' be the unique viscosity solution of 
the G-equation (HH) with u'^\t=Q = (fi. We take h = R'^'', H = Ci.Ltp{M.'^'^) 
and uj = {x,y) £ M^'^. The corresponding sublinear expectation E[-] is defined 
by E[^] = u'''(l,0, 0), for each ^ G "H of the form ^(w) = ('/'(a;, 2/))(x,y)6R2'i & 
Ci.LipiR'^'^)- The monotonicity and sub-additivity of u'^ with respect to (p are 
known in the theory of viscosity solution. For reader's convenience we provide a 
new and simple proof in Appendix C (see Corollary 12 .41 and Corollarv l2.5l) . The 
constant preserving and positive homogeneity of E[-] are easy to check. Thus 
the functional E[-] : "H — > M forms a sublinear expectation. 

We now consider a pair of d-dimensional random vectors (X, rj)(uj) = {x, y). 
We have 

nviXM = ^i''(l,0,0) for ^ e G.l,p(R2'^). 
In particular, just setting pr>{x, y) = \ {Ax, x) + {p, y), we can check that 

{t, X, y) = G{p, A)t + i {Ax, x) + {p, y) . 

We thus have 

E[i [ax,x) + {p,lf)] - (1,0,0) = G{p,A), {p,A) e R^ X §(d). 

We construct a product space 

(17,-H,E) = (fi X il,-H«)-H,i«)E), 
and introduce two pairs of random vectors 

{X,Vl){u}l,U}2) =U}i, (X,^)(wi,W2) = W2, (Wi,52) G f2 X f2. 

By Proposition 13 . 1 5 1 in Chapter HI {X,ri) ^ {X,r{) and {X,fi) is an independent 
copy of {X, rj). 
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We now prove that the distribution of {X, rj) satisfies condition (11.51) . For 
each (fi G C;.Lip(R^'') and for each fixed A > 0, {x, y) G K^'', since the function v 
defined by v{t, x, y) := u'^^Xt, x + -v/Ax, y + Xy) solves exactly the same equation 
(|1.6p . but with Cauchy condition 

v\t=o = (fix + \/A X •,y + A X •). 

Thus 

E[ip{x + VXX, y + Xt])] = v{l, 0, 0) = u'^(A, x, y). 

By the definition of E, for each t > and s > 0, 

E[<f{ViX + V^X, tij + sf])] = E[E[<f{Vix + ^fsX, ty + sf))] (:,,a)=(x,„)] 

= EK(s, VtX, try)] = ^"''('''■'■'(i, 0, 0) 
= w'^(t + s,0,0) 
= E[ip{^/tT~sX, {t + s)ri)]. 

Namely {^/iX + ^/sX, trj + sf]) = (V* + sX, {t + 5)77). Thus the distribution of 
{X, 77) satisfies condition (|1.5I) . 

Remark 2.1 From now on, when we mention the sublinear expectation space 
(QyHyE), we suppose that there exists a pair of random vectors {X, 77) on E) 
such that {X, 77) is G-distributed. 

Exercise 2.2 Prove that E[X^] > for X = N{{0} x [ff^a^]) with < . 

It is worth to point that E[ip{X)] not always equal to swp^2^^^^2 Ea-[ip{X)] 
for if G Ci,Lip(R), where E„ denotes the linear expectation corresponding to the 
normal distributed density function N[0,a'^). 

§3 Law of Large Numbers and Central Limit 
Theorem 

Theorem 3.1 (Law of large numbers) Let {Yi}'^^ be a sequence of W^- 
valued random variables on a sublinear expectation space (fij'HjE). We assume 

that li+i = Yi and Yi^i is independent from {Yi, • • • , Yi} for each i — 1,2, ■ ■ ■ . 
Then the sequence {Sn}'^=i defined by 

1 " 

Sn '■— / Yi 

n ^-^ 

i=l 

converges in law to a maximal distribution, i.e., 

lim E[(^(^„)] =E[<^(77)], (3.17) 
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for all functions (f € C(M'') satisfying linear growth condition < C(l + 

|a;|)), where rj is a maximal distributed random vector and the corresponding 
sublinear function g : M'' —> R is defined by 

g{p) :=E[(p,ri>], peR". 

Remark 3.2 When d = 1, the sequence {>S'n}^i converges in law to N([fj,, p] x 
{0}), where p — E[Yi] and ^ = — E[— Yi]. For the general case, the sum 
■^'^i^iYi converges in law to N{Q x {0}), where Q C M'* is the bounded, 
convex and closed subset defined in Example \1.14 If we take in particular 
(p{y) = dQ{y) = inf{|a; — y\ : x G &}, then by p.l7p we have the following 
generalized law of large numbers: 

1 " 

lim E[de(-V sup de(6i) 0. (3.18) 

// Yi has no mean-uncertainty, or in other words, Q is a singleton: Q = {9}, 
then p.lSp becomes 

1 " 

lim E[|- Vy, -6l|] =0. 

n^oo 77, ^- — ^ 
1=1 

Theorem 3.3 (Central limit theorem with zero-mean) Let {Xi}^-^ be a 
sequence of W'' -valued random variables on a sublinear expectation space {il,H, E). 

We assume that Xi^i ^ Xi and Xi^i is independent from {Xi, ■■ ■ ,Xi}foreach 
i — 1, 2, • • • . We further assume that 

E[Xi] = E[-Xi] = 0. 

Then the sequence {Sn}^=i defined by 

1 " 

Sn ■= —F= / Xi 

V" ^ 

^ z— 1 

converges in law to X, i.e., 

lim E[^(5„)] =E[v.(X)], 

for all functions ip £ C(U.'^) satisfying linear growth condition, where X is a 
G-normal distributed random vector and the corresponding sublinear function 
G : §(d) ^ R is defined by 

G{A) :=E[i(AXi,Xi)], AeS{d). 

Remark 3.4 When d = 1, the sequence {Sn\^=i converges in law to N{{{)\ x 
[ct^,ct^]), where a'^ = E[X^] and = — E[— Xj^]. In particular, ifa'^ = , then 
it becomes a classical central limit theorem. 
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The following theorem is a nontrivial generalization of the above two theo- 
rems. 

Theorem 3.5 (Central limit theorem with law of large numbers) Let 

{{XijYi)}'^^ be a sequence of W^' x R'^-valued random vectors on a sublin- 

ear expectation space (f2,H,E). We assume that (Xj+i,yi+i) = (Xi^Yi) and 
li+i) is independent from {(^i, Yi), • • • , {Xi, Yi)} for each i = 1, 2, • • • . 
We further assume that 

E[Xi] = E[-Xi] = 0. 
Then the sequence {Sn}^=i defined by 

^ Vn n 

2—1 ^ 

converges in law to X + rj, i.e., 

lim E[<fi{Sn)] = E[ip{X + Tj)], (3.19) 

n— foo 

for all functions cp € C(R'*) satisfying a linear growth condition, where the pair 
{X, rj) is G -distributed. The corresponding sublinear function G ; x §((i) — )• M 
is defined by 

G(p,A) :=E[(p,Fi) + i(AXi,Xi)], Ae^d), pen". 

Thus E[(/3(X + 77)] can be calculated by Corollaru ll.lSl 

The following result is equivalent to the above central limit theorem. 

Theorem 3.6 We make the same assumptions as in Theorem \3.5l Then for 
each function (p e G (M'' x M.'^) satisfying linear growth condition, we have 

hm Eb(^^,5]^)]=E[(p(X,^)]. 
n->oo — ' ,/n ^ — ' n 

i=l V i=l 

Proof. It is easy to prove Theorem [33] by Theorem 13.61 To prove Theorem l3.6l 
from Theorem 13.51 it suffices to define a pair of 2(i-dimensional random vectors 

X, = {X,,Q), = (0,y,) fori = l,2,--- . 

We have 

lim E[<p(X]^,f^^)]= lim n^{j2{^ + -))]=n^{x+v)] 

i=l V i=l 2=1 V 

-E[^(X,r;)] 

with X = {X, 0) and f] = (0, 77). □ 
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To prove Theorein l3.5[ we need the following norms to measure the regularity 
of a given real functions u defined on Q = [0, T] x K'': 

\W\\c«.'HQ) = sup \u{t,x)\, 

d 

ll"llci'i(Q) ^ ll"llco.O(Q) + ll'9t"llcO'0(Q) + ^ ll^2;,u||(70.0(Q) , 

i=l 



d 

W r) iiW 

C°.°(Q) • 



h3 = l 

For given constants a,(3 £ (0, 1), we denote 



11 II \u(s, x) — u(t, y)\ 

d 

II'"llci+°'i+'3(Q) = ll"llc°''3(Q) + ll^t^llc".'3(Q) + X/ ll^^i" 



C°./3(Q) ' 
i=l 

d 



If, for example, ||u||pi+a,2+/3(Q) < oo, then u is said to be a C^+^-'^+^-function 
on Q. 

We need the following lemma. 

Lemma 3.7 We assume the same assumvtions as in Theorem ] 3. 5\ We further 
assume that there exists a constant /3 > such that, for each A, A £ §(d) with 
A> A, we have 

V\{AXi,Xi)] - ¥.[{AXi,Xi)] > l3tr[A - A]. (3.20) 
Then our main result (j3.19p holds. 

Proof. We first prove ((3?T9)) for tp e Cb.LipO^-'')- For a small but fixed h > 0, 
let V be the unique viscosity solution of 

dtV + G{DV, D^V) = 0, (i, x) e [0, 1 + /t) x M'', V\t=i+h ^ V- (3.21) 

Since {X,r]) satisfies p.Sp . we have 

V{h,Q) =E[^{X + Tj)], V{l + h,x) = ip{x). (3.22) 

Since (I3.2ip is a uniformly parabolic PDE and G is a convex function, by the 
interior regularity of V (see Appendix C), we have 

||T^|lci+<:./2.2+„([o,i]xR'') < °° for some a e (0, 1). 
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We set (5 = i and So = 0. Then 

ri-l 

V{1, Sn) - V{0, 0) = Y.^V{{t + 1)6, - V{tS, S,)} 

i=0 

4=0 

i=0 

with, by Taylor's expansion, 

JJ = atl^(z<5, 5,)<5+i (i?V(z(5, S,)X,+i,X,+i) 6+(^DViiS, 5,), X.+i + Y,+id 
11 = 5 I [dtV{{i + - atF(z,5, S,+i)]d/3 + [dtV{i5, S,+^) - 9*^(1,5, S0]<5 



with 
Thus 

n— 1 n— 1 n— 1 n— 1 

J|] - E[- ^ /]] < E[F(1, 5„)] - ^(0, 0) < Jl] + ^s]- (3.23) 

1=0 1=0 4=0 1=0 

We now prove that E[X;r=ro^ Jl] = 0. For J], note that 

E[(^DVii6, S,), X,+iVs)] = E[- (^DViiS, S,),X,+iVs)] = 0, 
then, from the definition of the function G, we have 

E[J]] = E[dtV{iS, 5,;) + G{DV{i6, Si),D^V{i5, Si))]5. 

Combining the above two equahties with dtV + G{DV, D^V) = as well as the 
independence of (X^+i, Yi+i) from {{Xi,Yi), ■ ■ ■ , {Xi,Yi)}, it follows that 



n-2 



E[^ Jl] = E[^ JJ] ^ . . . = 0. 

1=0 1=0 

Thus (|3.23p can be rewritten as 



-E[- Y n] < nV{l, Sn)] - V{0, 0) < E[^ II]. 

i=0 i=0 
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But since both dtV and D^V are uniformly ^-holder continuous in t and a- 
holder continuous in x on [0, l]xR'^, we then have 

It follows that 

E[|/l|] < C(5i+"/2(i+E[|Xi|2+" + |yi|2+«]). 

Thus 

+ E[|Xi|2+" + < E[y(i,5„)] - 1/(0,0) 

n 

< C(-)"/2(l+E[|Xi|2+" + |ri|2+«]). 

n 

As 71 — > cx), we have 

lim E[y(l, S-n)] = V{0, 0). (3.24) 

n— >-oo 

On the other hand, for each t, f G [0,1 + /i] and x G M*^, we have 

\Vit,x)-V{t',x)\<C{V\t^\ + \t-t'\). 
Thus |T/(0,0) - V{h,0)\ < C{Vh + h) and, by (HH]), 

|E[-t^(l,5„)]-E[(^(5„)]| = |E[1/(1,5„)] -E[V(l + /i,5„)]| < C(Vh + h). 
It follows from (1X22)) and ((3?24| that 

limsup|E[(p(5„)] -E[(/7(X + r;)]| < 2C{y^+h). 

n^oo 

Since h can be arbitrarily small, we have 

lim E[ip{Sn)]=E[ip{X + r])]. 

n— >-oo 

□ 

Remark 3.8 From the proof we can check that the main assumption of identical 
distribution of {XijYi}'^^ can be weaken to 

nip, Yi) + \ {AX,, X,)] = G{p, A), I = 1, 2, • • • , p g R"^, A e S(d). 

Another essential condition is E[|X,p+'5] +E[|yi|i+'^] < C for some 5>Q. We 
do not need the condition E[|Xi|"] +E[|yi|"] < oo for each 7i G N. 

We now give the proof of Theorem 13.51 
Proof of Theorem [3751 For the case when the uniform elliptic condition p. 201) 
does not hold, we first introduce a perturbation to prove the above convergence 
for £ Cf,.Lip(IR'^). According to Dcfinition l3.14l and Proposition l3 . 1 5l in Chap I, 



§3 Law of Large Numbers and Central Limit Theorem 



31 



we can construct a sublinear expectation space (fi, %, E) and a sequence of three 

random vectors {{Xi, Yi, Ki)}"!^ such that, for each n = 1, 2, • • • , {{Xi, Yi)}^^-^ = 
{{X,,Yt)}f^i and (X„+i, is independent from {(X^, Y^, Ki)}"^]^ and, 

moreover, 

E[iPiX,, r„ R,)] = (27r)-'^/2 f E[V;(X„ K„ x)]e-^'=^" /^dx for G Cl.L^piM.'''"'). 

We then use the perturbation = + eRi for a fixed e > 0. It is easy to 
see that the sequence {{Xf, Yi)}^-^ satisfies all conditions in the above CLT, in 
particular, 

Ge(p,A) -.^El^iAXtat) + {P,Y^)] = Gip,A) + ^tT[A]. 
Thus it is strictly elliptic. We then can apply Lemma 13.71 to 




and obtain 

lim ElifiS'J] = E[ip{X + fj + eR)], 
where ((X,/^), (?7,0)) is (5-distributed under E[-] and 

G{p,A):^E[^{A{X,,R,f,{X^,R,f) + {p,{YuOf)], AeS{2d), peR^''. 

By Proposition 11.61 it is easy to prove that {X + eR, fj) is Ge-distributed and 
{X,fj) is G-distributed. But we have 

\E[ip{S^)] - EMS^M = \E[^{S'^ - eJ„)] - E[^(5^)]| 
< eGE[| J„|] < C'e 

and similarly, 

\E[ip{X + r;)] - E[ip{X + fj + eR)]\ = \E[ip{X+fj)] - E[ip{X+f] + eR)]\ < Ce. 
Since e can be arbitrarily small, it follows that 

lim E[(p(5„)] = E[ip{X + ry)] for ^ G Cb.Mp{^'')- 

n— f C30 

On the other hand, it is easy to check that supnEljS'np] + E\\X + ryp] < oo. 
We then can apply the following lemma to prove that the above convergence 
holds for (^eG(M'') with linear growth condition. The proof is complete. □ 

Lemma 3.9 Let (fi, %, E) and (f2, %, E) be two sublinear expectation spaces and 
let Yn (z H and Y ^H, n — 1,2, ■ ■ ■ , be given. We assume that, for a given p > 
1, sup„E[|r„|P] +E[|y|P] < oo. If the convergence \imn^ooHyy{Yn)] = E[ip{Y)] 
holds for each ip E Cb.Lip(R'^), then it also holds for all functions ip G G(IR.'') 
with the growth condition \^{x)\ < G(l + 
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Proof. Wc first prove that the above convergence holds for ip G Cb(IR.'^) with 
a compact support. In this case, for each £ > 0, we can find a, (p €: Cb.npi^'^) 
such that sup^gRd |<^(x) — (p{x)\ < |. We have 

|E[^(r„)] -E[^(r)]| < |E[(^(y„)] -E[^(y„)]| + |E[^(r)] -E[^(y)]| 
+ |E[^(y:„)] - < e + |E[^(r„)] - n^{Y)]\. 

Thus hmsup„_^Q^3 |E[<p(y„)] — E[(^(y)]| < e. The convergence must hold since e 
can be arbitrarily small. 

Now let be an arbitrary C(R'^)-function with growth condition < 
C(l + |a;|P~^). For each TV > we can find (pi,ip2 G C{W^) such that = '^1 + ^2 
where ipi has a compact support and ^2{x) = for < A'', and |</'2(a;)| < \<fi{x)\ 
for all X. It is clear that 

< '^^i^ + l^l") for X e R''. 

Thus 

|E[<p(y„)] - E[<p(y)]| = nMYn) + MYn)] - nMY) + My)]\ 

<\E[MYn)]-nMy)]\+n\MYn)\]+n\MY)\] 

< \E[MYn)]-nMY)]\ + — (2 + E[|y„n +E[|yr]) 

<|E[<^i(y„)]-E[<^i(F)]| + ^, 

where (7 = 2(7(2+sup„E[|F„|P]+E[|y|f]). We thus have lim sup„^^ |E[(p(F„)]- 
IE[(^(y)]| < ^. Since N can be arbitrarily large, E[(p(y^)] must converge to 

E[^(r)]. □ 

Exercise 3.10 Let Xi £ H,i = 1,2, ■■■ , be such that Xi+i is independent from 
{Xi, ■ ■ ■ , Xi}, for each i = 1,2, ■ ■ ■ . We further assume that 

E[X,] = -E[^X,] = 0, 

lim E[Xf] = <oo, hm -~E[-Xf] = a^, 

i—¥oo i— >-oo 

E[|Xip+*] < M for som,e 5 > and a constant M. 
Prove that the sequence {Sn}^=i defined by 



1 " 



^'^ i=i 

converges in law to X, i.e., 

lim E[^{Sn)] = E[ip{X)] for ip € Cb,iip{R), 

n^oo 

where X N({0} x [ct^ct^j) 

In particular, if = 21^ , it becomes a classical central limit theorem. 



§3 Law of Large Numbers and Central Limit Theorem 
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Notes and Comments 

The co nten ts of this chapter are mainiy from Peng (2008) tlQSj (see also Peng 
(2007) [io3). 

The notion of G-normal distribution was fi rstly introduced by Peng (2006) 

|l00l | for l-dimensional case, and Peng (2008) [l04j for multi-dimensional case. 
In the classical situation, a distribution satisfying equation is said to be 

stable (see Levy (1925) [t^ and (1965) [ll])- In this sense, our G-normal dis- 
tribution can be considered as the most typical stable distribution under the 
framework of sublinear expectations. 

Marinacci (1999) jsH used different notions of distributions and indepen- 
dence via capacity and the corresponding Choquet expectation to obtain a law 
of large numbers and a central limit theorem for non-additive probabilities (see 
also Maccheroni and Marinacci (2005) {83| ). But since a sublinear expectation 
can not be characterized by the corresponding capacity, our results can not be 
derived from theirs. In fact, our results show that the limit in CLT, under 
uncertainty, is a G-normal distribution in which the distribution uncertainty is 
not just the parameter of the classical normal distributions (see Exercise I2.2p . 

The notion of viscosity solutions plays a basic role in the definition and 
properties of G-normal distribution and maximal distribution. This notion was 
initially introduced by Crandall and Lions (1983) [29J. This is a fundamentally 
important notion in the theory of nonlinear parabolic and elliptic PDEs. Read- 
ers are referred to Crandall, Ishii and Lions (1992) [s^ for rich references of the 
beautiful and powerful theory of viscosity solutions. For books on the theory of 
viscosity solutions and the related HJB equations, see Barles (1994) [§], Fleming 
and Soner (1992) [H as weh as Yong and Zhou (1999) flii ]. 

We note that, for the case when the uniform elliptic condition holds, the vis- 
cosity solution (|1.10p becomes a classical G^+'^'^'*'"-solution (see Krylov (1987) 
ll and the recent works in Cabre and Caffarelh (1997) [i3| and Wang (1992) 
ligf ). In l-dimensional situation, when > 0, the G-equation becomes the 
following Barenblatt equation: 

dtu + j\dtu\ = Au, I7I < 1. 

This equation was first introduced by Barenblatt (1979) (see also Avellaneda, 
Levy and Paras (1995) Q). 
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Chapter III 



G-Brownian Motion and 
Ito's Integral 

The aim of this chapter is to introduce the concept of G-Brownian motion, to 
study its properties and to construct Ito's integral with respect to G-Brownian 
motion. We emphasize here that our definition of G-Brownian motion is con- 
sistent with the classical one in the sense that if there is no volatility uncer- 
tainty. Our G-Brownian motion also has independent increments with identical 
G-normal distributions. G-Brownian motion has a very rich and interesting 
new structure which non-trivially generalizes the classical one. We thus can 
establish the related stochastic calculus, especially Ito's integrals and the re- 
lated quadratic variation process. A very interesting new phenomenon of our 
G-Brownian motion is that its quadratic process also has independent incre- 
ments which are identically distributed. The corresponding G-Ito's formula is 
obtained. 

§1 G-Brownian Motion and its Characterization 

Definition 1.1 Let (ri,H,E) be a sublinear expectation space. {Xt)t>Q is called 
a d- dimensional stochastic process if for each t > 0, Xt is a d- dimensional 
random vector in %. 

Let G(-) : §(c?) — > M be a given monotonic and sublinear function. By 
Theorem 12.11 in Chapter [H there exists a bounded, convex and closed subset 
S C §+(d) such that 

G(^) = i sup Ae§(d). 

2 seE 

By Section[2]in Chapter HIl we know that the G-normal distribution A^({0} x S) 
exists. 

We now give the definition of G-Brownian motion. 
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Definition 1.2 A d-dimensional process {Bt)t>a on a suhlinear expectation 
space (rij'HjE) is called a G-Brownian motion if the following properties 
are satisfied: 

(i) BoH =0; 

(ii) For each i, s > 0, the increment Bt+s — Bt is -/V({0} x sT,)- distributed and is 
independent from {Bt^ , Bt^,- ■ ■ , Bt^), for each n e N and < ti < ■ ■ ■ < tn < t. 

Remark 1.3 We can prove that, for each to > 0, (Bt+to ~ -Bto)t>o is a G- 

Brownian motion. For each A > 0, {\~^ Bxt)t>o is also a G-Brownian motion. 
This is the scaling property of G-Brownian motion, which is the same as that 
of the classical Brownian motion. 

We will denote in the rest of this book 

Bf = (a, Bt) for each a = (ai , • • • , a^)'^ G R^. 

By the above definition we have the following proposition which is important 
in stochastic calculus. 

Proposition 1.4 Let {Bt)t>o &e a d-dimensional G-Brownian motion on a 
suhlinear expectation space (r2,'H,E). Then {Bf)t>o is a 1-dimensional Ga- 
Brownian motion for each a eM"^, where Ga{a) = |(o'^aTQ!"'" — c^^a'^Q!"), 
al^^ = 2G(aa^) = E[(a,Bi)2], a\^^ = -2G(-aa^) = -E[-(a, Si)^]. 

In particular, for each t,s>0, Bf^^ — Bf = N{{Q} x [scr^aaf' ■^"'aaT])- 

Proposition 1.5 For each convex function , we have 

1 r°° x^ 

E[^(Bf+, - Bf )] = / ^{x) exp(-— 

For each concave function y and c^^^t > 0, we have 

1 1'°° x^ 

n^{Bt+s - Bt)] = , ^ / ^{x) exp(-— ^ )dx. 

^2nsa\^^ J-oo '^sa_^^^ 

In particular, we have 

E[{Bt - Bf)'] = al^At - s), E[(Bf - Bf)^] = 3a^,.(t - s)\ 

E[-(Bf - Bff] = -a'_,Mt - s), E[-(Bf - Bff] = -3al,,.(t - sf. 

The following theorem gives a characterization of G-Brownian motion. 

Theorem 1.6 Let {Bt)t>o he a d-dimensional process defined on a suhlinear 
expectation space (il, E) such that 

(i) -Bo(w)- 0; 

(ii) For each t, s > 0, Bt+s — Bt and Bg are identically distributed and Bt+s — Bt 
is independent from {Bt^ , Bt^ ,■ ■■ , Bt^), for each n e N and <ti < ■ ■ ■ <tn < 
t. 

(ill) E[Bt] = E[-Bt] = and limt^o E[\B,\^]t-^ = 0. 

Then {Bt)t>o is a G-Brownian motion with G{A) = ^E[(ABi,Bi)], A e S{d). 



§1 G-Brownian Motion and its Characterization 



37 



Proof. We only need to prove that Bi is G- normal distributed and Bt — \/tBi . 
We first prove that 

¥\{ABt,Bt)] = 2G{A)t, A e 

For each given A £ S{d), we set b{t) =E[{ABt, Bt)]. Then &(0) = and \b{t)\ < 
|A|(E[|Bt|3])2/3 ^ as t ^ 0. Since for each t,s>0, 

h{t + s)^ ¥.[{ABt+s.Bt+s)] - t[{A{Bt+s - B, + Bs),Bt+s - B^ + S,)] 

= E[{A{Bt+s - Bs), [Bt+s - B,)) + {ABs,Bs) + 2{A{Bt+s ~ B,),Bs)] 
= b{t) + b{s), 

we have b{t) = b{l)t ^2G{A)t. 

We now prove that Bi is G-normal distributed and Bt == y/tBi. For this, 
we just need to prove that, for each fixed tf G Cb.Lipi^"^), the function 

u{t,x) ■.= E[Lp{x + Bt)], {t,x) e [0,oo) X E'^ 

is the viscosity solution of the following G-heat equation: 

dtU-G{D^u) = 0, u\t=o = ^- (1.1) 

We first prove that u is Lipschitz in x and i-Holder continuous in t. In fact, 
for each fixed t, u{t, •) &Cb.Lip{^'^) since 

\uit, x) - uit, y) I = \E[ip{x + Bt)]~ E[(^(y + Bt)] \ 
<E[\ip{x + Bt) - ip{y + Bt)\] 
<C\x~yl 

where G is Lipschitz constant of ip. 

For each (5 e [0, t], since Bt — Bg is independent from Bs, we also have 

u{t, x) = E[(p{x + Bs + {Bt - Bs)] 

= E[E[^{y + {Bt-Bs))]y=,+Bs], 

hence 

u{t,x) =E[u{t-S,x + Bs)]. (1.2) 

Thus 

\u{t,x)-u{t~S,x)\ = \E[u{t-6,x + Bs)~u{t-6,x)]\ 

< E[\u{t -S,x + Bs) - u{t - S, x)\] 

< E[C\Bs\] < Cy/2G{I)VS. 

To prove that u is a viscosity solution of (jl.ip . we fix {t, x) £ (0, oo) x M'' and 
let V e G^'^([0,oo) X R'^) be such that v > u and v{t,x) = u{t,x). From (IL^ 
we have 

v{t, x) = E[u{t -5,x + Bs)] < E[v{t -~5,x + Bs)]. 
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Therefore by Taylor's expansion, 

< E[v{t ~S,x + Bs)- v{t, x)] 
= E[v{t -5,x + Bs) - v{t, X + Bs) + {v{t, x + Bg) - v{t, x))] 

= E[-dtv{t, x)6 + {Dv{t, x), Bs) + ^{D^v{t, x)Bs,Bs) + h] 
< -dtv{t, x)S + ^E[{D^v{t, x)Bs,Bs)] + E[Is] 



where 



= -dtvit, x)6 + GiD^v{t, x))6 + E[Is] , 

h = I ~[dtv{t- pS,x + Bs)-dtv{t,x)]Sdl3 
Jo 

»i 



{{D^v{t, X + aPBs) - D^v(t, x))Bs, Bs)adpda. 



[(-U'V(i, X -h apns ) — l'^' 

IQ Jo 

With the assumption (iii) we can check that hm^j^o ~ 0, from which we 

get dtv{t,x) — G{D^v{t,x)) < 0, hence it is a viscosity subsolution of We 
can analogously prove that u is a viscosity supersolution. Thus u is a viscosity 
solution and (i?t)t>o is a G-Brownian motion. The proof is complete. □ 

Exercise 1.7 Let Bt be a 1- dimensional Brownian motion, and Bi ^ iV({0} x 
[ct'^jCT^]). Prove that for each m gN, 



n\Bt 



2(to - /y/2Tr m is odd, 

(m — iy.\a"^t~ m is even. 



§2 Existence of C-Brownian Motion 

In the rest of this book, we denote by 17 = Cq (R+) the space of all R'^-valued 
continuous paths {tjJtjteM+i with — 0, equipped with the distance 

oo 

i—l ' 

For each fixed T G [0, oo), we set Q,t {i^ at : ^ G We will consider the 
canonical process Bt{uj) = ujt, t & [0,oo), for w £ il. 
For each fixed T G [0, oo), we set 

L^pirtr) -.^WiBt.^Tr-- ^Bt^^r) : ne N,ti,-- - ,t„ £ [0,oo), ^ € Cl.L^p{^''''") }• 

It is clear that Lip{nt)CLip{ilT) , for t <T. We also set 

oo 

Lip{^) ■— (J Lipiyin). 

n=l 



§2 Existence of G-Brownian Motion 
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Remark 2.1 It is clear that C(.Lip(]R'^^"), Lip{ilT) and Lip(f2) are vector lat- 
tices. Moreover, note that (p,ip G Q.LipiM.'^^"') imply (p ■ i(j G Ci.Lip(IR''^"), 
then X, Y &Lip{Q,T) imply X ■ Y GLip{ClT)- In particular, for each t G [0, oo), 
Bt e Lip{ft). 

Let G(-) : S{d) — > R be a given monotonic and sublinear function. In the 
following, we want to construct a sublinear expectation on {Q,Lip{Q)) such 
that the canonical process (Bt)t>o is a G-Brownian motion. For this, we first 
construct a sequence of d-dimensional random vectors (Ci)i^i on a sublinear 
expectation space (O.'HjE) such that is G-normal distributed and ^j+i is 
independent from (^i , • • • , ^i) for each i = 1,2, ■ ■ ■ . 

We now introduce a sublinear expectation E defined on Lip{Cl) via the fol- 
lowing procedure: for each X G Lip{Cl) with 

X = ip{Bti - Bto , Bt^ - , • • • , Bt^ - Bt„_^) 

for some G G;.Lip(R'^^") and = to < < • • • < in < oo, we set 

E[^(B(, - Bt, ,Bt,-Bt„-- - , Bt^ - J] 

:= Wj[ip{\/t\ — io^l, • • • , \Jtn — tn-\iv)\- 

The related conditional expectation of X = p>(Btx > -^ta — -Bti ,■ ■■ , Bt^ — 
Btn-i) under Slf. is defined by 

t[X\nt^]=E[ip{Bt„Bt,-Bt„--- ,Bt^-Bt„_,)\nt^] (2.3) 
■.= ij{Bt„--- ,Bt,-Bt,_,), 

where 

■4>{X1,--- ,Xj) =E[ip{xi,--- ,Xj,^tj+l -tj^j+i,--- , ^/tn - tn-l^n)]- 

It is easy to check that E[-] consistently defines a sublinear expectation on 
Lip{n) and (i3i)t>o is a G-Brownian motion. Since Z/jj,(f2T)CLjp(f2), ¥,[■] is also 
a sublinear expectation on Lip{nT). 

Definition 2.2 The sublinear expectaiion E[-]: Ljp(17) — > M defined through the 
above procedure is called a G -expectation. The corresponding canonical process 
{Bt)t>o on the sublinear expectation space {ft, Ljp(f2), E) is called a G-Brownian 
motion. 

In the rest of this book, when we talk about G-Brownian motion, we mean 
that the canonical process {Bt)t>o is under G-expectation. 

Proposition 2.3 We list the properties ofE[-\dt] that hold for each X, Y Sijp(f2) 
(i) IfX>Y, then E[X|17t] > E[y|^lt]. 
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(ii) E[77|r2t] = T], for each t e [0, oo) and rj €zLip(flt)- 

(iii) E[x\nt] - E[Y\nt] < E[x - Y\nt]- 

(iv) E[r]X\nt] = V^E[X\nt] + V^t[-X\nt] for each rj £ L.piflt). 

(v) E[E[X|l]t]|l]^] = E[X|f7tAs], in particular, E[E[X\nt]] = E[X]. 

For each X € Lip{VL^), E[X\D,t\ = E[X], where Lip{Q.^) is the linear space of 
random variables with the form 

^{Bt2 - Bt^^Bt^ - Bt^, - ■ ■ , Bt^^^ - Bt^ ) , 
n = l,2,--- , LP eCl.L^p{M.'^'"''), ti,--- ,t„,t„+i e [t,(X)). 

Remark 2.4 (ii) and (iii) imply 

E[X + r]\nt] = E[X\nt] + 7] for r] G UpiVlt). 

We now consider the completion of sublinear expectation space {Vt, Lip{il), E). 

We denote by Lq{V,), p > 1, the completion of Lip{Vl) under the norm 
:= (E[\X\P])^/P. Similarly, we can define i:^(f2T), i^G(^T) and £^(1^*). 
It is clear that for each < t < T < oo, L^{flt) C L^(r2T) ^ L^ci^)- 

According to Sec H] in ChapU E[-] can be continuously extended to a sub- 
linear expectation on (ri,Lg.(f2)) still denoted by E[-]. We now consider the 
extension of conditional expectations. For each fixed t < T, the conditional 
G-expectation E[-|r2t] : Lip{flT) — >■ Lip{flt) is a continuous mapping under ||-||. 
Indeed, we have 

E[x\nt] - E[Y\nt] < E[x - Y\nt] < e[\x - Y\\nt], 

then 

\E[x\nt] - E[Y\nt]\ < E[\x - Ywrtt]. 

We thus obtain 

|]E[X|r!t]-]E[r|f]t]|| <\\X^Y\\. 

It follows that E[-|r2t] can be also extended as a continuous mapping 

E[-\nt]:Ll,{nT)^LU^t). 

If the above T is not fixed, then we can obtain ]E[-|rit] : ^^(i^) — > L]j(yit)- 

Remark 2.5 The above proposition also holds for X,Y ^ L]^(^l). But in (iv), 
77 G Lg.(nt) should he bounded, since X,Y E L'q{^1) does not imply X - Y £ 

In particular, we have the following independence: 

E[x\nt] = E[x], yx e Ll;{n*). 



We give the following definition similar to the classical one: 
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Definition 2.6 An n- dimensional random vector Y e (^^(rj))" is said to be 
independent from Qt for some given t if for each ip e Cf,.Lip(M") we have 

E[^(r)|a] = iE[^(r)]. 

Remark 2.7 Just as in the classical situation, the increments of G-Brownian 
motion {Bt+s — Bt)s>o are independent from Qf 

The following property is very useful. 

Proposition 2.8 Let X,Y e i^(f^) be such that E[Y\nt] = -E[-Y\nt], for 
some t G [0,r]. Then we have 

E[X + Y\nt] = E[X\nt] + E[Y\flt]. 

In particular, ifE[Y\nt] = EG[->"|r2t] 0, then E[X + Y\nt] = E[X\nt]- 

Proof. This follows from the following two inequalities: 

E[x + Y\nt] < E[x\nt] + E[Y\nt], 

E[x + Y\nt] > E[x\nt] - E[-Y\nt] = E[x\nt] + iE[y|a]. 

□ 

Example 2.9 For each fixed a gM'*, s < t, we have 

E[Bt - Bfm = 0, Eh(i?f - Bf)\ns] = 0, 
E[{Bt - BfriQs] = al^At - s), E[-(i3f - Bf^in,] = -a^^^t s), 
E[{Bt - Bf)'\ns] = 3ai^rit - sf, E[-(Sr - Btfin,] = -Sat^^A* ~ ^f, 

where a'^^T — 2G(aa'^) and cr^^g^r = — 2G(— aa^). 

Example 2.10 For each a gR'^, neN, 0<t<T, Xe L}.{nt) and ip G 
C'i.Lip(K), we have 

t[Xip{B^ - Bf)\nt] = X+E[(p{B^ - Bf)\nt] + X-E[-ip{B^ - )|f7t] 
= X+t[p!{B^ - Bf)] + X-E[-ip{B^ - Bf)]. 

In particular, we have 

E[XiB^ - Bf)\nt] = X+E[{B^ - Bf)] + X-E[-{B^ - )] = 0. 
This, together with Proposition \2.8\. yields 

E[Y + X{B^ - Bf)\nt] - E[Y\nt], Y G Ll^ifl). 
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We also have 

E[X{B^ - B^f\flt] = X+]E[(Bf - Bf )2] + X-]E[-(Bf - B^f] 



E[X(BJ - Bff''-^\flt] = X+E[{B^ - )^"-^] + X-E[-{B^ - Bff''-'^] 
= |X|E[(B»_,)2n-i]. 

ExEimple 2.11 Since 

n2Bt{B^ - B^)\n,] = n-2Bf{Bf - Bf)\n,] = 0, 

we have 

EliBtf - {BfflQs] = E[(i?r - + Bff - {Bffm 

= E[{Bf - Bf)' + 2{Bf - Bf)Bf\n,] 

Exercise 2.12 Show that if X e Lip{nT) and E[X] = -E[-X], then E[X] 

/ \ where P is a Wiener measure on fi. 

Exercise 2.13 For each s,t>Q, we set Bf := Bt+s — Bg. Let rj = {Vij)^^^! € 
L^(f^a;§(d)). Prove that 

E[{7^BlBI)\ng] = 2G{7j)t. 



§3 Ito's Integral with G— Brownian Motion 

Definition 3.1 For T e M+, a partition ttt of [0,T] is a finite ordered subset 
""T = {iojiij • • • ,tN} such that = to <ti < ■ ■ ■ < tN = T. 

^(ttt) := max{|fi+i - tj| : i = 0, 1, • • • , A'' - 1}. 

We use TT^ = {tQ , t^ ,■■ ■ , t%} to denote a sequence of partitions of [0, T] such 
that limjv-).oo /^(""^ ) = 0. 

Let p > 1 be fixed. We consider the following type of simple processes: for 
a given partition ttt = {to, ■ ■ ■ , tN} of [0, T] we set 

N-l 

= X! ?fc('^)i[t<=,t<=+i)W, 

fe=0 

where ^fc G Lg,(ritj^), fc = 0, 1, 2, • • • , AT — 1 are given. The collection of these 
processes is denoted by Mg'°(0,r). 
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Definition 3.2 For an rj e Afg'°(0,T) with r]t{uj) = ELo ^kii^)h,,t,+,){t), 
the related Bochner integral is 



T N-l 

/ 'nt{uj)dt := ^ ^fe(w)(ife+i - tk)- 
For each j? e Mg'°(0, T), we set 

EtN := / 7?tdf] = ^E[^ ^kioj)itk+i - tk)]. 



It is easy to check that E^ : Mq°(0, T) — >• M forms a subhnear expectation. We 
then can introduce a natural norm ||f]l|Mg(o t); under which, Mq°(0,T) can be 
extended to M^(0,T) which is a Banach space. 

Definition 3.3 For each p > 1, we denote by Mq{0,T) the completion of 
MGf'°(0,T) under the norm 



lhllMS(o,T) ■■= IVtfdt] 



i/p 



It is clear that Mg(0,T) D M^(0,T) for 1 < p < g. We also use Mg(0,r;M") 
for all n-dimensional stochastic processes r]t = {vh'" ^Vt)-: t>Q with r^l G 
M^{0,T),i = l,2,--- ,n. 

We now give the definition of Ito's integral. For simplicity, we first introduce 
Ito's integral with respect to 1-dimensional G-Brownian motion. 

Let {Bt)t>o be a 1-dimensional G-Brownian motion with G(a) = ^{a'^a'^ — 
^a~), where < o; < ct < oo. 

Definition 3.4 For each r] € Mg'°(0, T) of the form, 

N-l 

vt{(^) = ^3(^)hi,tj+i){t), 

3=0 

we define 

„T N-l 

Lemma 3.5 The mapping I : Mq'^{0,T) — )■ LQ(f2T) is a continuous linear 
mapping and thus can be continuously extended to I : Mq{0,T) ^^(flT'). 
We have 

i-T 

VtdBt] = 0, (3.4) 



E[/ 

Jo 



E[(/ r]tdBt)^] < a^E[ [ rj^dt]. (3.5) 

^0 ^0 
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Proof. From Example 12. lOi for each j, 
We have 

E[ [ TjtdBt] - E[ / " ' rjtdBt + ^N-i{Btr, - J] 
Jo Jo 

Jo 

= E[ T^tdBt]. 
Jo 

Then we can repeat this procedure to obtain p.4p . 

We now give the proof of (|3.5p . Firstly, from Example I2.10[ we have 

E[(^ r,tdBtf]^t[(^jJ 'vtdBt+^N^Bt^ -Bt^_,ij ] 



tN-l \2 

2 I r> D \2 



= E[^ ritdBt^ +^2^_^^Bt,-Bt,_,) 
+ 2 (^^ " ' ritdBt^ ^N-i{Bt^ - St„_ J] 



N-l 



= ...=E[Y,inBu^,^Bun 



i=0 

Then, for each i = 0,1, ■ ■ ■ , N — 1, we have 



E[CfiBt^^,-BtJ^-a'^fiU+,-U)] 

1:2 I'D T3 \2 -=2c2f 



=E[E[e/(i?,,^, - Bty^a'^f{U+i - t,)\nt,]] 
=E[a''^UU+i - U) - a^i1{U+i - U)] = 0. 

Finally, we have 

E[( f rudB.f] = - Bu)'] 

JV-1 AT-l JV-1 

<E[^ ^f{Bt^^, - Buf - ^'^'(^^+1 - t^)\+nY. ^'^'(^'+1 - ^0] 

i=0 1=0 i=0 

JV-1 N-l 



< mfiBu+r - Buf - - m + E[E ^^i'tiu+i - u 

N-l T 

=E[E a^e^{U+l - U)] - / ,7,2 di]. 
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□ 

Definition 3.6 We define, for a fixed rj € Mq(0,T), the stochastic integral 



[ i^tdBt I{ri). 
Jo 



It is clear that dM]) and (jM]) still hold for r] e Mg(0,T). 

We list some main properties of Ito's integral of C?-Brownian motion. We 
denote, for some < s < t < T, 



Js Jo 



Proposition 3.7 Let 'n,9 e M^{0,T) and let < s < r < t < T . Then we 
have 

(i) ijudBu = £ TjudBu + £ TjudBu- 

(ii) J^{ariu + 9u)dBu — a rjudBu + /J O^dBu, if a is hounded and in Lq{VIs). 

(iii) t[X + r]udBup,] = t[X\n,] for X e L}.{n). 

We now consider the multi-dimensional case. Let G(-) : §((i) — > R be a 
given monotonic and sublinear function and let {Bt)t>o be a d-dimensional G- 
Brownian motion. For each fixed a gR'', we still use Bf := {a,Bt). Then 
(i?f)t>o is a 1-dimensional Ga-Brownian motion with Ga{a) — ^{cr^g^rO:^ — 
tQ!~), where cr^ t = 2G(aa^) and cP' t = — 2G(— aa-^). Similar to 1- 

— 3.3. ' 33 ^ ' — 33 ^ ' 

dimensional case, we can define Ito's integral by 



/ 77tdBf, for77e Af2(0,T). 
Jo 



m 

We stiU have, for each ^ e M^(0, T), 

Jo 

E[(/ TjtdBf)'] < al^rnl Vtdt]. 



Furthermore, Proposition 13.71 still holds for the integral with respect to Bf 
Exercise 3.8 Prove that, for a fixed rj e Mq{Q,T), 

a% [ riidt] < E[{ [ rjtdBtf] < a^E[ [ if^dt], 
Jo Jo Jo 

where ct^ = t[Bl] and ^ = -E[-Bf]. 

Exercise 3.9 Prove that, for each rj e Mq{0,T), we have 

nTivti^'dt] < rniijtHdt. 
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§4 Quadratic Variation Process of G— Brownian 
Motion 

We first consider the quadratic variation process of 1-dinicnsional G-Brownian 

motion {Bt)t>o witli Bi = N{{0} x [ct^ct^]). Let 7rf, iV = 1, 2, • • • , be a 
sequence of partitions of [0,t]. We consider 



N-l 
3=0 

N-1 N~l 

E2B^N{B^N - B^n) + {B^N - B^N)'^. 

3=0 3=0 

As /i(7rf^) — >■ 0, the first term of the right side converges to 2 BgdBg in Lg(r2). 
The second term must be convergent. We denote its limit by (i?)^, i.e., 

(B), lim ^ - B,.f ^ B^ - 2 B^dB,. (4.6) 

By the above construction, {{B)^)t>o is an increasing process with (-B)q = 0. 
We cah it the quadratic variation process of the G-Brownian motion B. 
It characterizes the part of statistic uncertainty of G-Brownian motion. It is 
important to keep in mind that {B)f. is not a deterministic process unless ct=(t, 
i.e., when {Bt)t>o is a classical Brownian motion. In fact we have the following 
lemma. 

Lemma 4.1 For each 0<s<t<oo, we have 

n{B)^-{B)JQ,]=a^t~s), (4.7) 
E[-{{B)^-{B)J\n,]^-a'{t-s). (4.8) 



Proof. By the definition of (B) and Proposition 13. 71 (iii). 

E[{B)^ - (B)^ \n,] - t[Bf - - 2^ B^dB^lils] 

= E[B^ - B^.lfls] ^a^{t-s). 

The last step follows from Example 12.111 We then have ()4.7p . The equality 
(|4.8p can be proved analogously with the consideration of E[—{Bf — B'^)\ils] — 
-^it-s). □ 
A very interesting point of the quadratic variation process {B) is, just like 
the G-Brownian motion B itself, the increment (B) — (B)^ is independent 
from fig and identically distributed with {B)^. In fact we have 
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Lemma 4.2 For each fixed s,t > 0, (-B)^^^ — {B)^ is identically distributed with 
(-B)j and independent from ilg. 

Proof. The results follow directly from 

- (B), = - 2 / ^ BrdBr - - 2 / BrdBr] 

Jo Jo 

fS+t 

= {B,+t - B,f - 2 / {Br- Bs)d{Br - B,) 

J s 

where (i?*) is the quadratic variation process of the G-Brownian motion B^ = 
Bs+t -B„t>0. □ 
We now define the integral of a process ri G Mq{0,T) with respect to {B). 
We first define a mapping: 

T N-1 



Lemma 4.3 For each r] G M^°{0,T), 

E[\Qo,Tm < 'J'n r \Vt\dt]. (4.9) 

JO 

Thus Qo,T ■■ M^'°(0,T) ^L1;(1]t) is a continuous linear mapping. Conse- 
quently, Qo,T co-n be uniquely extended to Mq{0,T). We still denote this map- 
ping by 

[ r,td{B)^:=Qo,T{^) for r, e M^{0,T). 
Jo 

We still have 

E[\ r 7itd{B\\]<a^E[r \nt\dt] forr]GM^{0,T). (4.10) 
Jo Jo 



Proof. Firstly, for each j = 1, • • • , A'' — 1, we have 

mAmt,^^ - {B)t,)-a%\{t,+i - 1,)] 

=E[mMB)t,+^ - {B)t,)\nt^]-a^mtj+i-tj)] 
=mj\a^tj+, - tj) - a'lmj+i - tj)] = 0- 
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Then (I4.9P can be checked as foUows: 



j=0 j=0 

j=o j=o 



□ 

Proposition 4.4 Let < s < t, ^ e L'^i^s), X €^^(1^). T/ien 

IE[X + - )] = E[X + i{Bt - Bsf] 

= E[X + im,-{B)^)]. 

Proof. By (|4.6p and Proposition 13.71 (iii) , we have 

E[X + - Bl)] = E[X + e((B>t - {B)^ + 2 / B^dB^)] 

J S 

= E[X + a{B)t-{B)^)]. 

We also have 

E[X + ^B^ - B^)] = t[X + ^{{Bt - Bsf + 2{Bt - Bs)Bs)] 
^E[X + ^{Bt-Bsf]. 

□ 

We have the following isometry. 
Proposition 4.5 Let 77 e Af|(0,T). Then 

E[( f ij,dB,f]=t[f 7f,d{B)^]. (4.11) 
Proof. We first consider 77 G Mq'^{0,T) of the form 

j=0 
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and then rjtdBt = J2f=o^ ^ji^tj+i — Bt-). From Proposition I3.7| we get 

t[X + 2ej(Bt,+, - Bt^)UBu+^ - Bt,)] = nX] for X e Ll;{n), i ^ j. 
Thus 

„T N-1 N-1 

E[( / T^tdBt)'] = IE[(^ ^,{Bt^^, - Bt^)f] = IE[ ^ e,{Bt^^^ - Bt,)']- 
From this and Proposition 14.41 it follows that 

E[(^ ^jtdBtf]^E[J2^]i{B),^^^-{B),^)]=nl v^d{B),]. 



Thus (P1T|) holds for 77 € M^''{0,T). We can continuously extend the above 
equahty to the case 77 e M^{0, T) and get (|4lT|l . □ 
We now consider the multi-dimensional case. Let {Bt)t>a be a d-dimensional 
G-Brownian motion. For each fixed a SR'^, (_Bf )t>o is a 1-dimensional Ga.- 
Brownian motion. Similar to 1-dimensional case, we can define 

N-l 

(B^)t , 1™ E (^t\, - ^r«)' = (Btr - 2 / BtdBt, 

where {B'^) is called the quadratic variation process of B^. The above 
results also hold for (B^). In particular, 

E[| f md{B^)t |] < ^Lrti f \^Mt] for 7? e Ml (0,T) 
Jo Jo 

E[(/ i^tdBtf]^t[[ iid{B^)^] forr;eM2(0,r). 
JO Jo 

Let a = (fli, • • • , a^)-^ and a = (ai, • • • , Od)-^ be two given vectors in M''. We 
then have their quadratic variation processes of {B^) and {B^). We can define 
their mutual variation process by 

(5^5')t ■■-\[{B- + B~-)t-{B--B^)t\ 

= \\{B^^~^)t-{B''-~^)t]- 

Since {B""-^) = {B^-"") = {-B''-^), we see that {B'',B^)^ = {B^,B'^)^. In 
particular, we have {B'^,B^) = (B^). Let tt^, = 1, 2, • • • , be a sequence of 
partitions of [0,t]. We observe that 

N-l N-l 

y: (Bt^^^ - Br«)(5f^, - ^f^) = i E - Btrn' - {Bt:. - Btrn 

k=0 k=a 
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Thus as /u(7r^) we have 

N-l 

Urn Y{B% -B%){B% -B%) = (B'',B^)^ 

fe=0 

We also have 



= ][{Bf+y-2f Bf+^dBf+^ - {B^-^'f + 2 f Bf-'^dBf-^] 
4 Jo Jo 

= BfBf - f BfdBf - f BfdBf. 
Jo Jo 

Now for each rj € Mq{0, T), we can consistently define 

£ vtd <i?^ B-\ = \ £ {B^+-\ - \ £ n,d (B'^-^ . 

Lemma 4.6 Let rj^ e M^'°(0, T), N = 1,2,- ■ ■ , be of the form 

N-l 
fe=0 

with niTT^) — and t]^ ^ r] in Mq{0,T), as N ^ oo. Then we have the 
following convergence in ^^.(fiT)- 

fc=o 

Proof. Since 

(B^,B^) ti — (B^ , B^) = {Bf-N — B%){B% — Bf-w) 



'-k+l _ f'k + l 



N k lu.N 

k '■k 



{Bf - B% )dBf - (Bf - B% )dBl 



we only need to prove 

Af-l 



nY.^^k f{ I ^'\b! - B%)dB^f] ^ 0. 



fe=0 

For each fc = 1, • • • , A'' — 1, we have 

Jtk 

{B- - B-.JdB^f\^,.] - Ci^mtk+i - t^f 

<nc{i^f{ti, - t^r - cient^^, - t^f] = o, 



§5 The Distribution of {B) 



51 



where C = (tI^t^-t/^. 
Thus we have 

N—l 
fe=0 

<E[x; (cf )'[( {Bf - Bt^)dBfr - c{t^^, - 

fe=0 ■'^k 

fe=0 

fe=0 

fe=0 "^0 

As — >■ 0, the proof is complete. □ 

Exercise 4.7 Let Bt be a 1-dimensional G-Brownian motion andip be a bounded 
and Lipschitz function on R. Show that 

^lim E[| J^^{B,.)[{B,.^^ - B,.f - ((B),. ^ - {B),.)]\] = 0, 

fe=0 

where if = kT/N, k = 0,2,- ■ ■ ,N - 1. 

Exercise 4.8 Prove that, for a fixed rj e Mq{0,T), 

a^t[[ \nt\dt]<E[[ \r]t\d{B)t]<a^E[[ \r,t\dt], 
Jo Jo Jo 

where u'^ = E[Bl] and ^ = -E{-B\\. 



§5 The Distribution of (S) 

In this section, wc first consider the 1-dimensional G-Brownian motion (-Bt)j>o 

with =iV({0} X [£2,^2]). 

The quadratic variation process (-B) of G-Brownian motion B is a very in- 
teresting process. We have seen that the G-Brownian motion i? is a typical 
process with variance uncertainty but without mean-uncertainty. In fact, {B) is 
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concentrated all uncertainty of the G-Brownian motion B. Moreover, (B) itself 
is a typical process with mean- uncertainty. This fact will be applied to measure 
the mean- uncertainty of risk positions. 

Lemma 5.1 We have 

t[{B)l] < lOaH^. (5.12) 

Proof. Indeed, 

E[{Bf,]^E[{Bt^-2 f BudBuf] 
Jo 

< 2E[Bf] + 8E[{ [ BudBuf] 

Jo 

<6aH^+8a^E[f Bjdu] 
Jo 

< 6aH'^ + 8ct^ / E[Bu'^]du 

Jo 

□ 

Proposition 5.2 Let {bt)t>o be a process on a sublinear expectation space (f2, H, E) 
such that 

(i) bo = 0; 

(ii) For each t, s > 0, bt+s — bt is identically distributed with bg and independent 
from {btj^,bt2, ■ ■ ■ , ) for each n G N and < i i , • • • ,tn < t; 

(iii) limt4oE[62]t-i = 0. 

Then bt is A^([/it, /Ji] x {0}) -distributed with JI = E\bi\ and fi = — E[— 

Proof. We first prove that 

]E[6t]=p and -E{-bt\^ ^d. 

We set if{t) := E[bt]. Then ip{0) — and limt_i_o (p{t) ~0. Since for each i, s > 0, 

^{t + s) = E[bt+s] = E[{bt+s - b,) + bs] 
^ip{t) + ip{s). 

Thus ip(t) is linear and uniformly continuous in t, which means that E[bt] = Jit. 
Similarly — ]E[— fet] = /it. 

We now prove that bt is N^^tjJIt] x {0})-distributed. By Exercise 11.171 in 
Chaplin we just need to prove that for each fixed (p e Cfc.Lip(R), the function 

u{t, x) := E[ip{x + bt)], {t, x) € [0, oo) x R 
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is the viscosity solution of the following parabolic PDE: 

dtu- g{dxu) = 0, u\t=o^'P (5.13) 

with g{a) = Jla'^ — fm^ . 

We first prove that u is Lipschitz in x and ^-Holder continuous in t. In fact, 
for each fixed t, u{t, •) eCh.Lip(E) since 

\E[^ix + bt)] ~ EMy + bt)]\ < EMx + bt) - ip{y + 6*)!] 

<C\x-y\. 

For each 5 £ [0, i], since bt ~ bg is independent from bg, we have 

u{t, x) = E[(y9(a; + 65 + {bt — bg)] 

= 'k\k[>p{v + {bt-bs))\y^,+b,i 

hence 

u{t,x) ^¥Jyu{t- 5,x + bg)]. (5.14) 

Thus 

\u{t, x) — u{t — 6,x)\ = \&(a{t — (5, X + bg) — u{t — 5, x)] \ 
< E[\u{t -6,x + bg)- u{t - 5, x)\] 
<E[C\bg\] <CiVS. 

To prove that u is a viscosity solution of the PDE ()5.13p . we fix a point 
{t,x) G (0,00) X R and let v G C^'^{[0,oo) x M) be such that v > u and 
v{t,x) — u{t,x). From (|5.14p . we have 

v{t, x) = E[u{t -S,x + bg)] < E[v{t ~S,x + bg)]. 

Therefore, by Taylor's expansion, 

0<E[v{t-S,x + bg)-v{t,x)] 
= E[v{t -5,x + bg)- v{t, X + bg) + {v{t, x + bg) ~ v{t, x))] 
= E[-dtv{t, x)5 + dMt, x)bs + Ig] 
< -dtv{t, x)6 + E[d^v{t, x)bg] + E[Ig] 
= -dtv{t, x)5 + g{dMt, x))S + E[Is] , 



where 



Ig=S [ [-dtv{t -I36,x + bg) + dtv{t, x)]d(3 
Jq 

+ bg [ [d^v{t,x + I3bg)- d^v{t,x)]dl3. 
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With the assumption that hmi4_oIE[fe^]i ^ = 0, we can check that 

limE[|/,|]ri-0, 

o-lO 

from which we get dtv(t, x) — g{dxv(t, a;)) < 0, hence u is a viscosity subsohition 
of ()5.13p . We can analogously prove that u is also a viscosity supersolution. It 
follows that bt is N{[iit,'pt] x {0})-distributed. The proof is complete. □ 
It is clear that {B) satisfies all the conditions in the Proposition 15. 2[ thus 
we immediately have 

Theorem 5.3 {B)^ is N{[q^t, a^t]x{0}) -distributed, i.e., for eachip G C/.Lip(M) 

sup ^{vt). (5.15) 

Corollary 5.4 For each < t < T < oo, we have 

a^{T -t) < {B)j, - <a^iT- t) in ^^(ll). 
Proof. It is a direct consequence of 



¥.[{{B)j,~ {B)^-a'{T -t)y]= sup {v - (t')+ {T - t) = Q 



and 



¥.[{{B)j,- {B)^-qf{T -t))-]= sup {v-qfy{T-t)^Q. 

□ 

Corollary 5.5 We have, for each t,s > 0, n €E N, 

m (s>,+, - {B)j^\n4 = n{B)-] = a^-r (5.16) 

and 

n-i {B)t+s - {B)sr\^s] - Eh (B)^] = -a^-r. (5.17) 

We now consider the multi-dimensional case. For notational simplicity, we 
denote by B^ := B^' the i-th coordinate of the G~-Brownian motion B, under a 
given orthonormal basis (ei, • • • , e^) of R'^. We denote 

m,),r-= {B\B=)^. 

Then {B)^, i > 0, is an S(c?)-valued process. Since 



E[{ABt,Bt)] = 2G{A)t for A e §(d), 
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we have 

d 



= E[J2 a^ABlBl - f BldBi - f BidBl)] 
Jo Jo 

d 

= E[ ^ a.jBlBl] = 2G{A)t for A e S(d), 



where (aij)fj^i = A. 

Now we set, for each ip e C;.Ljp(S(d)), 

v{t, X) t[ip{X + (t, X) e [0, oo) X §(d). 

Let S C S+(ii) be the bounded, convex and closed subset such that 

G(A) = i sup AeS(d). 

2 Bes 

Proposition 5.6 The function v solves the following first order PDE: 

dtv - 2G{Dv) = 0, v\t=o = 

where Dv — {dxijV)f j^i- We also have 

v{t,X) = sup(^(X + iA). 
Age 

Sketch of the Proof. We have 

v{t + 5,X)= t[ip{X + {B), + {B)^^, - {B),)] 
= t[v(t,X + {B),)]. 

The rest part of the proof is similar to the l-dimensional case. □ 

Corollary 5.7 We have 

e ffi := {t X 7 : 7 e S}, 
or equivalently, dt^i{B)^) = 0, where du(X) = inf{^(X -Y,X -Y) : Y eU}. 
Proof. Since 

E[dts((S)J] = sup dts(iA) =0, 

AGS 

it follows that dtY.{{B)^) = 0. □ 
Exercise 5.8 Complete the proof of Proposition \5.6l 
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§6 G-Ito's Formula 

In this section, we give Ito's formula for a "G-Ito process" X. For simplicity, 
we first consider the case of the function $ is sufficiently regular. 

Lemma 6.1 Let $ e C^{W) with d^:^<P, dl^^.^ € C6.Lip(M") for fJ,,i^ = 
1, • • • , n. Let s <E [0, T] be fixed and let X = {X^, • • • , X")^ be an n-dimensional 
process on [s, T] of the form 

Xr = X: + a'^it -s) + 7j-'\{B\B^\ - {B\ B^)J + r{Bl - Bi), 

where, for v = 1, • • • , n, i, j = 1, • • • ,d, , rf^^ and (3'^^ are bounded elements 
in L^(p.g) and Xs = {Xl, • • • , X")^ is a given random vector in Lq{^Is). Then 
we have, in L'Q{flt), 

^Xt)-<^{Xs)= [ d,.^XuW^dBi+ f d,.^Xu)a''du (6.18) 

J s J s 

+ j\d..^{X^W^' + ]^dl,,.^X^)rr]d{B\Bi)^. 

Here we use the , i.e., the above repeated indices fi.,u, i and j imply the sum- 
mation. 

Proof. For each positive integer N, we set S = {t — s) /N and take the partition 
<t] = {i^>if>--- ,t%} = {s,s + S,--- ,s + NS = t}. 

We have 

$(Xt) - $(X,) = J2 M^t-J - *(^t~)] (6.19) 

fe=0 
fe=0 

+ J[92,,.$(X,.)(Xf - - X^) + %^]}, 

where 

r?f = [a2,,.$(x,«+0fc(x,« -x,«))-a2,,.$(x,«)](x^'^, -x;„)(x;w -x;^„) 

fc fc+l k k ^k+1 *^fc "^fc+l ''k 

with Ok e [0,1]. We have 

nvi^f] = E[|[a^,,.$(x,^ + OkiX,.^^ - X,.)) - dl^^.^X,.)] 
X (xi^N — x'^)(xi^iv — xf'jv)!^] 

Ifc + 1 ^k + 1 ^k ' ' 

<cE[\X,.^^-X,.\<^]<C[d^ + S% 
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where c is the Lipschitz constant of {d^^^i,^}^ and C is a constant inde- 
pendent of k. Thus 

N-l 

k=0 k=Q 

The rest terms in the summation of the right side of (|6.19p are + C/^ with 

N-l 

+ r {Bl. - Bl„ )] + \dl,.^MXt-, {B\. - B\. ) {Bl, - Bl^ )} 

^k + l '■k 2 fc + l k ^k+1 '■k 

and 

N-l 
k=0 

+ 2K(i^+i-i^)+,7^^^((i?\i3^)^„ ~{B\B^)^^WiBl. -Bl^)}. 

'-k + l '■k fc+1 k 

We observe that, for each u £ [tk ^tk+i)^ 

N-l 
k=0 

= E[|9,.$(X„)-9,.$(X,«)|2] 
<c^m\X^- X,n\^] <C\6 + 5^], 

'k 

where c is the Lipschitz constant of {dx^^^Y^^i and C is a constant independent 
of k. Thus X;f=V5:r-*(^t^)I[t«,t^^j(-) converges to d^u^X.) in Ml{Q,T). 
Similarly, E^ro'52,,„ci>(X,«)I[,«,,«^^)(.) converges to in Af2(o,r). 

From Lemma l476l as well as the definitions of the integrations of dt, dBt and 
the limit of in L^iflt) is just the right hand side of (|08)) . By the 
next Remark we also have ^ in L%{Vlt). We then have proved (|6T8l) . □ 

Remark 6.2 To prove in LQ{V,t), we use the following estimates: for 

SMc/i t/iat lim^^oo mCt^t ) = and t[Y.k^o " ^k )] < for all N = 

1, 2, • • • , we have E[| Y!^'^ (^f+i - if ^ anrf, /or any fixed a, a SR'', 

AT-l N-l 
fe=0 fc=0 

k=0 
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N-1 

E[i E ^k{{B'')t^^^ - - if )n 

fe=0 
fe=0 

as well as 

fe=0 

<CE[E [^"'^(if+i - tk)\Bt^^^ - sr«n 

fc=0 

Af-l 

fe=0 

and 

AT-l 

fe=0 
Af-l 

fe=0 

AT-l 

fe=o 

□ 

We now consider a general form of G-Ito's formula. Consider 

Jo Jo Jo 

Proposition 6.3 Let $ G C^{W') with d^.<i>, dl^^.^ e Cb.iip(IR") /or /x,!/ = 
1, • • • , n. -Lef a", Z?''^ and if^^ , = 1, - ■ ■ ,n, i,j = 1, - ■ ■ ,d be bounded processes 
in Mq{0,T). Then for each t>0 we have, in ^^'(^t) 

<^{Xt) - <^{Xs) = I d,.^{XuWJdBi+ I (6.20) 

J s J s 
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Proof. We first consider the case where a, ij and /3 are step processes of the 
form 

N-l 

= &(^^)i[t,,t,+i)(*)- 

fc=0 

From the above lemma, it is clear that (|6.20p holds true. Now let 

Jo Jo Jo 

where , rj^ and are uniformly bounded step processes that converge to 
a, rj and j3 in Mq{Q,T) as N ^ oo, respectively. From Lemma l6. 11 



) - <i>(xf ) = / d^.m^)(3:^^''dBi+ / d,.<i>{x^)a:^''du (6.21) 

J s J s 



Since 



E[\xr-xn'] 

Jo 

where C is a constant independent of A^, we can prove that, in A/q(0,T), 

We then can pass to limit as X — > oo in both sides of (|6.2ip to get (|6.20p . □ 
In order to consider the general $, we first prove a useful inequality. 
For the G-expectation E, we have the following representation (see ChapfVTI: 

t[X] = sup Ep[X] for X e L^(r2), (6.22) 

Pev 

where is a weakly compact family of probability measures on {Q.,B{^1)). 

Proposition 6.4 Let /3 e Afg(0,r) with p > 2 and let a e be fixed. Then 
we have /3tdBf e Lq{D,t) o,nd 

E[| r PtdBt\P] < CpE[\ r pU{B'')tf\ (6.23) 
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Proof. It suffices to consider the case where /3 is a step process of the form 

N-l 



fc=0 

For each ^ e Lip{Vlt) with t E [0,T], we have 



From this we can easily get Ep[£^ (3sdBf] — for each P £ V, which imphes 
that (/o PsdBf )t 

go,T] is a P-martingale. Similarly we can prove that 

Mt:=if MBf)^ - f t € [0,T] 

Jo Jo 

is a P-martingale for each P £ V. By the Burkholder-Davis-Gundy inequalities, 
we have 

Ep[\ r MBfn<C,Ep[\ r l3^d{B'^)t\P/']<CpE[\ C (3^^d{B^)t\P/\ 

Jo JQ Jo 

where Cp is a universal constant independent of P. Thus we get (I6.23p . □ 
We now give the general G-Ito's formula. 

Theorem 6.5 Let $ fee a -function on W such that d^^^,.^ satisfy polyno- 
mial growth condition for fi^v — 1, ■ ■ ■ ,n. Let , fi^^ and r]'^^^ , v = 1^ ■ ■ ■ ,n, 
i,j = l,---,dbe bounded processes in Mq{0, T). Then for each t > we have 
in Ll,{Vlt) 

^{Xt) - ^{Xs) = ^ d,.^{X^)Pl^dBi + ^ (6.24) 
[a,.$(X„)C^' + \dl,^,MXuW[i''J]d{B\B')^ . 

Proof. By the assumptions on $, we can choose a sequence of functions $Af G 
C^(M") such that 



|<i>^(a;)-<i>(a;)| + |a,.<i>jv(x)-a,.$(a;)| + |a2^^„$^(x)-a^,,.<i>(x)| < ^(l+|x|'=), 

where Ci and k are positive constants independent of N . Obviously, satisfies 
the conditions in Proposition [631 therefore, 

^N{Xt) - ^n{Xs) = ^ d,.^N{Xu)p''JdBi + ^ d^.^NiXuWu (6.25) 
+ f\d^.<^N{Xu)Vu' + \dl.,.<fN(XXP:^]d{B\B^)^ 
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For each fixed T > 0, by Proposition 16.41 there exists a constant C2 such that 

E[|Xt|2'=] < C2 for t £ [0,T]. 
Thus we can prove that $Ar(Xt) ^ <I>(Xt) in L'^Q{^t) and in Af^(0,r), 

d^.<^>N{x.)r ^d.,.<^>{x.)r- 

We then can pass to limit as — > 00 in both sides of (|6.25p to get (|6.24l) . □ 

Corollary 6.6 Let ^ be a polynomial and a, a^eR'* be fixed for v = 1, ■ ■ • ,n. 

Then we have 

where Xt = {Bf , ■ ■ • , Bf . In particular, we have, for k = 2,3, ■ ■ ■ , 

[Bff = k j\B^f-'dB^+^-^^^^ j\Bff-^d{B^),. 

If E becomes a hnear expectation, then the above G-Ito's formula is the 
classical one. 

§7 Generalized G-Brownian Motion 

Let G : M"^ X S((i) — )• R be a given continuous sublinear function monotonic in 
A S §(d). Then by Theorem 12. ll in ChapU there exists a bounded, convex and 
closed subset S c R'' x S+ {d) such that 

G{p, A) = sup [\iT[AB] + {p, q}] for {p, A) e R'^ x §{d). 

By Chapter ini we know that there exists a pair of d-dimensional random vectors 
{X, Y) which is G-distributed. 

We now give the definition of the generalized G-Brownian motion. 

Definition 7.1 A d- dimensional process {Bt)t>Q on a sublinear expectation 
space {il,H,K) is called a generalized G-Brownian motion if the follow- 
ing properties are satisfied: 

(i) Boioj) = 0; 

(ii) For each t, s > 0, the increment Bt+s — Bt identically distributed with 
y^X + sY and is independent from {Bt^,Bt2, ■ ■ ■ ,Bt^), for each n G N and 
< ti < ■ ■ ■ < tn < t, where {X, Y) is G-distributed. 
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The following theorem gives a characterization of the generalized G-Brownian 
motion. 

Theorem 7.2 Let {Bt)t>o be a d- dimensional process defined on a sublinear 
expectation space (fij'HjE) such that 

(i) BoiLo)= 0; 

(ii) For each t, s > 0, Bt+s — Bt and Bg are identically distributed and, iJj+s — Bt 
is independent from {Bt^, Bt^ , • • • i -Bf „ ) , for each n € N and < ti < ■ ■ ■ < tn < 
t. 

(iii) limt;oE[|Bt|3]f-i =0. 

Then (i?t)/>o is a generalized G-Brownian motion with G{p, A) = lim^j^o Bi)-\- 

\{ABs,Bs)]5-^ for {p,A) x S{d). 

Proof. Wc first prove that lim^j^o -Ba) + ^{ABs, Bs)]S~^ exists. For each 
fixed {p,A) eWL"^ X S{d), we set 

fit) :=E[{p,Bt) + ^{ABt,Bt)]. 

Since 

\f{t + h)- f{t)\ < E[{\p\ + 2\A\\Bt\)\Bt+h - Bt\ + \A\\Bt+h - B^p] ^ 0, 
we get that f{t) is a continuous function. It is easy to prove that 
E[{q,Bt)]=E[{q,B,)]t for g e M'^. 

Thus for each t,s > 0, 

\f{t + s)-f{t)-f{s)\<CE[\Bt\]s, 

where C = |^|E[|Bi|]. By (iii), there exists a constant > such that 
E[|Bt|3] < f for t < ^0- Thus for each fixed i > and A'' e N such that 
Nt < 6o, we have 

\f{Nt)-Nf{t)\<^C{Nt)'/K 

From this and the continuity of /, it is easy to show that limtj^o f{t)t^^ exists. 
Thus we can get G{p,A) for each {p,A) s R'' x It is also easy to check 

that G is a continuous sublinear function monotonic in A S 

We only need to prove that, for each fixed y e Cb.upi^'^)-, the function 

u{t, x) := E[ip{x + Bt)], (t, x) e [0, oo) x R'^ 

is the viscosity solution of the following parabolic PDE: 

dtu - G{Du, D'^u) = 0, u\t=o = V- (7.26) 
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We first prove that u is Lipschitz in x and ^-Holder continuous in t. In fact, 
for each fixed t, u{t, •) &Cb.Lip{^'^) since 

\E[^ix + Bt)] - E[ip{y + Bt)]\ < E[|v?(x + Bt) - ip{y + Bt)\] 

<C\x-y\. 

For each 5 £ [0,i], since Bt — Bs is independent from Bg, 

u{t, x) = E[^{x + Bs + {Bt - Bs)] 

= E[E[ip{y + {Bt - Bs))]y=:.+B,]. 

Hence 

u{t,x)=E[u{t-S,x + Bs)]. (7.27) 

Thus 

\u{t,x) - u{t-S,x)\ = \E[u{t-S,x + Bs) -u{t-S,x)]\ 

< E[\u{t -S,x + Bs) - u{t - 6, x)\] 

< E[C\Bs\] < Cy/G{0,I) + lVS. 

To prove that u is a viscosity solution of (|7.26l) . we fix a {t,x) G (0,oo) x R'^ 
and let v € C^''^([0, cx)) x M!^) be such that v > u and v{t,x) — u{t,x). From 
(fr27|) . we have 

v{t, x) = E[u{t -S,x + Bs)] < E[v{t -S,x + Bs)]. 

Therefore, by Taylor's expansion, 

< E[v{t -S,x + Bs) - v{t, x)] 
= E[v{t -S,x + Bs) - v{t, X + Bs) + {v{t, x + Bs) - v{t, x))] 

= E[-dtv{t, x)5 + {Dv{t, x), Bs) + ^{D\{t, x)Bs,Bs) + Is] 
< -dtv{t, x)6 + E[{Dv{t, x), Bs) + ^{D^v{t, x)Bs,Bs)] + E[Is], 



where 



1 

--[dtv{t -I36,x + Bs) - dtv{t, x)]Sdp 







{{D^v{t, X + ajSBs) - D^v{t, x))Bs,Bs)ad(3da. 



Jo 



With the assumption (iii) we can check that lim^j^o = 0, from which 

we get dtv{t,x) — G{Dv{t, x), D^v{t, x)) < 0, hence m is a viscosity subsolution 
of (|7.26p . We can analogously prove that -u is a viscosity supersolution. Thus u 
is a viscosity solution and {Bt)t>o is a generalized G-Brownian motion. □ 



64 



Chap. Ill G-Brownian Motion and Ito's Integral 



In many situations we are interested in a generalized 2(i-dimensional Brow- 
nian motion {BtM)t>o such that t[Bt] = -t[-Bt] = and E[|6t|2]/i 0, as 
i J, 0. In this case B is in fact a G-Brownian motion defined on Definition 2.1 of 
Chapter [Til Moreover the process b satisfies properties of Proposition 5.2. We 
define u{t, x, y) = 'E[Lp{x + Bt,y + bt)]. By the above proposition it follows that 
u is the solution of the PDE 

dtu = G{DyU, Dl^u), u\t=o = ve C/.L,p(R2'^). 

where G is a sublinear function of {p, A) G M"*, defined by 

G{p,A) ■.^E[{p,bt) + {ABt,Bt)]. 

Here (•, •) = {■r)^^. 

§8 G-Brownian Motion under a Nonlinear Ex- 
pectation 

We can also define a G-Brownian motion on a nonlinear expectation space 

Definition 8.1 A d- dimensional process {Bt)t>Q on a nonlinear expectation 
space (rij'HjE) is called a (nonlinear) G-Brownian motion if the following 
properties are satisfied: 

(i) BoH =0; 

(ii) For each t, s > 0, the increment Bt+s — Bt identically distributed with Bs 
and is independent from {Bt^, Bt^, • • • , Bt^), for each n gN and < ti < ■ ■ ■ < 
tn < t; 

(iii) \ixiitionW]t-^ =Q. 

The following theorem gives a characterization of the nonlinear G-Brownian 
motion, and give us the generator G of our G-Brownian motion. 

Theorem 8.2 LetK be a nonlinear expectation andE be a sublinear expectation 
defined on (U,,'H). lefR be dominated by E, namely 

E[X]-E[y] <E[X-r], X,Y gU. 

Let {Bt,bt)t>o be a given M.^'^-valued G-Brownian motion on (ri,H,E) such 
that E[Bt] = E[-Bt] = and limt^o E[|6tp]/t = 0. Then, for each fixed ip G 
Cb. 

Lip(J^ )j function 

uit, X, y) E[ip(x + Bt,y + bt)], {t, x, y) G [0, cx)) x M^d 
is the viscosity solution of the following parabolic PDE: 

dtU-G{DyU,Dlu)^{), u\t=Q^ip. (8.28) 
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where 

G(p, A) -E[(p,6i> + {p,A) e R'^ X §(d). 

Remark 8.3 Let G{p,A) := E[{p,bi) + ^{ABi,Bi)]. Then the function G is 
dominated by the sublinear function G in the following sense: 

G(p,A)-G(p',A')<G(p-p',A-A'), (y,A')eM^'xS(d). (8.29) 

Proof of Theorem We set 

/(<) = fAAt) ■■= n{P,bt) + ^{ABt,Bt)], t > 0. 

Since 

\f{t + h)- f{t)\ < E[(\p\ + 2\A\\Bt\)\Bt+h ~Bt\ + \A\\Bt+h - StP] ^ 0, 

we get that f{t) is a continuous function. Since E[Bt] = ]E[— i?t] = 0, it follows 
from Proposition O that E[X + {p,Bt)] = E[X] for each X e H a.nd p e W^. 
Thus 

fit + h) = E[{p,bt+k-bt) + {p,bt) 

+ ^{ABt+h - Bt,Bt+h ~Bt) + ^ {ABt,Bt}] 

= E[{p, bh) + ^{ABh, Bh)] + E[{p, bt) + \{ABu Bt)] 
= f{t) + f{h). 

It then follows that f{t) = f{l)t = G{A,p)t. We now prove that the function 
u is Lipschitz in x and uniformly continuous in t. In fact, for each fixed 
u{t,-) eCb.Lzpi^'^) since 

\E[^{x + Buv + bt)] - Ey{x' + Bu y' + h)]] 

< E[\ip{x + Buy + bt) - ^{x' + Bt,y' + 601] < C{\x ~ x'\ + |y - y'\). 
For each (5 G [0, t], since (Bt ~ Bs,bt — bg) is independent from {Bs, bs), 
u(t, x) = E[ip(x + Bs + {Bt - Bs), y + bs + {bt - bs)] 

= E[E[v3(x + {Bt -Bs),y + {bt ~ bs))]s=x+Bs.y=y+bA- 

Hence 

u{t,x) ^E[u{t - S,x + Bs,y + bs)]. (8.30) 

Thus 

\iL{t,x,y) - u{t - 5,x,y) \ = \E[u{t - 5,x + Bs,y + bs) - u{t - 5, x, y)]\ 

< E[\u{t -6,x + Bs,y + bs) - u{t - 6, x, y)\] 
<E[C{\Bs\ + \bs\)]. 



66 



Chap. Ill G-Brownian Motion and Ito's Integral 



It follows from (iii) of Definition 18. Il that u{t,x,y) is continuous in t uniformly 
in {t,x) G [0,oo) X R2<i_ 

To prove that u is a viscosity solution of (|8.28p , we fix a (t, x,y) E (0, oo) x M^"* 
and let v € C^''^([0, oo) x R^'*) be such that v > u and v{t,x,y) = u{t,x,y). 
From (|8.30p . we have 

v{t,x,y) ^E[u{t-S,x + Bs,y + bs)] < E[v{t - d,x + Bs^y + bs)]. 

Therefore, by Taylor's expansion, 

< E[u(< - S,x + Bs^y + bs) - v{t, x)] 
= E[v{t - 5,x + Bs,y + bs) - v{t, x + Bs,y + bs) 
+ {v{t,x + Bs,y + bs) -v{t,x,y)] 

= E[-dtv{t, X, y)5 + {Dyv{t, x, y)M) + {dxv{t, x, y), Bs) + ^{Dl^v{t, x, y)Bs,Bs} + Is] 
< -dtvit, X, y)6 + E[{Dyvit, x, y),bs) + ^{Dl^t, x, y)Bs, Bs)] + E[Is], 



where 



Is = / -[dtv{t-S-f,x + Bs,y + bs)-dtv{t,x,y)]Sdj 

JQ 

+ / {dyv{t,x + -fBs,y + '^bs) - dyv{t,x,y),bs) d'^ 
Jo 

+ / {dxv{t,x,y + jbs) - dxv{t,x,y),Bs} 







1 /.I 

X + ajBs, y + ^bs) - Dl^v{t, x, y))Bs,Bs)^d^dc 



^0 



With the assumption (iii) we can check that lim^j^o IE[|/5|](5^^ — 0, from which 
we get dtv{t,x) — G{Dv{t,x),D^v{t,x)) < 0, hence m is a viscosity subsolution 
of (I8.28p . We can analogously prove that m is a viscosity supersolution. Thus u 
is a viscosity solution. □ 

§9 Construction of G-Brownian Motions under 
Nonlinear Expectation 

Let G(-) : M"* x S{d) — > M be a given sublinear function monotonic on A e S{d) 
and G{-) : M'' x §{d) M be a given function dominated by G in the sense of 
(|8.29l) . The construction of a R^''-dimensional G-Brownian motion {Bt,bt)t>o 
under a nonlinear expectation E, dominated by a sublinear expectation E is 
based on a similar approach introduced in Section 2. In fact we will see that 
by our construction (i?f,fef)(>o is also a G-Brownian motion of the sublinear 
expectation E. 
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We denote hy Q — Cq'^{R'^) the space of all R^'^-valued continuous paths 
(i^t)tGR+- For each fixed T G [0,cx)), we set Qt ■— {w.at : € We will 
consider the canonical process {Bt,bt){uj) — ujt, t £ [0,oo), for uj £ ft. We also 
follow section 2 to introduce the spaces of random variables Lip{Q.T) and Lip{^) 
so that to define E and E on (fi, Lip{Q)). 

To this purpose we first construct a sequence of c?-dimensional random vec- 
tors {Xi,r]i)f^^ on a sublinear expectation space (ri,?^,E) such that {Xi,r]i) is 
G-distributed and (X^+i, 77^+1) is independent from ((Xi,r]i), ■ ■ ■ ,{Xi,'r]i)) for 
each z = 1, 2, • • • . By the definition of G-distribution the function 

u(t,x,y) ■.= E[(p{x + ViXi,y + trii)], t>0, x,y eR'^ 

is the viscosity solution of the following parabolic PDE, which is the same as 
equation p.6p in Chap. II. 

dtu - G{DyU, Dl^u) = 0, u\t=Q = (pe CupiM.^'^). 

We also consider the PDE (for the existence, uniqueness, comparison and dom- 
ination properties, see Theorem 12.61 in Appendix C). 

dtu ~ G{DyU, Dl^u) = 0, u\t=Q = ipe GL.p(K2d), 

and denote by Pt [</3] (x, y) — u(t, x, y). Since G is dominated by G, it follows from 
the domination theorem of viscosity solutions, i.e., Theorem l3.5l in Appendix C, 
that, for each (p,tpe Cb.Lipi^^'''), 

Pt [if] ix,y)- Pt [^] ix,y)< niip -^){x + VtXi ,y + tr^i)]. 

We now introduce a sublinear expectation E and a nonlinear E defined on 
Lipip.) via the following procedure: for each X £ Lip{il) with 

X = f(Bt^ -Bto,bt, -bt„,--- ,Bt, - Btg,bt„ -h^_^) 
for (p £ Ci,Lip{R?'^^'^) and = to < < • ■ • < < 00, we set 

t[p{Bt,- Bt^^,bt^-bt,r■■ ,Bt^- Bt^_,Mr.-K-i)] 
:= E[(^(Vii - taXi, {ti - to)T]i, • • • , ^<„ - t„_iX„, (t„ - i„_i)77„)]. 

and 

E[^{Bt,- Bt„bt,-bt,r-- ,St„ -6t_J] =(^„(0,0) 

where (pn £ Cb.Lipi^'^'^) is defined iteratively through 

'Pl{xi,yi,- ■ ■ ,Xn-l,yn-l) = ^'t„-t,._i [</?! (2:^1, 2/1, ■■ • , Xn~l , Vn-l , ■)]{0, 0) , 



fn-i{xi,yi) 
Vn{xi,yi) 



= Pt2-ttVPn^2{xi,yi,-)]{Q,Q), 

= Pt2-ttW7i.-i{-)]{xi,yi). 
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The related conditional expectation of X =ip{Bt^ — Bt„,bt^ — bt„, - ■ ■ ,Bt^ — 
Bt„_i,bt^ — bt^_-^) under Qtj is defined by 

E[X\nt,] ^n^iBt, ~ Bt,,bt, - bt„, - ■ ■ ,Bt„~Bt^_,,bt,^-bt„_,) \nt^] (9.31) 
:= ■0(5*1 - Bt„,bt, - bt„, ■■■ ,Bt^ - Bt.^,,bt^ - bt._,), 

where 

■ip{xi,--- ,Xj) = E[(p{xi,--- , Xj, v^^j+i - tjXj+i, (ii-io)?yj+i, ■ • • , ^yt~^t^Xn, {ti-to)r]n)]- 
Similarly 

E[X\nt,]^^n-jiBt,~Bt,,bt,-bt,r-- .Bt^-Bt^_,,bt^-bt^_,). 

It is easy to check that E[-] (rcsp. E) consistently defines a sublinear (resp. 
nonlinear) expectation and E[-] on {il, Lip{il.)). Moreover {Bt,bt)t>o is a G- 
Brownian motion under E and a G-Brownian motion under E. 

Proposition 9.1 We also list the properties o/E[-|rit] that hold for each X, Y ELip{il): 

(i) IIX>Y, thenE[X\nt] > E[Y\nt]. 

(ii) E[X + T]\nt]^E[X\nt]+ri, for each t > and rj eLip{nt). 

(iii) E[X\nt] - E[Y\qt] < E[X - Y\nt]. 

(iv) E[E[X|f]t]|f]s] = E[X\ntj\s], m particular, E^[X\Q.t]] = E[X]. 

(v) For each X S Lip(yi*), E[X\Q,t\ = E[X], where Lip{n^) is the linear space 
of random variables with the form 

n = l,2,--- , ^ e C,.Lip(M'*''"), h,--- ,tn,t„+i e [t,oo). 

Since E can be considered as a special nonlinear expectation of E dominated 
by its self, thus E[-|ilf] also satisfies above properties (i)-(v). Moreover 

Proposition 9.2 The conditional sublinear expectation E[-\Qt] satisfies (i)-(v). 
Moreover E[-\^t\ itself is sublinear, i.e., 

(vi) E[X\nt]-E[Y\^t]<E[X -Y\nt], . 

(vii) E[r)X\nt] = ri+E[X\Vlt] + ri-E[-X\Vlt] for each rj € L.piQt)- 

We now consider the completion of sublinear expectation space {ft, Lip{il),E). 

We denote by Lg,(51), p > 1, the completion of Lip{il) under the norm 
:= (E[|X|P])i/P. Similarly, we can define ^^(I^t), L^^i^tj,) and £^(f^*). 
It is clear that for each < t < T < oo, L^i^lt) ^ Lg(^t) C L^i^l). 

According to Sec|l]in Chap|ll ]£[•] can be continuously extended to {ft, L^(57)). 
Moreover, since E is dominated by E, thus by Definition l4.4l in Chap. I, (J7, Lq^H.), E) 
forms a sublinear expectation space and (51, LQ^fl), E) forms a nonlinear expec- 
tation space. 
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We now consider the extension of conditional G-expectation. For each fixed 
t < T, the conditional G-expectation E[-|ilt] : Liplflr) Lip{Vlt) is a continu- 
ous mapping under ||-||. Indeed, we have 

nxm - nY\^t] < nx - Y\nt] < nix - rp*], 

then 

\E[x\nt] - E[Y\nt]\ <K[\x -Y\\nt]. 

We thus obtain ^ ^ 

|E[X|r!t]-E[y|r!t]| < \\X-Y\\. 

It follows that E[-|r2t] can be also extended as a continuous mapping 

E[-\nt]:LUnT)^L},{nt). 

If the above T is not fixed, then we can obtain E[-|ri(] : ^^(ri) LQ{ilt). 

Remark 9.3 The above proposition also holds for X,Y LQ{fl). But in (iv), 
77 G L^ijlt) should be bounded, since X,Y Q LQ(fl) does not imply X ■ Y 

In particular, we have the following independence: 

E[x\nt] = E[x], yx e 

We give the following definition similar to the classical one: 

Definition 9.4 An n- dimensional random vector Y G (L^(r2))" is said to be 
independent from Qt for some given t if for each ip £ Cb.Lipi^"') we have 

E[^(y)|r!t] = E[^(r)]. 
Notes and Comments 

Bachelier (1900) proposed Brownian motion as a model for fluctuations of the 
stock market, Einstein (1905) (4^ used Brownian motion to give experimental 
confirmation of the atomic theory, and Wiener (1923) 121] gave a mathemati- 
cally rigorous construction of Brownian motion. Here we follow Kolmogorov's 
idea (1956) [t^ to construct G-Brownian motion by introducing infinite di- 
mensional function space and the corresponding family of infinite dimensional 



sublinear distributions, instead of linear distributions in 74 1. 

The notions of G-Brownian motion and the r elate d stochastic calculus of 
Ito's type were firstly introduced by Peng (2006) [Too'l for 1-dimensional case 
and then in (2008) ^104] for multi-dimensional situation. It is very interesting 
that Denis and Martini (2006) [s^ studied super-pricing of contingent claims 
under model uncertainty of volatility. They have introduced a norm on the 
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space of continuous paths il — C([0, T]) which corresponds to our Lg-norm and 
developed a stochastic integral. There is no notion of nonlinear expectation 
and the related nonlinear distribution, such as G-expectation, conditional G- 
expectation, the related G-normal distribution and the notion of independence 
in their paper. But on the other hand, powerful tools in capacity theory enable 
them to obtain pathwise results for random variables and stochastic processes 
through the language of "quasi-surely" (see e.g. Dellacherie (1972) [32], Del- 
lacherie and Meyer (1978 and 1982) [H, Feyel and de La Pradelle (1989) fisl) 
in place of "almost surely" in classical probability theory. 

The main motivations of G-Brownian motion were the pricing and risk mea- 
sures under volatility uncertainty in financial markets (see Avellaneda, Levy and 
Faras (1995) d and Lyons (1995) It was well-known that under volatil- 

ity uncertainty the corresponding uncertain probabilities are singular from each 
other. This causes a serious problem for the related path analysis to treat, 
e.g., path-dependent derivatives, under a classical probability space. Our G- 
Brownian motion provides a powerful tool to such type of problems. 

Our new Ito's calculus for G-Brownian motion is of course inspired from 
Ito's groundbreaking work since 1942 65| on stochastic integration, stochastic 
differential equations and stochastic calculus through interesting book s cited 



in Chapter IIVI Ito's formula given by Theorem 16.51 is from Peng [1Q0|, 104 1. 
Gao (2009) [5J] proved a more general Ito's formula for G-Brownian motion. 
An interesting problem is: can we get an Ito's formula in which the conditions 
correspond the classical one? Recently Li and Peng have solved this problem in 

Using nonlinear Markovian semigroup known as Nisio's semigroup (see Nisio 
(1976) tM]), Peng (2005) ^981] studied the processes with Markovian properties 
under a nonlinear expectation. 



Chapter IV 

G- martingales and Jensen's 
Inequality 

In this chapter, we introduce the notion of G-martingales and the related 
Jensen's inequality for a new type of G-convex functions. Essentially differ- 
ent from the classical situation, "M is a G-martingale" does not imply that 
"— M is a G-martingale" . 

§1 The Notion of G-martingales 

We now give the notion of G-martingales. 

Definition 1.1 A process (Mt)t>o is called a G— martingale (respectively, G- 
supermartingale, G-submartingale) if for each t G [0,00), Mt G ^^(ilt) 
and for each s G [0,t], we have 

t[Mt\fls]= Ms (respectively, < M^, > M^). 

Example 1.2 For each fixed X e Lq{CI), it is clear that {E[X\Clt])t>o is a 
G-martingale. 

Example 1.3 For each fixed a S R*^, it is easy to check that {Bf)t>o and 
{—Bf)t>Q are G-martingales. The process {{B^)t — (^1^Tt)t>o is a G-martingale 
since 

n{B^)t - oi^Tt\n,] = n{B'')s - ^Lrt + {{B^)t - {B^)s)m 

= {B'^)s-al^Tt + t[{B'^)t-{B-),\ 
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Similarly we can show that {—{{B^)t — '^'^gTt))t>o is a G-submartingale. The 
process ((Sf )^)t>o is a G-submartingale since 

E[{Btf\n,] = E[{Bfr + (Bf - Bf)' + 2Bf{Bf - Bf)\n,] 
= iBfr+E[{B^-Bf)'\ns] 
= {Bfr + al^At-s)>{Bfr. 

Similarly we can prove that ((Bf)^ — cr^a^i)t>o o-nd {{Bf)"^ — {B^)t)t>o o-re 
G-martingales. 

In general, we have the following important property. 

Proposition 1.4 LetM^ e R, = (v'^j^Li e Af2(o,r;R'') andri = {r]'^)f^j^i e 
M^(0,T;§(d)) be given and let 

= Mo + ^JBi + £ {B\ B^)^ - 2G{rju)du for t e [0, T]. 

Then M is a G -martingale. Here we still use the , i.e., the above repeated 
indices i and j imply the summation. 

Proof. Since E[f^ ipidBi\n,] ^ E[- ipidBHils] = 0, we only need to prove 
that 

Mt= / v'Jd{B\B^) - [ 2G{r]u)du forte[0,T] 
Jo Jo 

is a G-martingale. It suffices to consider the case where rj G Mg'°(0, T; S(d)), 
i.e., 

N-l 

We have, for s G [tN-i,tN], 

E[Mtm =Ms+ E[{nt,_,,{B), - (S) J - 2G(7?t„_ J(i - 

= M, + E[{A, {B)^ - (S) J]^=,,,_^ - 2G(7?t,_ J(t - s) 
= Ms. 

Then we can repeat this procedure backwardly to prove the result for s € 

[0,tN-l]. □ 

Corollary 1.5 Let rj S Mq(0,T). Then for each fixed a S M.'^, we have 

^-^■rn r \Vt\dt] < E[ r \vt\d{B-)t] < al,^E[ C \'nt\dt]. (1.1) 
Jo Jo Jo 
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Proof. For each ^ e M^{0,T), by the above proposition, we have 

E[ / ^td{B^)t - / 2Ga(6)di] = 0, 
Jo Jo 

where Ga(a) = 5(0-^^"^ " ^-aa^"~)- Letting ^ = |?7| and ^ = -\r]\, we get 

Jo Jo 

E[- f \m\d{B-)t + cr\^T f \vt\dt] = 0. 
Jo Jo 

From the sub-additivity of G-expectation, we can easily get the resuh. □ 

Remark 1.6 It is worth to mention that for a G-martingale M , in general, 
—M is not a G -martingale. But in Proposition ] when rj = 0, —M is still a 
G -martingale. 

Exercise 1.7 (a) Let {Mt)t>o be a G-supermartingale. Show that {^Mt)t>o "is 
a G-submartingale. 

(b) Find a G~submartingale {Mt)t>o such that {—Mt)t>o not a G-supermartingale. 

Exercise 1.8 (a) Let {Mt)t>o o,nd (A^t)t>o be two G-supermartingales. Prove 
that {Mt + Nt)t>o is a G-supermartingale. 

(b) Let {Mt)t>o Oif^d {—Mt)t>o be two G -martingales. For each G-submartingale 
(respectively, G-supermartingale) {Nt)t>o, prove that {Mt + Nt)t>o is a G- 
submartingale (respectively, G-supermartingale). 



§2 On G-martingale Representation Theorem 

How to give a G-martingale repr esentation theorem is still a largely open prob- 
lem. Xu and Zhang (2009) [122l| have obtained a martingale representation for 
a special 'symmetric' G-martingale process. A more general situation have been 
proved by Soner, Touzi and Zhang (preprint in private communications). Here 
we present the formulation of this G-martingale representation theorem under 
a very strong assumption. 

In this section, we consider the generator G : S(d) — > E satisfying the uni- 
formly elliptic condition, i.e., there exists a /? > such that, for each A,A^ S{d) 
with A> A, 

G{A)-G{A) >/3tiiA~A]. 

For each ^ = (C^)j^=i € M2(0,T;M'^) and 77 = {v'^)tj^i € (0, T; S(d)), 
we use the following notations 

[ {^udBt):=J2 [ ^dBl- [ {ijud{B)t):=Y, [ 4=d{B\B^)^. 
Jo Jo Jo ^j^^ Jo 

We first consider the representation of (p{Bt — Bt^) for < ii < T < 00. 
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Lemma 2.1 Let = (p{Bt ^ Bt^), ip e Cb.L,;p(M''). Then we have the following 
representation: 

e-tK]+ f{PudB^)+ r{f^t,d{B)t)~ r2G{vt)dt. 
Jtl Jtl Jtl 

Proof. We know that u{t, x) — E[(p{x + Bt — Bf )] is the solution of the following 
PDE: 

dtu + G{D^u) = (i, a;) G [0, T] x ^(r, a;) = ip{x). 
For each e > 0, by the interior regularity of u (see Appendix C), we have 

||u|lci+o/2,2+c.([o,T-c]xB<i) < °° foi" some a £ (0,1). 

Applying G-Ito's formula to u{t,Bt — Bt-^) on [ti,T — e], since Du{t,x) is uni- 
formly bounded, letting e — >■ 0, we have 

C = IE[e]+ / dtuit,Bt- Bt,)dt+ [ {Du{t,Bt- Bt,),dBt) 
ti 



1 '-^ 



{D'u{t,Bt~Bti),d{B)t) 

T 1 i-T 



m+J^ {Duit,Bt^Bt,),dBt) + ^J^ {DM-t,Bt-BtJ,d{B)t) 

- [ G{D^u{t,Bt-BtJ)dt. 
Jtl 



□ 

We now give the representation theorem of ^ v(^ti j ~ ^ti i ' ' ' > ~ 

Theorem 2.2 Let ^ = ^{Bti, Bt, - Bti, ■ ■ ■ ,Bt^-Bt^_,), (ys e Cf,.L»p(R''^'^), 
0<ii<i2<'''<^A'=r<oo. T/ien we have the following representation: 

^ = m+ r {l3udBt)+ r {r^ud{B)t)- 2G{rjt)dt. 
Jo Ja Jo 

Proof. We only need to prove the case ^ = fiBti , Bt — Bti). We set, for each 
(x,y) eM2d^ 

uit.x.y) =t[ip{x,y + Bt - Bt)]] ipi{x) = E[(^(a;, St - J]. 
For each x e K'^, we denote ^ = f{x, Bt — Bti). By Lemma [2. 1[ we have 

^"=^i(x)+^ {Dyu{t,x,Bt~ Bti),dBt) + ^ {Dlu{t,x,Bt~ Bti),d{B)t) 



f G{Dlu{t,x,Bt-Bti))dt. 
Jtl 
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By the definitions of tlic integrations of dt, dBt and d{B)t, we can replace x by 
Bt^ and get 



C = ¥'i(StJ+ / {Dyu{t,Bt,,Bt- Bt,),dBt) 
Jfi 

+ {Dlu{t,Bt,,Bt-Bt,),d{B)t)- G{Dlu{t,Bt,,Bt-Bt,))dt. 



Applying Lemma [2. II to ipi{Btj), we complete the proof. □ 
We then immediately have the following G-martingale representation theo- 
rem. 

Theorem 2.3 Let {Mt)te[a,T] be a G-martingale with Mt = (p{Btj^, Bt^—Bt^, ■■ ■ 
Bt„^,), LP e Gb.L^p{M.'^'^^), 0<ti<t2<---<tN = T<oo. Then 

Mt=t[MT]+ I {l3s,dBs)+ I {T^s,d{B)s)~ f 2G(?7,)ds, t<T. 



Proof. For Mt, by Theorem 12.21 we have 

MT=nMT]+ [ {l3s,dBs)+ I {Tis,d{B)s)- I 2G{i^s)ds. 
Jo Jo Jo 

Taking the conditional G-expectation on both sides of the above equality and 
by Proposition [LJl we obtain the result. □ 



§3 G— convexity and Jensen's Inequality for G— 
expectations 

A very interesting question is whether the well-known Jensen's inequality still 
holds for G-expectations. 

First, we give a new notion of convexity. 

Definition 3.1 A continuous function /i : M — R «s called G-convex if for 
each bounded ^ G the following Jensen's inequality holds: 

E[hio] > him). 

In this section, we mainly consider G^-functions. 

Proposition 3.2 Let h G G'^(K). Then the following statements are equivalent: 

(i) The function h is G -convex. 

(ii) For each bounded ^ G LQ{fl), the following Jensen's inequality holds: 

nhiom > him\nt]) fort>o. 
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(iii) For each G C^{R'^), the following Jensen's inequality holds: 

t[h{<p{Bt))] > h{t[<p{Bt)]) for t > 0. 

(iv) The following condition holds for each {y, z, A) e R x R'' x S{d): 

G{h'{y)A + h"{y)zz^) - h'{y)G{A) > 0. (3.2) 
To prove the above proposition, we need the following lemmas. 
Lemma 3.3 Let $ : M"^ — )• S{d) be continuous with polynomial growth. Then 

lini]E[ r\HBs), d{B),)]6~' = 2E[G($(i3t))]. (3.3) 

<5iO 

Proof. If $ is a Lipschitz function, it is easy to prove that 

ft+S 

E[| ($(B,) - $(Bt), d{B),)\] < C,S^/\ 
where Ci is a constant independent of S. Thus 

rt+S 

nmE[ mB,),d{B),)]d-' = UmE[($(B,), {B)t+5 - {B)s)]5-' 

= 2E[G($(Bt))]. 

Otherwise, we can choose a sequence of Lipschitz functions : — > S{d) 
such that 

\^n{x)-^x)\< ^{1 + 1x1% 
where C2 and k are positive constants independent of A''. It is easy to show that 

ft+S p 

n\j^ mB,)-<^N{Bs),d{B),)\]<-s 

and 

E[|G(<i>(Bi))-G($Ar(St))l]<^, 
where G is a universal constant. Thus 

rt+S 

\E[ ($(B«), d{B),)]S-' - 2E[GmBtm 

{<!>N{B,),d{B),)]6-' - 2E[G(cI>jv(St))]| + ^. 



Then we have 



limsup|E[/ {^{B,),d{B),)]5-^ -2-k[G{^{B,))]\ < — 
54,0 Jt 



54,0 Jt 

Since A'^ can be arbitrarily large, we complete the proof. □ 
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Lemma 3.4 Let 5" fee a C^-function on S.'^ such that D^'i> satisfy polynomial 
growth condition. Then we have 

\im{E['f (Bs)] - *(0))(5-i = G(i:'2*(0)). (3.4) 

6-10 

Proof. Applying G-Ito's formula to ^'(-65), we get 

^{Bs) = ^{0)+ [ {D^iB,),dB,) + ^ [ {DH{B,),d{B)s). 
Jo ^ Jo 

Thus we have 

E[*(i?5)] - vI/(0) = \tG{ [\DHiBs),d{B)s)]. 
^ Jo 

By Lemma [3. 3[ we obtain the result. □ 

Lemma 3.5 Let h £ C^(M) and satisfy (|3.2p . For each tp e Cb,Lip{S.'^), let 
u{t, x) be the solution of the G-heat equation: 

dtu - G{D'^u) = [t, x) e [0, 00) X m(0, x) = Lp{x). (3.5) 

Then u{t,x) := h(u{t,x)) is a viscosity suhsolution of G-heat equation p.Sp with 
initial condition u(0,a;) — h{ip{x)). 

Proof. For each e > 0, we denote by u'^ the solution of the following PDE: 

atu^ - Ge(L'^u^) {t,x) e[Q,oo)xK^, {Q,x) ^ Lp{x), 

where Ge(yl) := G{A) + £tr[v4]. Since G^ satisfies the uniformly elliptic condi- 
tion, by Appendix C, we have E G^'^((0,oo) x M.'^). By simple calculation, 
we have 

dth{u') = h'{u')dtu' = h'{u')G,{D^u') 

and 

dth{u') - G,{D%{u')) = f,{t,x), h{u'{0,x)) = h{^{x)), 

where 

f,{t, x) = h'{u'')G{D^u^) - G{D^h{u^)) - eh"{u'')\Du''\^. 

Since h is G-convex, it follows that fe < —eh"{u'^)\Du'^\'^. We can also deduce 
that \Du'^\ is uniformly bounded by the Lipschitz constant of (p. It is easy to 
show that uniformly converges to u as e — >■ 0. Thus h{u^) uniformly converges 
to h{u) and h" [u'^) is uniformly bounded. Then we get 

dth{u') - Ge{D^h{u')) < Ce, h{u'{0, x)) = h{(p{x)), 

where G is a constant independent of e. By Appendix C, we conclude that h{u) 
is a viscosity subsolution. □ 
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Proof of Proposition [3T2l Obviously (ii)=>(i)=>(iii). We now prove 
(iii)=^(ii). For ^ e L^^in) of the form 

^ = (p{Bt^,Bt^ - Bt^, - ■ ■ , Bt^ - Bt^_^), 

where ^ G Cl{W^'"^), < < • • • < t„ < cx), by the definitions of IE[-] and 
]E[-|r2t], we have 

t[h{om > Hmirtt]), t > 0. 

We then can extend this Jensen's inequality, under the norm || • || E[| • |], to 
each bounded ^ € Lq{Q). 

(iii) =4>(iv): for each (p e C'^{M.'^), we have E[h{if{Bt))] > h{E[if{Bt)]) for each 
t > 0. By Lemma [331 we know that 

\hn{EMBs)] ~ m)S-' = G(i?V(0)) 

04,0 

and 

lim(E[;j(^(i3^))] - h{^{0)))S~' - G{D^h{^m). 

o-lO 

Thus we get 

GiD'hi^m) > h'{^{0))G{D\iO)). 

For each (y,z, A) e M x M'' x S(d), we can choose a. ip Q C^(M'') such that 
((^(0), £1(^(0), £12(^(0)) = Thus we obtain (iv). 

(iv) =^(iii): for each (/? e Cj(R''), u{t,x) = ]E[(/3(a; + i?()] (respectively, u(t,a;) = 
E,[h{ip{x + Bt))]) solves the G-heat equation (|3.5p . By Lemma [3.5( h{u) is a 
viscosity subsolution of G-heat equation (|3.5p . It follows from the maximum 
principle that h{u{t,x)) < u{t,x). In particular, (iii) holds. □ 

Remark 3.6 In fact, (i)-4=>(ii)-4=>(iii) still hold without the assumption h G 
C2(R). 

Proposition 3.7 Let h be a G-convex function and X G LQ{fl) be bounded. 
Then Yt — /i(]E[X|f7t]), t > 0, is a G-submartingale. 

Proof. For each s < t, 

t[Yt\ns] = E[hiE[x\nt])\ns] > h{E[x\ns]) = 

□ 

Exercise 3.8 Suppose that G satisfies the uniformly elliptic condition and h G 
C2(M). Show that h is G-convex if and only if h is convex. 
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Notes and Comments 

This chapter is mainly from Peng (2007) [102]. 

Peng (1997) introduced a filtration consistent (or time consistent, or 
dynamic) nonlinear expectation, called (^-expectation, via BSDE, and then in 

(1999) ilMj for some basic properties of the g-martingale such as nonlinear Doob- 
Meyer decomposition theorem, see also Briand, Coquet, Hu, Memin and P eng 

(2000) Q, Chen, Kulperger and Jiang (2003) [13, Chen and Peng (1998) 
and (2000) jpL Coquet, Hu, Memin and Peng ( 2001) [H], and (2002) [23, 
Peng (1999) [il, (2004) Peng and Xu (2003) [io3, Rosazza (2006) [lT2| . 
Our conjecture is that all properties obtained for (/-martingales must has its 
correspondence for G-martingale. But this conjecture is still far from being 
complete. Here we present some properties of G-martingales. 

The problem G- marti ngale representation theorem has been raised as a prob- 
lem in Peng (2007) |l02|. In Section 2, we only give a result with very regular 
random variables. Some very interesting developme nts t o this important prob- 
lem can be found in Soner, Touzi and Zhang (2009) [III and Song (2009) [1^. 

Under the framework of ^-expectation, Chen, Kulperger and Jiang (2003) 
[2^, Hu (2005) (6^, Jiang and Chen (2004) investigate the Jensen's in- 
equality for ^-expectation. Recently, Jia and Peng (2007) [ill introduced the 
notion of (/-convex function and obtained many interesting properties. Certainly 
a G-convex function concerns fully nonlinear situations. 
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Chapter V 

Stochastic Differential 
Equations 

In this chapter, we consider the stochastic differential equations and backward 
stochastic differential equations driven by G-Brownian motion. The conditions 
and proofs of existence and uniqueness of a stochastic differential equation is 
similar to the classical situation. However the corresponding problems for back- 
ward stochastic differential equations are not that easy, many are still open. We 
only give partial results to this direction. 

§1 Stochastic Differential Equations 

In this chapter, we denote by Mg(0, T; R"), p > 1, the completion of Mg'°(0, T; R") 

under the norm {J^^ E[\T]t\P]dtf^P . It is not hard to prove _that Mg(0,T;R") C 
(0, T; R"). Wc consider all the problems in the space Mg(0, T; R"), and the 
sublinear expectation space (rij'HjE) is fixed. 

We consider the following SDE driven by a rf-dimensional G-Brownian mo- 
tion: 

Xt = Xo+ f b{s, X,)ds+ f h,,{s, Xs)d {B\ B')+ f a,{s, X,)dBi, t G [0, T], 
Jo Jo Jo 

(1-1) 

where the initial condition Xo G R" is a given constant, and 6, hij, aj are given 
functions satisfying a;) , hij{-,x), (Tj{-,x) € Mq(0,T;R") for each a; € R" and 
the Lipschitz condition, i.e., \(j){t,x) — (p{t,x')\ < K\x — x'\, for each t g [0,r], 
X, x' G R", (/) = 6, hij and aj, respectively. Here the horizon [0,T] can be 
arbitrarily large. The solution is a process X G Mg.(0,T;R") satisfying the 
SDE (fTTD . 

Wc first introduce the following mapping on a fixed interval [0, T]: 
A. : M^(0,T;R") M^(0,T;R") 
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by setting At, t e [0,T], with 

AtiY)=Xo+ f b{s,Y,)ds+ f h,j{s,Ys)d{B\B=)^+ [ a,{s,Ys)dBl. 
Jo Jq Jo 

We immediately have the following lemma. 
Lemma 1.1 For each Y,Y' e A/q(0,T;M"), we have the following estimate: 

n\MY)-AtiY')\']<C fE[\Ys-Y:\^]ds, t e [0,T], (1.2) 

Jo 

where the constant C depends only on the Lipschitz constant K . 

We now prove that SDE (jl.ip has a unique solution. By multiplying e^^*^* 
on both sides of ()1.2p and integrating them on [0,r], it follows that 

r t[\At{Y)-At{Y')\^]e-^'''dt<C Te-^^* f tcWs -Y^\^]dsdt 
Jo Jo Jo 

= C / e-^^*diE[|y, - y;flds 



\ I (e-2^^-e-2c^^)]E[|i;-y;p]rf.. 

^ Jo 



2 

We then have 

^ t[\At{Y)~At{Y')\^]e-^'''dt<\ rE[\Yt-Y;\^]e'^^'dt. (1.3) 
We observe that the following two norms are equivalent on M^(0,T;]R"), i.e., 

( rn\Yt\W/' ~ ( r]E[|r,|2]e-2C*di)i/2. 
Jo Jo 

From (|1.3p we can obtain that A{Y) is a contraction mapping. Consequently, 
we have the following theorem. 

Theorem 1.2 There exists a unique solution X G Mq(0,T;M") of the stochas- 
tic differential equation (jl.ip . 

We now consider the following linear SDE. For simplicity, we assume that 
d = 1 and n = I. 

Xt - Xo+ f ibsXs+k)ds+ f {h,Xs+k)d{B)s+ f {asXs+a,)dBs, t G [0,T], 
Jo Jo Jo 

(1.4) 

where Xq G M is given, fe., /i., a, are given bounded processes in Mq{0, T\ R) and 
h,,h,, a , are given processes in Mq(0, T; M). By Theorem II. 2 [ we know that the 
linear SDE (II. 4p has a unique solution. 
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Remark 1.3 The solution of the linear SDE (|1.4p is 

Xt^T-\XQ + f kT,ds+ f {hs-asds)Tsd{B)s+ f a,T,dB,), te[0,T], 
Jo Jo Jo 

where Tt = exp(- /g bsds - ^^{hs - ^a^)d{B)s - (JsdBs). 

In particular, ij h,,h,,a, are constants and h,,h,,(j, are zero, then X is a 
geometric G-Brownian motion. 

Definition 1.4 We call X is a geometric G-Brownian motion if 

Xt^exp{at + /3{B)t+-fBt), (1.5) 
where a, /3, 7 are constants. 

Exercise 1.5 Prove that Mg(0,T;R") C Mg(0,T;R"). 
Exercise 1.6 Complete the proof of Lemma \l.l\ 



§2 Backward Stochastic Differential Equations 

We consider the following type of BSDE: 

Yt^n^ + J^ f{s,Y,)ds + J^ h,j{s,Y,)d{B\B^)jnt], te[0,T], (2.6) 

where ^ £ ^^(riT ;]&"■) is given, and /, /ly are given functions satisfying f{-,y), 
hij{-,y) S M^(0,T;R") for each y G R" and the Lipschitz condition, i.e., 
|(/.(t,y) - 0(<,y')l < K\y - y'l for each t G [0,r], y, y' £ R", (f> ^ f and 
hij, respectively. The solution is a process Y £ Mq{0,T;M.") satisfying the 
above BSDE. 

We first introduce the following mapping on a fixed interval [0,T]: 
A. : Af^(0,r;R") ^ A7^(0,r;R") 
by setting At, t £ [0,T], with 



At{Y)=E[i + J^ fis,Ys)ds + J^ h,,is,Y,)d{B\B^)JQt]. 

We immediately have 
Lemma 2.1 For each Y,Y' £ Mq{Q,T;W^), we have the following estimate: 

t[\At{Y)~At{Y')\]<cj\[\Y,-Y^\]ds, t£[0,T], (2.7) 
where the constant C depends only on the Lipschitz constant K . 
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We now prove that BSDE (|2.6p has a unique solution. By muhiplying e^*^* 
on both sides of (|2.7p and integrating them on [0,T], it follows that 



T rT 



[ E[\At{Y)-At{Y')\]e'^'dt<cf f 
Jo Jo Jt 



\Ys-Y^\]e^^*dsdt 



T 



= C / E[|r, - / e^^'dtds 
Jo Jo 

= ^£n\Ys-YMe''''-l)ds 
<^£Em^Y:\]e'^^ds. (2.8) 
We observe that the following two norms are equivalent on M^(0,T;R"), i.e., 

E[\Yt\]dt^ r E[\Yt\]e^'''dt. 



From (|2.8|) . we can obtain that A{Y) is a contraction mapping. Consequently, 
we have the following theorem. 

Theorem 2.2 There exists a unique solution (i^t)te[o.T] S Mg.(0,T;R") of the 
backward stochastic differential equation ()2.6p . 



Let Y^ , V — 1,2, be the solutions of the following BSDE: 

Yi' = tr+ r {f{s,Y:) + v:)ds+ r{h,,{s,Y:)+rr)d{B\B^)jn,]. 



Then the following estimate holds. 
Proposition 2.3 We have 

where the constant C depends only on the Lipschitz constant K . 
Proof. Similar to Lemma [2Tl we have 

- Y^\] < cij^nY} - Y^\]ds + E[\e - e\] 

By the Gronwall inequality (see Exercise 12. 5p . we conclude the result. □ 
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Remark 2.4 In particular, if = 0, i/j^ = — /(s, 0), -0s'''^ — —hij{s,0), ip] — 
0, ipl'''^ = 0, we obtain the estimate of the solution of the BSDE. Let Y be the 
solution of the BSDE (|2.6p . Then 

E[|rt|]<Ce^(^-*)(E[|e|]+^ E[|/(s,0)| + |/iy(s,0)|]ds), (2.10) 

where the constant C depends only on the Lipschitz constant K . 

Exercise 2.5 (The Gronwall inequality) Let u{t) be a nonnegative function 
such that ^ 

uit) <C + A u{s)ds forO<t<T, 
Jo 

where C and A are constants. Prove that u{t) < Ce^* for < t < T . 

Exercise 2.6 For each ^ S ^^(riTiR")- Show that the process (^[£,\^t])te[o,T] 
belongs to M^{0,T;M."). 

Exercise 2.7 Complete the proof of Lemma \2.1\ 



§3 Nonlinear Feynman-Kac Formula 

Consider the following SDE: 

r dXl'^=h{Xl-'^)ds + h,,{Xl'^)d{B\B^)^+a,{Xl'^)dBl [t,T], 

where ^ G LQ(51t;M") is given and 6, hij, aj : R" — > M" are given Lipschitz 
functions, i.e., \(t){x) — (t){x')\ < K\x — x'\, for each x, x' G M", (f) — h, hij and 

We then consider associated BSDE: 

= E[$(X^«) + f{X*/X'^}dr + J"^ g^,{X'r^X'^)d{B\B')^\ns], 

(3.12) 

where $ : M" — > M is a given Lipschitz function and /, gij : R" x R — > R are 
given Lipschitz functions, i.e., \(/){x,y) — (/){x',y')\ < K{\x — x'\ + \y — y'\), for 
each X, x' G R", y, y' £ R, (j) — f and gij. 
We have the following estimates: 

Proposition 3.1 For each ^, G iQ(f24;M"), we have, for each s G [t,T], 

E[\x'/ - x'/\^\nt] < c|e - (3.13) 

and 

nx'/\'m<c{i + \cf), (3.14) 

where the constant C depends only on the Lipschitz constant K . 



86 



Chap.V Stochastic Differential Equations 



Proof. It is easy to obtain 

E[\x'/ - xl'i'\'\nt] < Ci(|e - + nx'/^ - X'/'\'\nt]dr). 

By the Gronwall inequality, we obtain 

E[\X'/ - X'/\^\nt] < Cie^i^ie - • 
Similarly, we can get p.l4p . □ 
Corollary 3.2 For each £, e L%{nt;m."), we have 

n\X'4,-^\'\nt]<C{l + \e)S forSe[0,T~t], (3.15) 
where the constant C depends only on the Lipschitz constant K . 
Proof. It is easy to obtain 

i-t+S 

E[|x;;«, - ^na] < Ci (l + ]E[|X*^«na])ds. 

By Proposition 13. 1[ we obtain the result. □ 
Proposition 3.3 For each ^, ^' e iQ(f24;M"), we have 

|r/'«-r/'«'|<qc-c'| (3.16) 

and 

|y/'«|<C(l + |C|), (3.17) 
where the constant C depends only on the Lipschitz constant K . 

Proof. For each s E [0,r], it is easy to check that 

_ < c^E[\x'/ - x'/\ + j\\xl'^ - xl^^'\ + \y^'^ - y;^«'|)dr|i],]. 

Since 

t[\xl'^ - x'/'\\n,] < (E[|x*^« - x*^«'ni7,])V2, 

we have 

By the Gronwall inequality, we obtain p.l6p . Similarly we can get (I3.17p . □ 
We are more interested in the case when ^ = x G R". Define 

u{t, x) := r/'"^, {t, x) e [0, T] X M". (3.18) 

By the above proposition, we immediately have the following estimates: 

\u{t,x)-u{t,x')\<C\x~x'l (3.19) 



\u{t,x)\ < C(l + |a;|), 
where the constant C depends only on the Lipschitz constant K. 



(3.20) 
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Remark 3.4 It is important to note that u{t,x) is a deterministic function of 
{t,x), because Xl'^ andY^'^ are independent from ^It- 

Theorem 3.5 For each ^ G L'^{Vlt]W^), we have 

uit,O^Y,''i. (3.21) 

Proposition 3.6 We have, for S £ [0,T ~t], 

u{t,x) = t[u{t+5,Xl^^,)+ j f{Xy\Y::^^)dr+ j B^) J 

(3.22) 

Proof. Since = xT'""*'^' for s S [i + 5, T], we get Y^f^ = i"t+/^^*"'- By 
Theorem 13. 51 we have Y^^^ = u{t + 5,Xl'^g), which imphes the resuh. □ 
For each A £ §(n), p e R", r e K, we set 

F{A,p,r,x) GiB{A,p,r,x)) + {p,b{x)) + f{x,r), 

where B(A,p,r,x) is a d x c? symmetric matrix with 

Bij{A,p,r,x) := {Aa.i{x),aj{x)) + {p,hij{x) + hji{x)) + gij{x,r) + gji{x,r). 

Theorem 3.7 u{t, x) is a viscosity solution of the following PDE: 

dtu + F{D^u,Du,u,x) ^ 0, , . 

u{T,x) = ^{x). 

Proof. Wc first show that u is a continuous function. By (j3.19p wc know that u 
is a Lipschitz function in x. It foUows from ()2.10p and p.l4p that for s E [t, T], 
]E[|X*'^|] < C(l + |a;|). Noting (l3T5ll and jS^S), we get \u{t, x) - u{t + 6, x)\ < 
C{l + \x\){S'^^'^ +S) for 5 e[0,T-t]. Thus u is ^-Holder continuous in t, which 
imphes that it is a continuous function. We can also show, that for each p > 2, 

M^l+s ~ < C(l + \x\P)6P/^, (3.24) 

Now for fixed {t,x) € (0,r) x M", let V G C'^'^{[0,T] x M") be such that ip > u 
and il^{t,x) = u{t,x). By (|3.22p . p.24p and Taylor's expansion, it follows that, 
for ^ e (0,T-i), 

rt+S 

< Em + 5, Xlf,) - i^it, x)+ J f{Xi:-Xndr 

i-t+S 



g,,{Xl.-\Y^^-)d{B\B')^] 

Jt 

< ^E[iBiD^^it,x),D^it,x),^it,x),x),{B)t+s - (B)t)] 

+ (9tV(i, x) + {Di,{t, x), b{x)) + fix, i,{t, x)))6 + C{1 + \x\ + |xp + jx^js/^ 

< {dtip{t, x) + FiD^it, x), D^j{t, x),^j{t, x),x))S + C(l + \x\ + \x\^ + \xf)S^^^, 
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then it is easy to check that 

dtip{t, x) + F{D^'ip{t, x), Dip{t, x),ip{t, x),x)> 0. 

Thus u is a viscosity subsolution of p.23p . Similarly we can prove that u is a 
viscosity supersolution of p.23l) . □ 

Example 3.8 Let B = {B^,B^) he a 2 -dimensional G-Brownian motion with 

G(A) =Gi(aii) + G2(a22), 

where 

In this case, we consider the following 1-dimensional SDE: 

dXl-"" = fiXl^'^ds + vXl'''d{B^)^ + aXl'^dBl, Xl'"" = x, 

where ^, v and a are constants. 

The corresponding function u is defined by 

u{t,x) := E[ip{X*^'')]. 

Then 

u{t,x) ^E[u{t + S,Xl:^s)] 
and u is the viscosity solution of the following PDE: 

dtu + fixdxU + 2Gi{vxdxu) + a'^x'^G2{d1.^u) — 0, u{T,x) — (p{x). 

Exercise 3.9 For each ^ € L^(f^t;M") withp> 2, show that SDE ([3111) has a 

unique solution in M^it, T; R"). Furthermore, show that the following estimates 
hold. 

E[|X*'"-X*'"'|P] < C\x-x'\P, 
E[\Xl'T]<Cil + \x\n, 

E[\xi-,-xn<cii + \x\nsp/'. 

Notes and Comments 

This chapter is mainly from Peng (2007) |l02| . 

There are many excellent books on Ito's stochastic calculus and stochastic 
differential equations founded by Ito's original paper (65j . as well as on martin- 
gale theory. Readers are referred to Chung and Williams (1990) [2^, Dellacherie 
and Meyer (1978 and 1982) fsS^, He, Wangand Yan (1992) "sT^, Ito and McKean 
(1965) Ikeda and Watanabe (1981) [pi, Kallenberg (2002}J72,], Karatzas 
and Shreve (1988) [zl], 0ksendal (1998) |87|, P rotter (1990) [ll37R-evuz and 
Yor (1999) [m| and Yong and Zhou (1999)[l2Jl. 
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Linear backward stochastic differential equation (BSDE) was first introduced 
by Bismut in (1973) fr2i] and (1978) tlSj]. Bensoussan developed this approach 
in (1981) [l3] and (1982) [Il[. The existence and uniqueness theorem of a gen- 
eral nonlinear BSDE, was obtained in 1990 in Pardoux and Peng [s^l- The 
present version of the proof was based on El Karoui, Peng and Quenez (1997) 
|44| . which is also a very good survey on BSDE theory and its applications, spe- 
cially in finance. Comparison theorem of BSDEs was obtained in Peng (1992) 
(90| for the case when 5 is a C^-function and then in 44 1 when g is Lipschitz. 
Nonlinear Feynman-Kac formula for BSDE was introduced by Peng (1991) [s^ 
and (1992) '91]. Here we obtain the corresponding Feynman-Kac formula unde r 
the framework of G-expectation. We also refer to Yong and Zhou (1999) 12j|, 
as well as in Peng (1997) ^ (in 1997, in Chinese) and (2004) ^ for system- 
atic presentations of BSDE theory. For contributions in the developments of 
this theory, readers can be referred to the literatures listing in the Notes and 
Comments in Chap. HI 
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Chapter VI 



Capacity and Quasi-Surely 
Analysis for G-Brownian 
Paths 

In this chapter, wc first present a general framework for an upper expectation 
defined on a metric space B{^)) and the corresponding capacity to introduce 
the quasi-surely analysis. The results are important for us to obtain the pathwise 
analysis for G-Brownian motion. 

§1 Integration Theory associated to an Upper 
Probability 

Let be a complete separable metric space equipped with the distance d, 

the Borel cr-algebra of O and A4 the collection of all probability measures on 

(0,B(f2)). 

• iy°(ri): the space of all S(f2)-measurable real functions; 

• Bb{fl): all bounded functions in L°{Q); 

• C(,(f2): all continuous functions in Bb{Q). 

All along this section, we consider a given subset V C M.. 

1.1 Capacity associated to V 

We denote 

c{A) := sup P{A), A € B{n). 

Pev 

One can easily verify the following theorem. 



91 



92 



Chap. VI Capacity and Quasi-Surcly Analysis for G-Brownian Paths 



Theorem 1.1 The set function c{-) is a Choquet capacity, i.e. (see f22. Is^]), 

1. 0< c{A) < 1, VA c n. 

2. If Ac B, then c{A) < c{B). 

3. If {An)^^i is a sequence in B{ri), then c(U^„) < ^c(A„). 

4- If {An)^^i is an increasing sequence in B{VL): An t ^ = Uv4„, then 
c{\jAn) = lini„^oo ciAn). 

Furthermore, we have 

Theorem 1.2 For each A e B{n), we have 

c{A) = sup{c{K) : K compact K C A}. 
Proof. It is simply because 

c{A) ~ sup sup P{K) = sup sup P{K) = sup c{K). 

PG'P K compact K compact PGV K compact 

ifCA K(ZA KCA 

□ 

Definition 1.3 We use the standard capacity-related vocabulary: a set A is 
polar if c{A) — and a property holds "quasi- surely" (q.s.) if it holds outside 
a polar set. 

Remark 1.4 In other words, A £ B{fl) is polar if and only if P{A) — for 
any P CzV. 

We also have in a trivial way a Borel-Cantelli Lemma. 
Lemma 1.5 Let (j4„)„gN be o, sequence of Borel sets such that 

oo 

c{An) < OO. 

Then lim sup„_^g^ An is polar . 

Proof. Applying the Borel-Cantclli Lemma under each probability P E V. □ 
The following theorem is Prokhorov's theorem. 

Theorem 1.6 V is relatively compact if and only if for each e > 0, there exists 
a compact set K such that c{K'^) < e. 

The following two lemmas can be found in [g^. 

Lemma 1.7 V is relatively compact if and only if for each sequence of closed 
sets Fn i 0, we have c{Fn) 4 0. 
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Proof. We outline the proof for the convenience of readers. 
"=>" part: It follows from Theorem 11.61 that for each fixed e > 0, there exists 
a compact set K such that c{K'^) < e. Note that F„ n | 0, then there exists 
an > such that FnHK = $ for n > N, which implies lim„ c(F„) < e. Since 
e can be arbitrarily small, we obtain c(i^„) I 0. 

"<;=" part: For each e > 0, let (A'l^)°^^ be a sequence of open balls of radius 
1/fc covering Q. Observe that {U^^iA'^Y I 0, then there exists an Uk such that 
c((U"iiAf)^) < £2-'=. Set K = n°l7i4S7^- It is easy to check that K is 
compact and c{K'^) < e. Thus by Theorem 11.61 is relatively compact. □ 

Lemma 1.8 Let V he weakly compact. Then for each sequence of closed sets 
Fn i- F , we have c{Fn) i c(F). 

Proof. We outline the proof for the convenience of readers. For each fixed e > 0, 
by the definition of c(F„), there exists a. Pn & V such that P„(F„) > c(F„) — e. 
Since V is weakly compact, there exist P„j. and P G V such that P„j. converge 
weakly to P. Thus 

P{Fm) > limsupP„JP™) > limsupF„jF„J > lim c(F„) - e. 

k^oo k^oa 

Letting m — )■ oo, we get P{F) > lim„^oo c(P„) — e, which yields c(F„) J, c{F). 
□ 

Following [H^ (see also [3l,[53) the upper expectation of V is defined as follows: 
for each X £ L°(f2) such that Ep[X] exists for each P eV, 

E[X] = E^[X] := sup Ep[X]. 

It is easy to verify 

Theorem 1.9 The upper expectation E[-] of the family V is a suhlinear expec- 
tation on Bi,{fl) as well as on Cb{^), i-e., 

1. for all X,Y in Bb{fl), X >Y =^ E[X] > E[Y]. 

2. for all X,Y in Bbin), E[X + Y]< E[X]+E[Y]. 

3. for all X>0, X e Bb{n), E[\X] = \E[X]. 

4. for allceR, X e Bbin) , E[X + c] = E[X\ + c. 
Moreover, it is also easy to check 

Theorem 1.10 We have 

1. LetE[Xj^ andEE^^iX„] he finite. T/ien EE^^, X„] < ^^^^ E[X„]. 

2. Let X,, t X and E[X„], E[X] he finite. Then E[X„] t E[X]. 

Definition 1.11 The functional E[-] is said to be regular if for each {X„]^^i 
in Cfc(r2) such that Xn i on D,, we have E[Xn] \. 0. 
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Similar to Lemma 11.71 we have: 

Theorem 1.12 E[-] is regular if and only ifV is relatively compact. 

Proof. part: For each sequence of closed subsets F„ J, such that 

Fn, n = 1,2,- ••, are non-empty (otherwise the proof is trivial), there exists 
{gn}^=i C Cb(r2) satisfying 

< 5n < 1, gn ^ ^ on Fn and (;„ = on {oj e 51 : d(uj, Fn) > —}■ 

n 

We set /„ — Af^igt, it is clear that /„ e Cb{^) and 1f„ < fn i 0. E[-] is 
regular implies E[/„] l and thus c{Fn) i 0. It follows from Lemma [1.71 that V 
is relatively compact. 

part: For each {Xn}^^^ C Cf,(il) such that Xn i 0, we have 

poo poo 
E[Xn] - sup Ep[Xn] = sup / P({X„ > t})dt < / c{{X„ > t})dt. 

Pev Per Jo Jo 

For each fixed t > 0, {Xn > i} is a closed subset and {Xn > i} i as n f oo. 
By Lemma [TTTl c({X„ > t}) 10 and thus c({X„ > t})dt | 0. Consequently 

nxn] ; 0. □ 

1.2 Functional spaces 

We set, for p > 0, 

• CP := {X e L°{n) : E[\X\p] = supp^-p Ep[\X\p] < oo}; 

• AfP := {X e L^\n) : E[\X\p] = 0}; 

• TV := {X e : X = 0, c-q.s.}. 

It is seen that CP and A/''' are linear spaces and AfP = J\f, for each p > 0. 

We denote := CP/Af. As usual, we do not take care about the distinction 

between classes and their representatives. 



Lemma 1.13 Let X eEP. Then for each a> 

\X\P] 



c{{\X\>a})< 



aP 



Proof. Just apply Markov inequality under each P E V. □ 
Similar to the classical results, we get the following proposition and the proof 
is omitted which is similar to the classical arguments. 

Proposition 1.14 We have 

1. For each p > 1, L'' is a Banach space under the norm \\X\\,^ := {E[\X\p])p . 
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2. For each p < 1, is a complete metric space under the distance 
d{X,Y) :=E[|X-y|f]. 

We set 

:= {X G L"{n) : 3 a constant M, s.t. \X\ < M, q.s.}; 
L°° C^/Af. 

Proposition 1.15 Under the norm 

\\X\\^:^M{M>0:\X\<M, q.s.}, 

L°° is a Banach space. 

Proof. From {\X\ > \\X\\^} = i {\X\ > \\X\\^ + 1} we know that \X\ < 
||X||^, q.s., then it is easy to check that \\-\\^ is a norm. The proof of the 
completeness of L°° is similar to the classical result. □ 
With respect to the distance defined on L^, p > 0, we denote by 

• the completion of i?b(ri). 

• LP the completion of Cb(ri). 
By Proposition II .141 we have 

LP c LP c LP, p>0. 
The following Proposition is obvious and the proof is left to the reader. 

Proposition 1.16 We have 

1. Lef p, g > 1, i + i = 1. Then X e LP andY e L"? implies 

XY e L^ and E[|Xr|] < (E[|X|p])p (E[|r|«])^ ; 

Moreover X eU^ and F e Lf implies XY e L;i; 

2. LPi C LP^ Lf C Lf , LPi C LP^ < P2 < Pi < oo; 

3. \\X\\^^ \\X\\^, for each X ^l.^ . 

Proposition 1.17 Letp e (0, oo] and (X„) he a sequence in LP which converges 
to X ini^P . Then there exists a subsequence (Xn^) which converges to X quasi- 
surely in the sense that it converges to X outside a polar set. 
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Proof. Let us assume p € (0, oo), the case p — co is obvious since the conver- 
gence in L°° implies the convergence in for all p. 
One can extract a subsequence {Xn^ ) such that 

We set for all k 

Ak^{\X~Xr,,\>l/k}, 

then as a consequence of the Markov property (Lemma I1.13P and the Borel- 
Cantelli Lemma fl.51 c(limfe_).oo^fe) = 0. As it is clear that on (limi;_j.oo^fc)^, 
(Xn^) converges to X, the proposition is proved. □ 
We now give a description of . 

Proposition 1.18 For each p > 0, 

hl = {XehP: hm E[\X\n^^x\>n}] - 0}. 

Proof. We denote Jp = {X e hP : lim„^oc IE[|^|''l{|x|>n}] = 0}. For each 
X e Jp let X„ = (X A n) V (-n) G Bb{^l). We have 

E[\X - X„|f] < E[\X\n[\x\>n}] ^ 0, as n ^ (X3. 

Thus X GhP. 

On the other hand, for each X e L|^, we can find a sequence {^nj^i in Bi,{n) 
such that E[|X-F„|P] 0. Let y„ = sup^^^ l^nC^^)! and X„ = (XAy„)V(-?/„). 
Since \X - X„\ < \X-Y„\, we have E[\X - Xn\P] 0. This clearly implies that 
for any sequence (a„) tending to oo, lim„^oo E[|^ — (X A a„) V (— = 0. 
Now we have, for all n E N, 

E[|X|f l{|x|>„}] - E[{\X\ -n + n)n^^x\>n}] 

< (1 V 2P-') {E[{\X\ - nrii\x\>n}] + nPcilXl > n)) . 

The first term of the right hand side tends to since 

E[{\X\ ~ n)f l{|x|>„}] = E[\X -{XAn)V {~n)\P] 0. 

For the second term, since 

^lilXl>u} < i\X\ - ^rM\X\>n} < {\X\ f ri{|X|>t}, 

we have 

-c{\X\ > n) = _E[l{|^|>„j] < E[{\X\ - -ri{\x\>^}] ^ 0. 

Consequently X <E Jp. □ 

Proposition 1.19 Let X eE\. Then for each e > 0, there exists a 5 > Q, such 
that for all A G B{n) with c{A) < 6, we have E[\X\1a] < e. 
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Proof. For each e > 0, by Proposition 11.181 there exists an > such that 
E[|X|l||x|>^|] < §. Take 5^ ^. Then for a subset A G B{n) with c{A) < 6, 
we have 

E[|X|1^] < E[\X\lAl{\x\>N}]+n\X\lAl{\X\<N}] 

<E[\X\l[ix\>N}] + Nc{A) < e. 

□ 

It is important to note that not every element in satisfies the condition 
hm„_j>oo lE[|-''^|''l{|x|>n}] = 0- We give the following two counterexamples to 
show that and are diflferent spaces even under the case that V is weakly 
compact. 

Example 1.20 Let n = V = {Pn ■ n e N} where Pi({l}) = 1 and 
P„({1}) = 1 — i, P„({n}) — for n = 2,3, • • • . V is weakly compact. We 
consider a function X on N defined by X{n) = n, n Cz N. We have E[|^|] = 2 
but E[\X\l[\x\>n}] = 1 7^ 0. In this case, X eE^ but X El. 

Example 1.21 Let 17 = N, P = {P„ : n e N} where Pi({l}) = 1 and 
P„({1}) = 1- Pni{kn}) = ^, k = l,2,...,n,forn = 2,3,---. V is weakly 
compact. We consider a function X on N defined by X{n) — n, n £ N. We 
have E[\X\] - f and nE[l[\x\>n}] - ^ ^ 0, but E[\X\l{\x\>n}] = ^ + ^^0. 
In this case, X is in E^ , continuous and nE[l^^x\>n}] ^ 0, but it is not in Lj. 

1.3 Properties of elements in 

Definition 1.22 A mapping X on D, with values in a topological space is said 
to be quasi- continuous (q.c.) if 

Ve > 0, there exists an open set O withc[0) < e such that X\o<: is continuous. 

Definition 1.23 We say that X : O — > M has a quasi- continuous version if 
there exists a quasi- continuous function y : i7 — > R with X ~Y q.s.. 

Proposition 1.24 Let p > 0. Then each element inE^ has a quasi- continuous 
version. 

Proof. Let (Xn) be a Cauchy sequence in Ch(ri) for the distance on L^*. Let 
us choose a subsequence {Xnk)k>i such that 

E[\X„,^,-X^J]<2~'\ Vfc>l, 

and set for all k, 

oo 
i—k 

Thanks to the subadditivity property and the Markov inequality, we have 

oo oo 
c{Ak) < J2 - ^n. I > 2~^/P) < = 2-'^+!. 

i=k i=k 
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As a consequence, linife_s.oo c{Ak) = 0, so the Borel set A = Hfe";! ^fe is polar. 
As each Xn^ is continuous, for all fc > 1, Ak is an open set. Moreover, for all 
k, (Xn-) converges uniformly on A'j, so that the limit is continuous on each A^. 
This yields the result. □ 
The following theorem gives a concrete characterization of the space LJ?. 

Theorem 1.25 For each p > 0, 

L*! ~ {X CzU' : X has a quasi- continuous version, lim E[|X|^l|ix|>n|] = 0}- 
Proof. We denote 

Jp — {X eLP : X has a quasi-continuous version, lim E[|X|p1||J(: ~ 0}. 

Let X E IJ^, we know by Proposition ll.24l that X has a quasi-continuous version. 
Since X e h^, we have by Proposition 11.181 that lim„^oo Eil-'^l^ldx^n}] = 0. 
Thus X e Jp. 

On the other hand, let X e Jp be quasi-continuous. Define Yn — {X An)\/ (— n) 
for all n G N. As E[\X\Pl{\x\>n}] ^ 0, we have E[\X - r„|P] ^ 0. 
Moreover, for all n € N, as F„ is quasi-continuous , there exists a closed set F„ 
such that c(F^) < -^^^ and F„ is continuous on F„. It follows from Tietze's 
extension theorem that there exists Zn € Cb{^) such that 

\Zn\ < n and Z„ = F„ on i^„. 

We then have 

E[|F„-Z„n<(2nfc(i^„=)<i|g^. 

So E[|X - Z„|f] < (1 V 2P-i)(E[|X - y„|P] +E[|y„ - Z„|P]) ^ 0, and X e L^. 
□ 

We give the following example to show that L*? is different from even under 
the case that V is weakly compact. 

Example 1.26 Let $7 = [0,1], V = {5^ ■ x € [0,1]} is weakly compact. It is 
seen that L^! — Cb(r2) which is different from L^. 

We denote := {X E L°° : X has a quasi-continuous version}, we have 
Proposition 1.27 is a closed linear subspace o/L°°. 

Proof. For each Cauchy sequence {A'„}^-^ of under we can find 

a subsequence {Xni}°Zi such that ||X„;^j — || < 2^'. We may further 
assume that each X„ is quasi-continuous. Then it is easy to prove that for each 
e > 0, there exists an open set G such that c(G) < e and |^Tii+i — \ ^ 2~* 
for all i > 1 on G^, which implies that the limit belongs to L^. □ 
As an application of Theorem 1 1.25[ we can easily get the following results. 
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Proposition 1.28 Assume that X : il — )• M has a quasi- continuous version 
and that there exists a function f : — > satisfying lim4_^oo = oo and 
E[f{\X\)] < oo. Then X € LP. 

Proof. For each e > 0, there exists an > such that ^f- > i, for ah t > N. 
Thus 

n\X\n{\x\>N}] < eE[f{\X\)l{\x\>N}] < eE[f{\X\)]. 
Hence hm^v^oo E[|X|Pl{|x|>Ar}] = 0. From Theorem [TH] we infer X e hP. □ 

Lemma 1.29 Let {Pn}^^i C V converge weakly to P E V . Then for each 
X e Lj, we have Ep^ [X] -> Ep[X]. 

Proof. We may assume that X is quasi-continuous, otherwise we can consider 
its quasi-continuous version which does not change the value Eq for each Q Cz V. 
For each e > 0, there exists an TV > such that E[|X|l{|x|>7v}] < §• Set 
Xn — {X AN)\/ {—N). We can find an open subset G such that c(G) < ^ and 
Xn is continuous on C. By Tietze's extension theorem, there exists Y e Cb(ri) 
such that |r| < and y = Xa, on Obviously, for each Q G P, 

\Eq[X] ~ Eq[Y]\ < Eq[\X -Xn\]+ Eq[\Xn - Y\] 

< £ + 2N— = e. 
- 2 4iV 

It then follows that 

hm sup Ep,^ [X] < hm Ep,^ [Y] + e = Ep [Y] + e < Ep [X] + 2e, 

and similarly liminf„^oo Ep^ [X] > Ep[X]--2e. Since e can be arbitrarily small, 
we then have Ep^ [X] Ep[X]. □ 

Remark 1.30 For continuous X, the above lemma is Lemma 3.8.7 in f/l /. 

Now we give an extension of Theorem 11.121 

Theorem 1.31 Let V he weakly compact and let {Xn}"^^-^ C Lj be such that 
Xn i X, q.s.. Then E[X„] ; ¥.[X]. 

Remark 1.32 It is important to note that X does not necessarily belong to Lj. 

Proof. For the case E[X] > — oo, if there exists a. S > Q such that E[X„] > 
E[X] + S, n = 1,2,- ••, we then can find a Pn £ V such that Ep^[Xn] > 
E,[X] + S — i, n ~ 1,2,---. Since V is weakly compact, we then can find a 
subsequence {Pni}°^i that converges weakly to some P E V. From which it 
follows that 

Ep[Xn,] = lim Ep^,[Xn^] > limsup£;p„. [X„J 

> limsup{E[X] +5 - —} = E[X] + S, i = 1, 2, • • • . 
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Thus Ep[X] > E[X] + S. This contradicts the definition of E[-]. The proof for 
the case E[X] = — oo is analogous. □ 
We immediately have the following corollary. 

Corollary 1.33 Let V be weakly compact and let {^n}^i be o, sequence in 
decreasingly converging to q.s.. Then E[X„] \, 0. 



1.4 Kolmogorov's criterion 

Definition 1.34 Let L be a set of indices, {Xt)t£i and {Yt)tei be two processes 
indexed by L . We say that Y is a quasi-modification of X if for all t £ I , 
Xt=Yt q.s.. 

Remark 1.35 In the above definition, quasi-modiflcation is also called modifi- 
cation in some papers. 

We now give a Kolmogorov criterion for a process indexed by with d G N. 

Theorem 1.36 Let p > Q and {Xt)ti=[Q^i]d be a process such that for all t G 
[0, 1]'', Xt belongs to II . A ssume that there exist positive constants c and £ 
such that 

n\xt-XsY']<c\t-s\''+'. 

Then X admits a modification X such that 



E 



sup 



\Xt - X, 



s^t \t- 



< oo, 



for every a G [0,e/p). As a consequence, paths of X are quasi-surely Hoder 
continuous of order a for every a < e/p in the sense that there exists a Borel 
set N of capacity such that for all w G N'^ , the map t — > X{w) is Hoder 
continuous of order a for every a < e/p. Moreover, if Xt G for each t, then 
we also have Xt G L^. 



Proof. Let D be the set of dyadic points in [0, 1] 



d. 



I "> Or? ' 



Let a G [0,£/p). We set 



,^<iG{0,l,••• ,2"} 



M ■ 



\Xf 
sup — 

s,t<£D,s^t F " 



X, 



Thanks to the classical Kolmogorov's criterion (see Revuz-Yor |lllj ). we know 
that for any P E V, Ep[MP] is finite and uniformly bounded with respect to P 
so that 

E[MP] = sup Ep[MP] < oo. 
Per 
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As a consequence, the map t ~> Xt is uniformly continuous on D quasi-surely 
and so we can define 

yt e [0, l]'*, Xt = lim X,. 

s^t.seD 

It is now clear that X satisfies the enounced properties. □ 

§2 G- expect at ion as an Upper Expectation 

In the following sections of this Chapter, let ft — Cg (M+) denote the space of 
all E**— valued continuous functions (wt)(gR+, with ojq = 0, equipped with the 
distance 

oo 

p{u;\co') := ^2"M(max \u;l ~ A 1], 

and let fl = (R'')['^'°°) denote the space of all M''— valued functions (wt)tgR+. Let 
B{il) denote the cr-algebra generated by all open sets and let B{fl) denote the 
(T-algebra generated by all finite dimensional cylinder sets. The corresponding 
canonical process is Bt{uj) = uit (respectively, Bt{Cj) = Wt), t G [0, oo) for w S 57 
(respectively, cZi £ fJ). The spaces of Lipschitzian cylinder functions on fi and 
are denoted respectively by 

L,p(r!) :={(^(Bi,,Bt,,--- ,StJ :Vn> l,ti,-- - , i„ € [0, oo), G Q.L,p(IR'*><")}, 

L,p(0) :={(^(Bt,,^t,,--- ,StJ :V7i> l,ti,-- - , i„ G [0, (»), G C,.L.p(K'^><")}. 

Let G'(-) : §((i) — > M be a given continuous monotonic and sublinear function. 
Following Sec. 2 in Chap Hill we can construct the corresponding G-expectation 
E on (yi^Liyf^t)). Due to the natural correspondence of Lip(f2) and Lip(O), we 
also construct a sublinear expectation E on (f2, L,ip(r2)) such that (i3t((I'))t>o is 
a G-Brownian motion. 

The main objective of this section is to find a weakly compact family of (cr- 
additive) probability measures on (fi, Bi^l)) to represent G-expectation E. The 
following lemmas are a variety of Lemma 13.31 and 13.41 

Lemma 2.1 Let <ti < h <_■■■< < oo and C Gi.Lip(M'^^™) 

satisfy ipn i 0. Then E[(p„(St^ , , • • • , Bt^)] | 0. 

We denote T := {i = (ii, . . . ,im) : Vm G N, < ti < i2 < ■ ■ • < im < oo}. 

Lemma 2.2 Let E be a finitely additive linear expectation dominated by E on 
Lip{ri). Then there exists a unique probability measure Q on (ri,B{Q)) such 
that E[X] = Eq[X] for each X G L,pln). 

Proof. For each fixed t= {ti,. . . ,tm) G T, by Lemma [2.11 for each sequence 
Wn}^=i C Cl.L^p{R''''"~) satisfying i 0, we have E[ipn{Bt, , Bt,, ■ ■ ■ , BtJ] i 
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0. By Daniell-Stone's theorem (see Appendix B), there exists a unique probabil- 
ity measure on (R'^><™, S(R'*^'")) such that Eq^[i^] = E[ip{Bt, ,Bt^,--- , Bt,J] 
for each (p e C/.Lip(M'^^™). Thus we get a family of finite dimensional distribu- 
tions {Qt : t £ T}. It is easy to check that {Qt : i S T} is consistent. Then 
by Kolmogorov's consistent theorem, there exists a probability measure Q on 
{fl,B{Cl)) such that {Qt : t e T} is the finite dimensional distributions of Q. 
Assume that there exists another probability measure Q satisfying the condi- 
tion, by Daniell-Stone's theorem, Q and Q have the same finite-dimensional 
distributions. Then by monotone class theorem, Q = Q. The proof is complete. 
□ 

Lemma 2.3 There exists a family of probability measures Ve on (f2,S(f2)) such 
that 

E[X] = max Eq[X], for X e L^p{n). 

Proof. By the representation theorem of sublinear expectation and Lemma[221 
it is easy to get the result. □ 

For this Ve, we define the associated capacity: 

c{A) := sup Q{A), AeB{n), 

and the upper expectation for each ;B(i7)-measurable real function X which 
makes the following definition meaningful: 

E[X]:= sup Eq[X]. 
Qen 

Theorem 2.4 For {B)t>o , there exists a continuous modification {B)t>o of B 
(i.e., c{{Bt ^ Bt}) = 0, for each t>0) such that Bq = 0. 

Proof. By Lemma [2731 we know that E = E on Lip{rt). On the other hand, we 
have 

E[\Bt-Bs\^] ^E[\Bt~Bs\^] ^d\t~s\^ for s,ie [0,oo), 

where d is a constant depending only on G. By Theorem ll.36[ there exists a 
continuous modification B of B. Since c{{Bq ^ 0}) = 0, we can set _Bo = 0. 
The proof is complete. □ 

For each Q e 7?, let Q o B^^ denote the probability measure on (17, S(f2)) 
induced by B with respect to Q. We denote V\ = {Q o B~^ : Q e 7?}. By 
Lemma [2741 we get 

E[\Bt - B,\^] ^E[\Bt ~ Bs\^] ^ d\t - s\\ys,t e [0,oo). 

Applying the well-known result of moment criterion for tightness of Kolmogorov- 
Chentsov's type (see Appendix B), we conclude that Vi is tight. We denote by 
V = Vi the closure of Vi under the topology of weak convergence, then V is 
weakly compact. 

Now, we give the representation of G-expectation. 
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Theorem 2.5 For each continuous monotonic and sublinear function G : §((i) — >■ 
M, /ei E be the corresponding G-expectation on Lip{Q)). Then there exists a 
weakly compact family of probability measures V on {U,,B{^)) such that 

t[X] = max£;p[X] for X G Lip{Q.). 

Proof. By Lemma 12.31 and Lemma 12. 4[ we have 

t[X] = max Ep[X] for X G Up{VL). 

For each X G L,p{n), by Lemma O we get W(\X - {X A N) V {-N)\] i as 
iV — T' oo. Noting that V = Vi, by the definition of weak convergence, we get 
the resuh. □ 

Remark 2.6 In fact, we can construct the family V in a more explicit way: Let 
{Wt)t>Q = (W7)iLi t>o ^6 0, d-dimensional Brownian motion in this space. The 
filtration generated by W is denoted by JF^ . Now let T be the bounded, closed 
and convex subset in M''^'^ such that 

G{A) = sup M^77^], A G S{d), 
■yer 

(see U.lSp in Ch. II) and At the collection of all Q-valued {J-^)t>Q- adapted 
process [0,oo). We denote 

B] / -isdWs, t>0, jeAr- 
Jo 

and Vo the collection of probability measures on the canonical space {^,B{Q)) 
induced by {B'^ : 7 G ^r}- Then V = Vo (see ^37] for details). 



§3 G-capacity and Paths of G-Brownian Motion 



According to Theorem 12.51 we obtain a weakly compact family of probability 
measures V on {Q,,B{Vl)) to represent G-expectation ]£[•]. For this V, we define 
the associated G-capacity: 

c{A) := sup P{A), A G B{Vt) 

Pev 

and upper expectation for each X G LP{^) which makes the following definition 
meaningful: 

E[X] := sup Ep[X]. 
Pev 

By Theorem 12.51 we know that E = E on Lip{il), thus the E[| • |]-completion 
and the E[| • |]-completion of Lip{n) are the same. 
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For each T > 0, we also denote by fir = Cq{[0,T]) equipped with the distance 

We now prove that Lq{^) = Lj, where is defined in Sec.l. First, we need 
the following classical approximation lemma. 

Lemma 3.1 For each X e Cb{^) and n = 1, 2, • • • , we denote 

:= ij^i^{X{uj') +n\\u- c^'||c,^([o,„])} for co e n. 

Then the sequence 

{^'"^}^=i satisfies: 

(i) -M < < <---<X,M = sup^^a I^MI; 

(ii) |X(")(a;i)-X(")(w2)| <n||a;i-a;2||c„.([o,„]) for , ^2 e fi; 

(iii) X(")(a;)tX(a;) for w e fi. 

Proof, (i) is obvious. 
For (ii), we have 

X(")(wi)-X(")(a;2) 
< nllwi -W2||c„^([o,„]) 

and, symmetrically, X^"'\uj2) — X^'^\uJi) < n\\uJi — ^2\\c^(^[o Thus (ii) fol- 
lows. 

We now prove (iii). For each fixed w € fi, let w„ G f2 be such that 

X{un) +n\\u)- a;„||c„^([o,„]) < X^^\u>) + ^. 

It is clear that n \\uj - Wn||c^([o,„]) < 2M+1 or ||a; - uJn\\cd,^[o,n]) ^ ^^IT^- ^^^'^^ 
X e Cb{fl), we get X{iOn) — X{u)) as n — )• oo. We have 

X{oj) > X(")(w) > X{oJn) + n\\uj- a;„||c„^([o,„]) - K 
thus ^ 

n llw - (^n\\cg{lO,n]) < - + -• 

We also have 

X(u;„) - X(c.) + n lie. - c.n|lc„-([o,„]) > X^^Hco) - X(u;) 

> X(w„) - X(w) + n ||W - Wn|lc„^([o,„]) - ^- 
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From the above two relations, we obtain 

< 2{\X{LUn) - X{uj)\ + -) ^ as n -> oo. 
n 

Thus (iii) is obtained. □ 

Proposition 3.2 For each X E Cf,(ri) and e > 0, there exists Y G Lip{il.) such 
thatE[\Y ~X\] < e. 

Proof. We denote M — sup^^gfj By Theorem 11.121 and Lemma |3.1[ 

we can find > 0, T > and X e Cbinr) such that E[|X - X\] < e/3, 
sup^gn < M and 

|X(a;) — X{uj')\ < fi\\uj ~ ^'\\ci^{[o t]} ^'^^ w,a;' e fl. 
Now for each positive integer n, we introduce a mapping cj 

u;(")(c.)(t) = ^ -j±:^^±lL-[(t-^^ - t)u{tl) + {t~ tl)^{tl^,)] + l[T,oo)(iVW, 

where = ^, fc 0, 1, • • • , n. We set := X(a;(") (w)), then 

<M sup |w(")(w)(t)-cj(")(cj')(i)| 

t6[0,T] 

= /z sup \u{1^{,)-uj'{tl)\. 

We now choose a compact subset iiT C f2 such that ^[l^^-c] < e/QM. Since 
sup^g^ supjg[o,T] -> 0, as n — cx), we can choose a sufRciently 

large rtg such that 

sup \X{lu) ~ X("«)(a;)| = sup \X{uj) ~ X(cj("°)(w))| 
<^sup sup |a;(t) - 
< e/3. 

Set y it follows that 

E[|^->"|] <]E[|X-X|]+E[|X-X("'')|] 

< E[|X - X|] + ¥.[1k\X - + 2ME[l^c] 

< £. 

The proof is complete. □ 
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By Proposition 13.21 we can easily get L^{fl) = L;!,. Furthermore, we can get 
LP.{n) = HP, yp > 0. 

Thus, we obtain a pathwise description of Lq^^I) for each p > 0: 

L^{n) ~ {X e L^{n) : X has a quasi-continuous version and hm ]E[|X|P/r|x|>,i>] — 0}. 

Furthermore, E[X] = E[X], for each X e 
Exercise 3.3 Show that, for each p > 0, 

L^{Qt) — {X £ L'^{Qt) ■ X has a quasi- continuous version and hm ]E[|X|P/r|j(-|>„i ] — 0}. 



Notes and Comments 

The results of this chapter for G-Brownian motions were mainly obtained by 
Denis, Hu and Peng (2008) ;33] (see also Denis and Martini (2006) [si] and 
the related comments after Chapter 111). Hu and Peng (2009) [Hsl then have 
introduced an intrinsic and simple approach. This approach can be regarded 
as a combination and extension of the original Brownian motion construction 
approach of Kolmogorov (for more general stochastic processes) and a sort of 
cylinder Ligschitz functions technique already introduced in Chap. Illll Section 
1 is from |37| and Theorem 12.51 is firstly obtained in [s^] , whereas contents of 
Sections 2 and 3 are mainly from [58]. 

Choquet capacity was first introduced by Choquet (1953) jl^l, see also Del- 
lacherie (1972) 32] and the references therein for more properties. The ca- 
pacitability of Choquet capacity was first studied by Choquet_|24| under 2- 
alternating case, see Dellacherie and Meyer (1978 and 1982) [33|, Huber and 
Strassen (1972) [62] and the references therein for more general case. It seems 
that the notion of upper expectations was first discussed by Huber (1981) [6l| 
in robust statistics. Recently, it was rediscovered in mathematical finance, es- 
pecially in risk measure, see Delbaen (1992, 2002) [sj,!!^, FoUmer and Schied 
(2002, 2004) ^5^ and etc.. 



Appendix A 

Preliminaries in Functional 
Analysis 



§1 Completion of Normed Linear Spaces 

In this section, we suppose "H is a linear space under the norm || • ||. 

Definition 1.1 {x„} E H is a Cauchy sequence, if {xn} satisfies Cauchy's 
convergence condition: 

lim \\Xn - Xm\\ = 0. 

Definition 1.2 A normed linear space H is called a Banach space if it is 
complete, i.e., if every Cauchy sequence {x„} of H converges strongly to a 
point Xqo of H : 

lim \\x„ - a;oo|| = 0. 

n— fcxD 

Such a limit point Xoo, if exists, is uniquely determined because of the triangle 
inequality \\x — a;'|| < \\x ~ x„|| + ||a;„ — a;'|| . 

The completeness of a Banach space plays an important role in functional anal- 
ysis. We introduce the following theorem of completion. 

Theorem 1.3 Let % be a normed linear space which is not complete. Then % 
is isomorphic and isometric to a dense linear subspace of a Banach- space %, 
i.e., there exists a one-to-one correspondence cc i of T-L onto a dense linear 
subspace of % such that 

x + y — x-\-y, ax — ax, ||x|| = ||a;|| . 

The space % is uniquely determined up to isometric isomorphism. 

For a proof see Yosida [HI (1980, p.56). 
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§2 The Hahn-Banach Extension Theorem 

Definition 2.1 Let Ti and T2 be two linear operators with domains D{Ti) and 
D{T2) both contained in a linear space %, and the ranges R{Ti) and R{T2) both 
contained in a linear space Ai. Then Ti = T2 if and only if D{Ti) ~ D[T2) 
and Tix — T2X for all x € D{Ti). If D{Ti) C D(T2) and Tix — T2X for all 
X G D{Ti), then T2 is called an extension of Ti, or Ti is called a restriction 

0fT2. 

Theorem 2.2 (Hahn-Banach extension theorem in real linear spaces) 

Let % be a real linear space and let p{x) be a real-valued function defined on H 
satisfying the following conditions: 

p{x + y) < p{x) + p{y) (subadditivity); 

p{ax) — ap{x) for a > (positive homogeneity). 

Let C be a real linear subspace of % and fo be a real-valued linear functional 
defined on C : 

foiax + I3y) = afo{x) + (ifo{v) for x,y e C and a, /3 G M. 

Let fo satisfy fa{x) < p{x) on C Then there exists a real-valued linear func- 
tional F defined on % such that 

(i) F is an extension of fo, i.e., F{x) = fo{x) for all x € £. 

(ii) F{x) <p{x) forx G H. 



For a proof see Yosida [125| (1980, p.l02). 

Theorem 2.3 (Hahn-Banach extension theorem in normed linear spaces) 

Let % be a normed linear space under the norm \\ -W, C be a linear subspace of 
T-L and let fi be a continuous linear functional defined on C Then there exists 
a continuous linear functional f , defined on %, such that 

(i) / is an extension of fi . 

(ii) ll/ill = 11/11- 



For a proof see for example Yosida |125{ (1980, p. 106) 



§3 Dini's Theorem and Tietze's Extension The- 
orem 

Theorem 3.1 (Dini's theorem) Let T-L be a compact topological space. If a 
monotone sequence of bounded continuous functions converges pointwise to a 
continuous function, then it also converges uniformly. 

Theorem 3.2 (Tietze's extension theorem) Let C be a closed subset of a 
normal space % and let f : C ^ M. be a continuous function. Then there exists 
a continuous extension of f to all ofH with values in M. 



Appendix B 

Preliminaries in Probability 
Theory 



§1 Kolmogorov's Extension Theorem 

Let X be a random variable with values in M" defined on a probability space 

{n, T, P). Denote by B the Borcl cr-algcbra on K". Wc define X's law of distri- 
bution Px and its expectation Ep with respect to P as follows respectively: 



where B G B. 

In fact, we have Px{B) = Ep[1b{X)]. 

Now let {Xt}teT be a stochastic process with values in R" defined on a prob- 
ability space {fl,T,P), where the parameter space T is usually the halfiine 
[0,-hoo). 

Definition 1.1 The finite dimensional distributions of the process {Xt}t^T 
are the measures Uti,-- M defined on R"'^, /c = 1, 2, • • • , by 

fit,,- ,t, (Bi X • • • X Bk) := P[Xt, GBi,--- ,Xt,€ Bk], ti € T, i = 1,2,- ■■ ,k, 

where Bi G B,i = 1,2, - ■■ ,k. 

The family of all finite-dimensional distributions determine many (but not all) 
important properties of the process {Xt}teT- 

Conversely, given a family {t'd,... ,tk ■ ti € T,i = 1,2, - ■ ■ , fc. A; G N} of probabil- 
ity measures on M"*^ , it is important to be able to construct a stochastic process 
{Yt)t£T with i'ti,---,tk as its finite-dimensional distributions. The following fa- 
mous theorem states that this can be done provided that {i^ti,--- ,tk} satisfy two 
natural consistency conditions. 



Px{B) := P(w : X{ij) e B); Ep[X] := 
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Theorem 1.2 (Kolmogorov's extension theorem) For all ti, t2, ■■■ , tk, 

fc S N, let iyti,---,tk be probability measures on M"*^ such that 

i't„(i),---,t„(fc)(-Bi X • • • X Bk) = lyti,--- ,tkiBTr-i{i) X • • • X 
for all permutations tt on {1, 2, • • • , fc} and 

ut,,.. {BiX---xBk) = ut,,- MM+i,- .t.+™ (Bi X • • • X Bfc X M" X • • • X M") 

for all m G N, where the set on the right hand side has a total of k + m factors. 
Then there exists a probability space (f2, P) and a stochastic process {X^) on 
n, Xt:n~>-W^, such that 

yt,,-M{Bi X ••• X Bk) = P[Xt, eBi,--- ,Xt,e Bk] 

for all ti G T and all Borel sets Bi, i = 1, 2, • • • , fc, fc G N. 

For a proof see Kolmogorov f?! (1956, p.29). 

§2 Kolmogorov's Criterion 

Definition 2.1 Suppose that (Xt) and (Yf) are two stochastic processes defined 
on (f2,J^, P). Then we say that (Xt) is a version of (or a modification of) 

(Yt) tf 

: Xtiuj) = Yt{uj)}) = 1 for all t. 

Theorem 2.2 (Kolmogorov's continuity criterion) Suppose that the pro- 
cess X — {Xt}t>o satisfies the following condition: for all T > there exist 
positive constants a, P, D such that 

E[\Xt~ X,\'']<D\t^ s\^+f', 0<s,t<T. 

Then there exists a continuous version of X . 



For a proof see Stroock and Varadhan jll7l | (1979, p. 51) 



Let i5 be a metric space and B be the Borel ti-algebra on E. We recall a few 
facts about the weak convergence of probability measures on {E,B). If P is 
such a measure, we say that a subset ^ of is a P-continuity set if P{dA) = 0, 
where dA is the boundary of A. 

Proposition 2.3 For probability measures Pn{n G N) and P, the following 
conditions are equivalent: 

(i) For every bounded continuous function f on E, 

lim / /dP„ = / fdP- 



§2 Kolmogorov's Criterion 



HI 



(ii) For every hounded uniformly continuous Junction f on E, 



lim / fdP„ - / fdP; 




(iii) For every closed subset F of E, lira sup, 

(iv) For every open subset G of E, liminf„_ 



ooPn{F)<P{F); 
Pn(G) > P{G); 



(v) For every P -continuity set A, lim„^oo Pn{A) ~ -P(^)- 

Definition 2.4 // P„ and P satisfy the equivalent conditions of the preceding 
proposition, we say that (P„) converges weakly to P. 

Now let TT be a family of probability measures on {E,B). 

Definition 2.5 A family n is weakly relatively compact if every sequence 
of TT contains a weakly convergent subsequence. 

Definition 2.6 A family tt is tight if for every e € (0, 1), there exists a compact 
set such that 



With this definition, we have the following theorem. 

Theorem 2.7 (Prokhorov's criterion) If a family n is tight, then it is weakly 
relatively compact. If E is a Polish space (i.e., a separable completely metrizable 
topological space), then a weakly relatively compact family is tight. 

Definition 2.8 // (X„)„gN o,nd X are random variables taking their values in 
a metric space E, we say that (X„) converges in distribution or converges 
in law to X if their laws Px„ converge weakly to the law Px of X . 

We stress the fact that the (X„) and X need not be defined on the same prob- 
ability space. 

Theorem 2.9 (Kolmogorov's criterion for weak compactness) Let{X"} 
be a sequence of M.'^-valued continuous processes defined on probability spaces 
(r2", such that 

(i) the family {-Pjjn} of initial laws is tight in R*^. 

(ii) there exist three .strictly positive constants a, l3, 7 .such that for each s, i G 1R+ 
and each n, 



then the set of the laws of the {Xn) is weakly relatively compact. 



P{K^) > 1 — e for every P € tt. 



Ep.[\Xl^-XlT]<l3\s^t\ 



For the proof see Revuz and Yor (1999, p. 517) 
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§3 Daniell-Stone Theorem 

Let (fi, T ^ /i) be a measure space, on which we can define integration. One 
essential properties of integration is its linearity, thus it can be seen as a lin- 
ear functional on L^(i7, J^, /i). This idea leads to another approach to define 
integral-Daniell's integral. 

Definition 3.1 Let $7 he an abstract set and H be a linear space formed by a 
family of real valued functions. % is called a vector lattice if 

feH^lfleUjAleU. 

Definition 3.2 Suppose that H is a vector lattice on and I is a positive linear 
functional onH, i.e., 

f,geH,a,^eR^ I{af + (ig) = «/(/) + pi{g); 
f eHJ>O^Iif)>0. 

If I satisfies the following condition: 

/„ eH,/„iO^/(/„) ^0, 

or equivalently, 

fn eH,fnt.f^H^ I{f) = lim /(/„), 

n— f oo 

then I is called a Daniell's integral on %. 

Theorem 3.3 (Daniell-Stone theorem) Suppose that H is a vector lattice 
on fl and I is a Daniell integral on %. Then there exists a measure fi £ J- , 
where T := a{f : f e U), such thatnc L^{n,T,fi) and I{f) = fi{f), V/ G n. 
Furthermore, «/ 1 G where H*^ :— {/ : 3/„ > 0, /„ G H such that fn t /}; 
then this measure ii is unique and is a-flnite. 

For the proof see Dellacherie and Meyer [sS] (1978, p.59), Dudley [4l| (1995, 
p.142), or Yan (1998, p.74). 



Appendix C 
Solutions of Parabolic 
Partial Differential 
Equation 



§1 The Definition of Viscosity Solutions 

The notion of viscosity solutions was firstly introduced by Crandall and Lions 
0981) 128| and (1983) [gg (see also Evans's contribution (1978) ^] and (1980) 
46 1) for the first-order Hamilton- Jacobi equation, with uniqueness proof given in 
2^ . The the proof of second-order case for Hamilton- Jacobi-Bellman equations 
was firstly developed by Lions (1982) [s^] and (1983) |81|] using stochastic control 
verification arguments. A breakthrough was achieved in the second-order PDE 
theory by Jensen (1988) ^6^. For all other important contributions in the 
developments of this theory we refer to the well-known user's guide by Crandall, 
Ishii and Lions (1992) For reader's convenience, we systematically interpret 
some parts of [3u| required in this book into it's parabolic version. However, 
up to my knowledge, the presentation and the related proof for the domination 
theorems seems to be a new generalization of the maximum principle presented 
in ^3d\. Books on this theory are, among others, Barles (1994) Fleming, and 
Soner (1992) ^a], Yong and Zhou (1999) :i24j. 
Let r > be fixed and let O c [0, T] x R^. We set 

USC{0) — {upper semicontinuous functions m : O — > K}, 

LSC{0) — {lower semicontinuous functions u : O — > K.}. 
Consider the following parabolic PDE: 

(E) dtu - G{t, X, u, Du, D^u) = on (0, T) x M^, , , 

(IC) u{0, x) = tp{x) for X e R^, ^ 

where G : [0, T] x R^ x R x R^ x E,{N) R, Lp e C(R^). We always suppose 
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that G is continuous and satisfies the following degenerate elliptic condition: 

G{t,x,r,p,X)>G{t,x,r,p,Y) whenever X > F. (1.2) 

Next we recall the definition of viscosity solutions from Crandall, Ishii and Lions 
[sBl. Let u : (0,r)xM^ ^ M and {t,x) G (0,r)xR^. We denote hyV'^'^ u[t,x) 
(the "parabolic superjet" of u at (<, x)) the set of triples (a,p, X) G M x x 
§(iV) such that 

u{s, y) < u(t, x) + a(s — t) + {p.y — x) 

+ \{X{y ^x),y-x)+ o{\s -~t\ + \y- xp). 

We define 

P^'+uit,x) ■.^{{a,p,X) e M X X S{N) : 3(t„, x„, a„,p„, X„) 
such that (a„,p„,X„) G V'^'^uitn, Xn) and 

{tn,Xn,u{tn,Xn),an,Pn,Xn) -> {t,X,u{t,x),a,p,X)}. 

Similarly, we define 'P'^'^u{t, x) (the "parabolic subjet" of u at {t,x)) by 
V'^'^u{t,x) := -V'^^+{-u){t,x) andV'^^-u{t,x) hy'P'^^'u{t,x) := -V^'+{-u){t,x). 

Definition 1.1 (i) A viscosity subsolution of (E) on (0, T) xR^ is a function 
u G USG{{0,T) X M^) such that for each {t,x) G (0,r) x M^, 

a- G{t,x,u{t,x),p,X) < for {a,p,X) G V^-+u{t,x); 

likewise, a viscosity supersolution of (E) on (0,T) X R^ IS a function v G 
LSG{{0,T) X M^) such that for each {t,x) G (0,T) x M^, 

a - G'(i,a;,w(t,a;),p,X) > for {a,p,X) G ■p^'^v{t,x); 

and a viscosity solution of (E) on (0,r) X R^ is a function that is simultane- 
ously a viscosity subsolution and a viscosity supersolution of (E) on (0, T) x R^. 
(ii) A function u G U SC{\Q,T) x R^) is called a viscosity subsolution of 
(fTTI) on [0, T) X R^ if u is a viscosity subsolution of (E) on (0,T) X and 
M(0,a;) < f{x) for x G R^; t/ie appropriate notions of a viscosity supersolution 
and a viscosity solution of (11.11) on [0,r) X R^ are then obvious. 

We now give the following equivalent definition (see Crandall, Ishii and Lions 
0). 

Definition 1.2 A viscosity subsolution of (E), or G- subsolution, on(0,T)xR^ 
IS a function u G USG{{0,T) x R^) smc/i that for all {t,x) G (0,T) x R^, 
G C2((0, T) X K^) suc/i that u{t, x) = <l>{t, x) and u < (j) on (0, T) x R^\{t, x), 
we have 

dt<i){t, x) - G{t, x, (j){t, x), D4){t, x),D'^4){t, x)) < 0; 
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likewise, a viscosity supersolution of (E), or G- super solution, on (0,T) x 

is a function v G LSC{{Q,T) x R^) such that for all [t,x) S (0,r) x R^, 

cj) e C2((0, T) X R^) such that u{t, x) = <i)(t, x) anduxj) on (0, T) x R^\(i, x), 

dtcj){t, x) - G{t, X, (pit, x), D(l){t, x),D^(f){t, x)) > 0; 

and a viscosity solution of (E) on (0, T) x R^ is a function that is simultaneously 
a viscosity subsolution and a viscosity supersolution of (E) on (0, T) x R^. The 
definition of a viscosity solution of (jl.ip on [Q,T) xR^ is the same as the above 
definition. 



§2 Comparison Theorem 

We will use the following well-known result in viscosity solution theory (see 
Theorem 8.3 of Crandall, Ishii and Lions [30|). 

Theorem 2.1 Let u, eUSC{{0,T) x R^') for i ^ 1, ■ ■ ■ ,k. Let ip be a func- 
tion defined on (0, T) x R^^"' such that {t,xi, . . . , Xk) — > ^p{t, xi, . . . , Xk) 
is once continuously difjerentiable in t and twice continuously differentiable in 
(xi, • • • ,Xfe) e R^i+'-'+^s^. Suppose that i £ (0,r), € R^' /or i = 1, • • • , fc 
and 



w{t,xi,--- ,Xk) := ui{t,xi)-\ \-Ukit,Xk) - ip{t,xi, - ■ ■ ,Xk) 

< w(i,xi, ■■ ■ ,Xk) 

for t e (OjT) and Xi £ R^' . Assume, moreover, that there exists r > such 
that for every M > there exists constant C such that for i — 1, ■ ■ ■ ,k. 



bi < C whenever {bi,qi,Xi) G V'^'^Ui{t, Xi), 

\x,-Xi\ + \t~i\<r and \ui{t,Xi)\ + \q,\ + \\Xi\\ < M. 

Then for each e > 0, there exist Xi £ S{Ni) such that 

{i) {b.„Dx^ip{t,xi,- ■ ■ ,Xk),Xi) £V ' Ui{t,Xi), i = l,---,k, 

(ii) 

~ Xi ■■■ 



{- + \\A\\)I< 







Xk 



<A + eA^ 



(iii) &i H h 6fe = dt(p{i,xi, ■ ■■ ,Xk), 

where A = Dlip{i, x) £ S{Ni H \-Nk). 



(2.3) 



Observe that the above condition (12.31) will be guaranteed by having each Ui be 
a subsolution of a parabolic equation given in the following two theorems. 
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In this section we will give comparison theorem for G-solutions with different 

functions G. 

(G) We assume that 

G : [0, T] X X M X X §{N) ^ R, i = 1, • • • , fc, 

are continuous in the following sense: for each t £ [Q,T), v £ R, x, y, p £ 
andXGS(Af), 

\Gi{t,x,v,p,X) - Gi{t,y,v,p,X)\ 

< w(l + (T - t)-' + \x\ + \y\ + \v\M\x -y\ + \p\ ■ \x - y\), 

where oj, w : M+ — >■ M+ are given continuous functions with ui{0) = 0. 

Theorem 2.2 (Domination Theorem) We are given constants /3j > 0, i = 
I,-- - ,k. LetUiG USCi[0, T]xR^) be subsolutions of 

dtu- Gi{t,x,u,Du,D^u) = 0, i = l,---,k, (2.4) 

on (0,T) X such that (^i^i PiUi{t,x)j 0, uniformly as \x\ oo. We 

assume that the functions {Gi}i = l'^ satisfies assumption (G) and the following 
domination condition holds for : 

k 

Y,PiGi{t,x,Vi,Pi,Xi)<Q, (2.5) 

i=l 

for each {t,x) e (0,T) x and {vi,pi,Xi) e M x x S{N) such that 

E-=i Pi^i > o> Eti ftK = 0, Eti ft^. < 0. 

Then a domination also holds for the solutions: if the sum of initial values 
EiLi /3jUi(0, •) is a non-positive function on , then Ei=i /3i'"i(^) < 0, for 
all t > 0. 

Proof. We first observe that for ^ > and for each 1 < i < fc, the functions 

defined by iii := Ui — d/{T — t) is a subsolution of 

dtUi - Gi{t, x, Ui, Dili, D'^Ui) < - 



(T - 1)2 ' 



where Gi{t,x,v,p,X) := Gi{t,x,v + S/(T — t),p,X). It is easy to check that 
the functions Gj satisfy the same conditions as Gj . Since 

Ei=i A^i < follows 

from Ei=2 < in the limit 5 4 0, it suffices to prove the theorem under the 
additional assumptions: 



dfUi — Gi{t,x,Ui,Dui, D'^Ui) < — c, where c = S/T^, 
and limt-^TU^{t,x) = — oo uniformly on [0, T) x 



,N (2.6) 
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To prove the theorem, we assume to the contrary that 

k 

sup N, x) — mo > 

(t,2;)e[0,T)xR" 

We will apply Theorem [2l1 for x = {xi, ■ ■ ■ ,Xfc) G M'^'^^ and 



For each large a > 0, the maximum of it; — (^q achieves at some (i", x") inside 
a compact subset of [0,T) x M.'^'^'^ . Indeed, since 

k 

M„ = ^Au,(r,xn-'^a(x") >mo, 

we conclude must be inside an interval [0, Tq], Tq < T and must be inside 
a compact set {x e M*^^^ : suptg[Qj^jj] w{t,x) > We can check that (see 

Isdl Lemma 3.1) 



(i) lim^^oc Vaix") = 0, 

(ii) limQ^ooA^a = lima_j.oo x?) H h /3A;Ufe(i", x^)) 

= sup(t_^|g[o.T)xK« [/3iui(t,a;) H h^feUfc(i,x)] 

= [/?iui(£,x) H |-/3fcUfe(i, x)] = mo. 



(2.7) 



where (t, x) is a limit point of (t", x"). Since Ui £ USC, for sufficiently large a, 
we have 

/3iui(t",x?) + • • • + /3,Ufe(r,x^) > ^. 

If £ = 0, we have limsup„_^^ X^iLi ) = I]*Li /3iUi(0, x) < 0. We 
know that t > and thus must be strictly positive for large a. It follows 
from Theorem [O that, for each e > there exist bf eR, X, € S{N) such that 



(6^/3-ii?,.^(x"),X,) e7^2,+^,(r,x,n, ^A6? = 0for* = l,... ,fc, (2. 
and such that 



i=l 



-(- + PII)/< 



f PiXi ... Q \ 

... Pk-iXk-i 

V ... PkXk J 



<A + eA^, (2.9) 
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where A = D'^Lpa{x°') G S{kN) is explicitly given by 



A = aJkN, where JkN = 



( In —In 
—In "21 n 



\ 



\ 



2In —In 
—In In / 



The second inequahty of p.9p impUcs l^i^i 1^ 0- Set 



Thus PiPf = 0- From this together with and dM]), it follows that 

By ((2J)) (i), we also have liniQ^oo bf | • |a;f - | ^ 0. This, together with the 
domination condition (|2.5p of G,;, implies 



i=l i— 1 i— 1 i— 1 

k 

> -^/3,G,(r,a;?,w,r,xn,P?,^») 
1=1 

- ^ A|G,(i", x^ xn,P?, - G,r , X?, , xn,Pr, 



4=1 



> -Y^PM^ + {T- To)-' + |x?| + Ixfl + H{e,xf)\) ■ u;{\xf - x?| 



z=l 



The right side tends to zero as a — >■ oo, which induces a contradiction. The 
proof is complete. □ 

Theorem 2.3 We assume that the functions Gi — Gi{t,x,v,p, X), i = 1, • • • ,k 
and G = G{t,x,v,p, X) satisfy assumption (G). We also assume G dominates 
{Gi}^^i in the following sense: for each {t,x) € [0, oo) x M.^ and {vi,pi, Xi) G 
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M X X S{N), 

J2G^{t,x,v„p„X,) < G(i,x,^w„^K,^X,). (2.10) 

i— 1 i—l i—1 i—1 

Moreover there exists a constant C such that, for each {t,r,x,p) G [0,oo) x 
M^xM^ and, Yx,Yi^ ^{N) such that Y2 > Yi 

G{t, X, r, p, Fa) - G{t, x,ri,p,Yi)>0 

and there exists a constant C such that 

\G{t, X, r,p, Xi) - G{t, X, r,p, Xa)! < C(|ri - ral + |Xi - X2I). 

Let Ui eUSCi[0,T] x R^) be a G, -sub solution and u eL5C([0,T] x M^) he a 
G-supersolution such that Ui and u satisfy polynomial growth condition. Then 
X^iLi ^) — u{t,x) on [0,T) X provided that X]i=i'"i|t=o ^ '"|t=o- 

Proof. For a fixed and large constant A > (7 + C we set ^{x) := (1 + |a;p)'/^ 
and 

u.i{t,x) -.^ u.i{t,x)£,{x)e~^\ i = l,-- - ,fc, Uk+i{t,x) := -u{t,x)£,{x)e~^\ 

where I is chosen to be large enough such that ^ |ui(t, a;)| — > uniformly as 
|x| — > 00. It is easy to check that, for each i = l,-- - ,fc + l, {tiisa subsolution 
of 

dtUi + Xui - Gi{t, X, Ui,Dui, D^Ui) = 0, 
where, for each i = 1, • • • , fc, the function Gi{t, x, v,p, X) is given by 

e-^*r^G,;(t, X, e^*<, e^\p + vD^),e^\X£, +p^D(,+D^^p + vD^S)), 

and Gfc+i(t, x, v,p, X) is given by 

-e-^*r^G(t, x, -e^*<, -e^*(p + vD^, -e^\X^+p(g) D^ + D^(g)p + vD^^))- 
Observe that 

Dax)^iax){l + \x\^)-'x, 

D^^ix) = ^ix)[l{l + \x\^)-^I + lil - 2)(1 + ® x]. 

Thus both ^^^(a::)|_D^(a;)| and ^^^(x)|-D^^(a;)| converges to zero uniformly as 
\x\ — > 00. 

From the domination condition (|2.10l) . for each {vi,pi,Xi) G R x R^ x S(A^), 

i = 1, • • • , fc + 1, such that Ya^I «j = 0> Ya^I Pi = 0, and X]i=^i = 0^ 
have 

fc+i 

Gi{t,x,Vi,pt,Xi) < 0. 
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For v,r eR,peR^ and X,Y e S{N) such that r > 0, y > and r > 0, smce 
G is still monotone in X, 

Gk+i{t,x,v,p,X) - Gk+i{t,x,v - r,p,X + Y) 

< Gk+i{t,x,v,p,X) - Gk+i{t,x,v - r,p,X) 

< c{C + Ci)r, 

where the constant Ci does not depend on {t,x,v,p,X). We then apply the 
above theorem by choosing f3i = 1, i — 1, ■ ■ ■ , k + 1. Thus ^i\t=o < 0. 

Moreover for each e M, p, e and X, e S{N) such that v = J^'^^^ > 0, 
Y.'i^i Pi = and X = Y.t=i Xi < we have 

fc+i fc+i 

,x,Vi,pi,Xi) 

i=l i=l 
k 

= -Aw + '^Gi{t,x,Vi,pi,Xi) + Gk+i{t,x,vk+i - v,pk+i,Xk+i - X) 
1=1 

+Gk+iit, X, Vk+i,Pk+i, Xk+i) - Gk+iit,x,Vk+i - v,pk+i,Xk+i - X) 

< -\V + Gk+lit,X,Vk+l:Pk+l,Xk+l) - Gfe+l(t, X, Vk+1 - V,Pk+l,Xk+l) 

< -Xv + c{C + C)v <0 

It follows that all conditions in Theorem 12.21 are satisfied. Thus we have 
E^=/ "^ < 0> or equivalently, J2i=i Ui{t,x) < u{t,x) for {t,x) € [0,r) x M^. □ 
The following comparison theorem is a direct consequence of the above domi- 
nation theorem. 

Theorem 2.4 (Comparison Theorem) We are given two functions G — G{t,x,v,p,X) 
and Gi = Gi{t,x,v,p^ X) satisfying condition (G). We also assume that, for 
each {t, X, v, p, X) € [0, oo) x x R x R^ and Y G S{N) such that X >Y , 

G{t,x,v,p,X) > Gi{t,x,v,p,X), (2.11) 
G{t,x,v,p,Y) > G{t,x,v,p,X). (2.12) 

We also assume that G is a uniform Lipschitz function in v and X , namely, for 
each {t,x) e [0,oo) x R^ and {v,p,X), {v',p,X') e R x R^ x S{N), 

\G{t,x,v,p,X) ~ G{t,x,v\pl ,X')\ < C{\v - v'\ + \X- X'\). 

Let ui eUSC{[0,T] x R^) be a Gi-suh solution and u &LSC{[0,T\ x R^) he 
a G- super solution on (0, T) x R^ satisfying the polynomial growth condition. 
Then u > ui on [0, T) X R^ provided that u\t=Q > iti|f=o- In particular this 
comparison holds for the case where G = Gi, which is a Lipschitz function in 
{v,X) and satisfies the elliptic condition \2.12]) . 

The following special case of the above domination theorem is also very useful. 
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Theorem 2.5 (Domination Theorem) We assume that Gi and G satisfy the 
same conditions given in the previous theorem except that the condition \2.11\) 
is replaced by the following one: for each {t,x) G [0,oo) x M.^ and {v,p,X), 
{v',p',X') e M X X S{N), 

Gi {t, X, v,p, X) - Gi {t, X, v',p', X') > G{t, x,v~v',p~ p', X - X') 

Let u eUSC{[0,T] x R^) be a Gi-subsolution and v eLSC{[0,T] x R^) be a 
G I- super solution on (0,T) x R^ and w is a G- super solution. They satisfy the 
polynomial growth condition. If [u — v)\t={) = w|t=o then u — v < w on [0, T) x 
R^. 

The following theorem will be frequently used in this book. Let G : R^ x §(iV) — !• 
R be a given continuous sublinear function monotonic in A e 'B{N). Obviously, 
G satisfies conditions (G) of Theorem l2.3l We consider the following G-equation: 

dtU-G{Du,D'^u) = Q, u{Q,x) ^ ip{x). (2.13) 

Theorem 2.6 Let G : R" x S(iV) ^ R 6e a given continuous sublinear function 
monotonic in A G S{N). Then we have 

(i) If u E USC([0,T] X R^) with polynomial growth is a viscosity subsolution 

of (I2.13P and v £ LSC([0,r] x R-'^) with polynomial growth is a viscosity 
super solution of (j2.13l) . then u < v. 

(ii) Ifu'^ e d{[0, T] xR^) denotes the polynomial growth solution of (12.13^ with 

initial condition (p, then u^"^ = Xu"^ for each A > and u^^^ < -\- u"^ . 

(iii) // a given function G : R"'^ x S(A^) i— > R is dominated by G, i.e. 

G{p,X)~G{p',X')<G{p-p',X~X'), forp,p'eR, X,X'eS{N), 

then for each ip G G(R^) satisfying polynomial growth condition, there 
exists a unique G-solution u'^{t,x) on [0,oo) x R^ with initial condition 
u'^\t=o = ip (see the next section for the proof of existence) , i.e., 

dtu'^ - G{Du'^, D^u'P) = 0, u'^|t=o = V- 

Moreover 

u'^{t,x)-u'''{t,x) <u'^-'^{t,x), fort>0, xeR^. 
Consequently, the following comparison holds: ij} > ip implies > u''^ . 



Proof. By the above corollaries, it is easy to obtain the results. 



□ 
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§3 Perron's Method and Existence 

The combination of Perron's method and viscosity solutions was introduced by 
H. Ishii [64]. For the convenience of readers, we interpret the proof provided in 
Crandah, Ishii and Lions [soj into its parabohc situation. 
We consider the foUowing parabohc PDE: 



where G : [0, cx)) x x M x x S{N) R, ip e C(R^). 

To discuss Perron's method, we wiU use the following notations: if m : O — > 

[— oo, oo] where O C [0, oo) x M^, then 

f u*{t,x) = limr4,oSup{u(s, y) : (s, y) G O and ^ J\s -t\ + \y - < r}, 
\ a;) = limr4.o inf{w(s, y) : (s, y) G O and ^J\s -t\ + \y- a^p < r}. 



One calls u* the upper semicontinuous envelope of u; it is the smallest upper 
semicontinuous function (with values in [—00,00]) satisfying u <u*. Similarly, 
u* is the lower semicontinuous envelope of u. 

Theorem 3.1 (Perron's Method) Let comparison hold for p.l4p . i.e., if w is 
a viscosity subsolution of p.l4p and v is a viscosity supersolution of p.l4p . 
then w < V. Suppose also that there is a viscosity .subsolution u and a viscosity 
supersolution u of (|3.14|) that satisfy the condition u^(0, a;) — u*[Q,x) — f^x) 
forxe R^. Then 

W(t,x) = sup{w(t,x) : u< w < u and w is a viscosity subsolution of (j3.14p } 
is a viscosity solution of p.l4p . 

The proof consists of two lemmas. For the proof of the following two lemmas, 
we also see flj. The first one is 

Lemma 3.2 Let J- be a family of viscosity subsolution of l|3.14p on (0,oo)xR^. 
Let w(t,x) = sup{u(i,x) : u € J-} and assume that w*{t,x) < 00 for (t,x) € 
(0,00) X R^. Then w* is a viscosity subsolution of (j3.14p on (0,oo) x R^. 

Proof. Let {t,x) e (0, 00) x R^ and consider a sequence s„, y„, Un (1 J- such 
that lim„^oo(sn, yn, w„(s„, y„)) = {t,x,w* {t,x)). There exists r > such that 
Nr = {(s, y) e (0, 00) X R^ : y/\s -t\ + \y- x^ < r} is compact. For </> e 
such that (j){t,x) = w*{t,x) and w* < (j) on. (0,oo) x R^\{t,x), let (t„,x„) be a 
maximum point of u„— over Nr, hence Un{s, y) < Unitn, Xn)+(t>is, y)—(f>{tn, Xn) 
for (s, y) G Nr. Suppose that (passing to a subsequence if necessary) (i„, Xn) 
{t,x) as n — >■ 00. Putting (s, y) = (s„,y„) in the above inequality and taking 
the limit inferior as n — ^ 00, we obtain 



{ 



dtu — G{t, X, u, Du, D^u) = on (0, 00) x R 
u(0,x) = (p{x) for X e R^, 



(3.14) 



(3.15) 



w*{t,x) < liminf u„(t„,x„) + (f>{t,x) — (f>{t,x) 



< w*(t,x) x) - <j){t,x). 



§3 Perron's Method and Existence 



123 



From the above inequalities and the assumption on 4>, we get lim„_>.oo(in, 
{t,x,w* {t,x)) (without passing to a subsequence). Since m„ is a viscosity sub- 
solution of (|3.14l) . by definition we have 

dt(l>{tn,Xn) - G{t (tn,Xn),D(j){tn,Xn),D'^(j){tn,Xn)) < 0. 

Letting n oo, we conclude that 

dt(t){t,x) - G{t,x,w*{t,x),D(f>{t,x),D^(j){t,x)) < 0. 

Thus w* is a viscosity subsolution of p.l4[) by definition. □ 
The second step in the proof of Theorem 13.11 is a simple "bump" construction 
that we now describe. Suppose that u is a viscosity subsolution of p.l4p on 
(0, cx)) X and is not a viscosity supersolution of (13.141) . so that there is 
{t,x) € (0, oo) X and 4> e with u,,{t,x) = 4>{t,x), > on (0,oo) x 
M.^\{t,x) and 

dtcj){t, x) - G{t, X, cj){t, x),D(l)[t, x),D^(f){t, x)) < 0. 

The continuity of G provides r,6i > such that Nr = {{s, y) : — t\ + \y — 
r} is compact and 

dt(j)-Gis,y,(b + S,D(j),D^(^) < 

for all s,y,5 £ Nr x [0,5i]. Lastly, we obtain ^2 > for which > (j) + $2 on 
dNr- Setting 6o — min((5i, (^2) > 0, we define 

jj _ i max(it, <p + Sq) on Nr 

1 u elsewhere. 

By the above inequalities and Lemma l3.2[ it is easy to check that C/ is a viscosity 
subsolution of (|3.14l) on (0,oo) x R^. Obviously, U > u. Finally, observe that 
Ut,{t,x) > max{u^,{t,x),(f>{t,x) + Sq) > u*(i,a;); hence there exists (s, j/) such 
that U{s,y) > u{s,y). We summarize the above discussion as the following 
lemma. 

Lemma 3.3 Let u be a viscosity subsolution of (|3.14p on (0, 00) x R^. // 
fails to be a viscosity supersolution at some point (s, z), then for any small k > 
there is a viscosity subsolution of (|3.14|) on (0, 00) x R-'^ satisfying 

( UK,{t,x) > u{t,x) and sup(C/k — m) > 0, 

\ [/^(t, x) = a;) /or y^T^^sj + l^^^^zp > K. 

Proof of Theorem 13. li With the notation of the theorem observe that < 
<W <W* <u* and, in particular, W^{0,x) = W{0,x) = W*{0,x) = ip{x) 
for X e R^. By lemma IX2l W* is a viscosity subsolution of p.l4p and hence, 
by comparison, W* < u. It then follows from the definition of W that W = W* 
(so is a viscosity subsolution). If W^, fails to be a viscosity supersolution at 
some point {s,z) E (0, 00) x R^, let be provided by Lemma [3.31 Clearly 
u < Wk, and Wk{0,x) = (p{x) for sufficiently small k. By comparison, < u 
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and since W is the maximal viscosity subsolution between u and u, we arrive at 
the contradiction Wk, < W. Hence is a viscosity supersolution of p.l4p and 
then, by comparison for (I3.14p . W* ~ W < W^, showing that W is continuous 
and is a viscosity solution of p.l4p . The proof is complete. □ 
Let G : x E>{N) — >■ IR. be a given continuous sublinear function monotonic in 
A G S{N). We now consider the existence of viscosity solution of the following 
G-equation: 

dtU-G{Du,D^u) ^0, u{0,x) ^ ip{x). (3.16) 

Case 1: If </? £ C^(R^), then u{t, x) = AH + ip{x) and u{t, x) ^ Mt + tp{x) are 
respectively the classical subsolution and supersolution of p.l6p . where M — 
infj-gRN G{Dip{x), D'^ip{x)) and M = sup^g^w G{Dip{x), D'^ip{x)). Obviously, u 
and u satisfy all the conditions in Theorem 13. II By Theorem 12. 6[ we know the 
comparison holds for (j3.16p . Thus by Theorem 13. 11 we obtain that G-equation 
(|3.16p has a viscosity solution. 

Case 2: If e Gf,(R^) with \im^^^^^ = 0, then we can choose a sequence 
^Pn G C'b (K^) which uniformly converge to if. For (y9„, by Case 1, there exists a 
viscosity solution u^'^ . By comparison theorem, it is easy to show that u"^"- is 
uniformly convergent, the limit denoted by u. Similar to the proof of Lemma 
13.21 it is easy to prove that u is a viscosity solution of G-cquation p. 161) with 
initial condition ip. 

Case 3: li Lp € G(R^) with polynomial growth, then we can choose a large 
/ > such that <^(x) = if{x)^~^(x) satisfies the condition in Case 2, where 
^(x) = (1 + It is easy to check that m is a viscosity solution of G- 

equation p.l6p if and only if u{t,x) = u{t,x)£^~^{x) is a viscosity solution of 
the following PDE: 

dtu - G{x, u, Du, D^u) = 0, u(0, x) = <p, (3.17) 
where G{x, v,p, X) — G{p + vi]{x),X +p (X) ri{x) + ri{x) ®p-\- vn{x)). Here 
Tl{x) ■.^r\x)Di{x)=l{l + \x\')-^x, 

n{x) ~ C^{x)D^£,{x) = /(I + l-Tp)-!/ + l{l ~ 2)(1 + \x\'^)-^x ® x. 

Similar to the above discussion, we obtain that there exists a viscosity solution 
of p.l7p with initial condition (p. Thus there exists a viscosity solution of G- 
equation p.l6p . 

We summarize the above discussions as a theorem. 

Theorem 3.4 Let ip G G(R"'^) with polynomial growth. Then there exists a 
viscosity solution of G-equation (I3.16P with initial condition ip. 

Theorem 3.5 Let the function G be given as in the previous theorem and let 
G{t,x,p,X) : [0,oo] X R^ x R^ x S{N) R 6e a given function satisfying 
condition (G) in Theorem \2.S[ We assume that G is dominated by G in the 
sense of 

G{t,x,p,X)-G{t,x,p\X') < G{p~p\X-X'), 

for each t,x, p,p' and X, X'. 
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Then for each given ip G C(M^) with polynomial growth condition the viscosity 
solution of dfU — G(t,x, Du, D^H) = with u\t=o — V exists and is unique. 
Moreover the comparison property also holds. 

Proof. It is easy to check that a G-solution with u\t=o = tf is & G-supersolution. 
Similarly denoting G^,{p,X) := —G{—p, ~X), a G*-solution with u\t=Q = (/s is a 
G-subsolution. 

Wc now prove that comparison holds for G-solutions. Let ui be a G-supersolution 
and U2-be a G-subsolution with ui|t=o = </'i G SC([0;r) x R^), M2|t=o = 
ip2 G PSC([0,T) X K.^), then by the above domination theorem we can ob- 
tain ui — U2 < u where u is a G-supersolution with u\t=o = — <^2- On 
the other hand it is easy to prove that, in the case when ipi < ip2, u — 
{(pi{x) — ip2ix))l[Qj{t) such type of G-supersolution. Consequently we have 
ui < M2- This implies that the comparison holds for G-equations. We then can 
apply Theorem 13. ll to prove that G-solution u with u|t=o = V exists. □ 



§4 Krylov's Regularity Estimate for Parabolic 
PDE 

The proof of our new central limit theorem is based on a powerful C'1+"/2:2+q_ 
regularity estimates for fully n onlin ear parabolic PDE obtained in Krylov [t^. 
A more recent result of Wang [119[_(the version for elliptic PDE was initially 
introduced in Cabre and Caffarelli [17|), using viscosity solution arguments, can 
also be applied. 

For simplicity, we only consider the following type of PDE: 

dtu + G{D^u,Du,u)^0, u{T,x) = ip{x), (4.18) 

where G : S{d) x R'' x R ^ M is a given function and ip e Gb{R'^). 

Following Krylov [zil, we fix constants K>e>0,T>0 and set Q = (0, T) x 

R''. Now we give the definition of g{s, K, Q) and g{e, K, Q). 

The following definition is according to Definition 5.5.1 in Krylov [76| . 

Definition 4.1 Let G : E>{d) x M'^ x R M 6e given, written it as G{uij,Ui,u), 
i,j — l,...,d. We denote G € g{e,K,Q) if G is twice continuously dijferen- 
tiable with respect to {uij,Ui,u) and, for each real-valued Uij — Uji, Uij — Uji, 
Ui, Hi, u, u and A*, the following inequalities hold: 

e\\? <Y.X'XJdu^^G<K\\\\ 

\G - ^ii„-9„.^G| < M<^{u){l + 
|5„G| + {l + Y. K\) ^ |9„.G| < Mf (u)(l -I- \u,\^ + J2 
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i.j i i,j 



where the arguments {uij, Ui, u) of G and its derivatives are omitted, rj = {uij, Ui,it), 
and 

^(riKri) — E ^^j^rsdl^^^^^G + 2^U,jUrdl^^^^G + 2^il,jildl^^^G 

i,j,r,s i-ij 

M^{u) and M^{u,Uk) are some continuous functions which grow with \u\ and 
UkUk and M2 > 1- 

Remark 4.2 Let el < A = (aij) < KI. It is easy to check that 

+ b,Ui + cu 

belongs to Q{e, K, Q) . 



The following definition is Definition 6.1.1 in Krylov 76[. 

Definition 4.3 Let a function G — G{uij,Ui,u) : S(d) x M'' x M — > M be given. 
We write G G Q{e,K,Q) if there exists a sequence Gn G Q{£tKtQ) converging 
to G as n ^ 00 at every point {uij,Ui,u) € S{d) x M'' x M such that 

(i) M^' = Mf ^ =•••=: Aff , i = 1, 2; 

(ii) for each n = 1, 2, . . the function Gn is infinitely differentiate with respect 

to {utj,Ui,u); 

(ill) there exist constants Sq Sq > and AIo —: AIq > such that 

Gn{u^j,0,-MQ) > So, G„(-Uy ,0,Mo) < 

for each n>l and symmetric nonnegative matrices (uij). 

The following theorem is Theorem 6.4.3 in Krylov fld\ , which plays important 
role in our proof of central limit theorem. 

Theorem 4.4 Suppose thatG £ G{£,LC,Q) andip € C6(R'') with sup ^^Tg^a \ip{x)\ < 
. Then PDE (j4.18l) has a solution u possessing the following properties: 

(i) u e G([0,T] X R''), |u| < on Q; 
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(ii) there exists a constant a G (0,1) only depending on d,K,s such that for 
each K > 0, 

II'^IIci+°/2.2+'»([o,t-k]xR'') <oo. (4-19) 

Now we consider the G-equation. Let G :M.'^ x S{d) — > M be a given continuous 
sublinear function monotonic in ^ S §(d). Then there exists a bounded, convex 
and closed subset E C M'' x S+(d) such that 

G{p, A) = sup [l-ti[AB] + {p, q)] for {p, A) eR'' x §(rf). (4.20) 
iq.B)e^ 2 

The G-equation is 

dtU + G{Du,D'^u) = 0, u{T,x) ip{x). (4.21) 

We set 

u{t,x) = e*~'^u(t,x). (4.22) 
It is easy to check that u satisfies the following PDE: 

dtu + GiDu,D^u) -u^O, u{T,x) = ip{x). (4.23) 

Suppose that there exists a constant e > such that for each A,A^ S{d) with 
A> A, we have 

G{0,A)^G{0,A)>etT[A~A]. (4.24) 

Since G is continuous, it is easy to prove that there exists a constant K > 
such that for each A,Ag S{d) with A> A, we have 

G{0,A)-G{0,A)<Ktr[A-A]. (4.25) 

Thus for each (q, B) e S, we have 

2el <B < 2KI. 

By Remark I4.2[ it is easy to check that G{uij,Ui,u) := G{ui,Uij) — u € 
Q{e,K,Q) and S^f — can be any positive constant. By Theorem 14.41 and 
(|4.22|) . we have the following regularity estimate for G-equation (|4.21|) . 

Theorem 4.5 Let G satisfy (|420)) and (|4?24t . ip e Cb{R'^) and let u be a 
solution of G-equation (|4.21[) . Then there exists a constant a G (0,1) only 
depending on d, G, e such that for each k > 0, 

||w||ci+o/2,2+„([o,T-K]xR<i) < oo. (4.26) 
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A Coherent acceptable set [T3] 

Bh{fl) Space of bounded functions on [HI] 

B C?-Brownian motion [551 

(B) Quadratic variation process of B HHl 

Cb(ri) Space of bounded and continuous functions on [HI] 

C6,iip(M") Space of bounded and Lipschitz continuous functions on R" [3] 

C;,iip(R") Space of locally Lipschitz functions [2] 

E Sublinear expectation [T] 

E G-expectation [39] 

H Space of random variables [1] 

Space of all S(r2)-measurable real functions [HI] 
Completion of -Bb(ri) under norm || • ||p [IS] 
Completion of C;,(r2) under norm || • ||p [95] 

Mg'°(0,T) Space of simple processes |42J 

Mg(0, T) Completion of M^'^{0, T) under norm {IE[/„^ \rit\Pdt]Y/P ^ 

M^iO,T) Completion of M^'^{0,T) under norm {J^ E[\Tjt\''dt]y/P M 

q.s. Quasi-surely [M] 

§{d) Space oi d x d symmetric matrices [TH| 

(d) Space of non- negative d x d symmetric matrices [H] 

p Coherent risk measure [T3| 

(fJj'HjE) Sublinear expectation space [D 

= Identically distributed [7| 

(a;, y) Scalar product of x, y e M" 

|x| Euclidean norm of x 

(A, B) Inner product (A, B) := tr[AB] 
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G-normal distribution, 18 
G-submartingale, 71 
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Distribution, 6 
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Extension, 108 
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Prokhorov's criterion. 111 
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